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ABSTRACT 


This  study  is  a  theoretical  investigation  of  the  polarization  properties 
of  the  paraboloidal  antenna  for  both  the  focal  feed  system  and  the  Cassegrain 
feed  system  when  the  feed  horn  is  either  circularly  or  linearly  polarized. 

The  properties  are  analyzed  in  terms  of  general  analytic  equations  and  also 
in  terms  of  quantitative  graphs  pertaining  to  the  specific  case  of  the  30-foot 
diameter  paraboloidal  antenna  which  is  part  of  The  MITRE  Corporation's 
Bedford  Radar  Installation. 
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SECTION  I 


INTRODUCTION  AND  SUMMARY 

SCOPE  OF  THE  STUDY 

The  scope  of  this  study  is  primarily  confined  to  the  theoretical  investi¬ 
gation  of  the  polarization  properties  of  the  paraboloidal  antenna.  Thus, 
answers  are  sought  to  the  questions:  What  is  the  ellipticity  (ratio  of  major  to 
minor  axes)  in  the  near-  and  far-fields  when  the  feed  is  circularly  polarized? 
What  sort  of  polarization  prevails  in  the  near-  and  far-fields  when  the  feed 
is  linearly  polarized? 

The  answers  are  given  in  terms  of  general  analytic  equations  (Section  II, 
for  the  focal  feed  system;  Section  III,  for  the  Cassegrain  system),  and  also 
in  terms  of  quantitative  graphs  pertaining  to  the  MITRE  Antenna  (Section  IV, 
for  the  focal  feed  system;  Section  V,  for  the  Cassegrain  system). 

Field  amplitude  distributions  (as  functions  of  various  parameters)  are 
given  considerable  attention  since  this  information  must  be  obtained  for 
polarization  purposes,  and  since  amplitude  comparisons  with  known  or 
experimental  values  serve  as  checks  on  the  equations  and  computations. 
Effects  of  paraboloid  imperfections  and  of  defocusing  both  in  the  focal  plane 
and  along  the  boresight  axis  are  also  considered. 

SUMMARY  OF  RESULTS 

Some  of  the  more  important  findings  of  the  study  are  as  follows: 

(1)  For  any  circularly  polarized  primary  illumination,  ellipticity  on 
boresight  is  0  db  (perfect  circularity). 
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(2)  Off  boresight,  ellipticity  is  strongly  dependent  on  the  primary 
illumination  distribution  whose  practical  extremes  are  the  ideal 
equibeam  and  the  dipole-like  distributions.  Within  the  3-db  beam- 
width,  the  ellipticity  may  be  practically  0  db  (as  with  the  equibeam 
distribution  in  the  focal  feed  system),  or  as  high  as  0.  8  db  (as  with 
the  dipole-like  distribution  in  the  focal  feed  system).  The 
Cassegrain  feed  system  is  relatively  less  exacting  regarding  the 
type  of  primary  illumination,  and,  within  the  3-db  beamwidth, 
ellipticity  is  about  0. 15  db  for  either  distribution. 

Ellipticity  increases  in  the  region  of  the  null  and  decreases  again 
in  the  region  of  the  sidelobe  peaks. 

The  magnitude  of  the  ellipticity  is  constant  with  azimuth  (at  a 
given  elevation  angle),  provided  that  the  source  can  be  considered 
as  infinitely  small  in  extent  (as  with  the  equibeam  and  dipole-like 
distributions,  but  not  writh  the  pseudo-equibeam  distribution,  which 
represents  another  extreme  case  of  illumination  imperfection). 

(3)  When  linear  polarization  is  applied  at  the  feed,  the  polarization  of 
the  near-  and  far-fields  suffer  a  deviation  angle  from  the  direction 
of  primary  illumination  which  depends  on  azimuth  and  increases 
(or  decreases)  as  would  ellipticity,  were  circular  polarization 
applied. 

(4)  For  the  focal  feed  system,  ellipticity  (at  a  given  elevation  angle) 
is  lower  as  the  main  dish  is  flatter  (longer  focal  length),  or  as 

the  flat  portion  of  the  dish  is  favored  at  the  expense  of  the  periphery 
region  (as  achieved  by  a  large  tapering  of  the  primary  illumination). 
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(5)  For  an  imperfect  primary  distribution,  such  as  the  dipole-like, 
the  Cassegrain  system  is  superior  insofar  as  the  subreflector 
subtends  a  smaller  angle  from  the  source  than  that  subtended  by 
the  main  dish  from  its  virtual  source.  The  overall  improvement 
over  the  focal  feed  system  is  comparable  to  that  afforded  by  the 
effective  "magnification"  of  focal  length  of  ray  optics  .  Another 
viewpoint  (Section  V)  is  that  illumination  defects  at  the  edge  of 
the  subreflector  are  not  magnified  at  the  edge  of  the  main  dish, 
although  their  angular  locations  are  magnified. 

(6)  Near-field  and  far-field  ellipticities  are  almost  equal  in  magnitude 
within  the  half-power  beamwidth,  but  show  large  differences  near 
the  nulls.  The  far-field  can  be  even  more  closely  simulated  at 
close  range,  even  with  regard  to  gain,  by  judicious  de-focusing 
along  the  boresight  axis. 

(7)  Ellipticity  is  hardly  affected  by  random  but  sizeable  (one  tenth  of 
a  wavelength)  imperfections  on  the  surface  of  the  paraboloidal 
main  dish  -  although  gain  (but  not  directivity)  is  lowered.  Semi- 
correlated  imperfections  of  comparable  magnitude  are  far  more 
detrimental. 

(8)  Ray  optics  methods  of  treating  the  Cassegrain  system  are  some¬ 
what  inexact;  yet  they  afford  worthwhile  qualitative  information. 

MECHANISMS  OF  POLARIZATION 

The  polarization  mechanisms  involve,  for  the  focal  feed  system,  the 
interplay  of  the  reflector  paraboloidal  surface  with  the  primary  illumination, 
and,  for  the  Cassegrain  system,  interplay  of  reflecting  surfaces  with 
primary  and  secondary  illuminations. 
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Dipole- Like  Distribution 


Consider  first  the  dipole-like  distribution.  *  This  type  of  illumination 
originates  from  a  single  point  source  (phase  center)  to  which  is  associated  a 
spherical  wavefront.  On  this  spherical  wavefront,  the  H-field  is  oriented 
along  parallel  "latitude"  lines  and  its  strength  varies  with  azimuth  (for  a 
given  elevation  angle)  so  that  its  magnitude  in  the  H-plane  is  different  from 
its  magnitude  in  the  E-plane  (Figure  la). 

On  the  paraboloidal  dish,  the  tangential  components  of  the  parallel 
"latitude"  H-lines  give  rise  to  current  lines,  which  are  no  longer  parallel, 
but  contain  "cross-polarized"  components.  Thus,  an  x-oriented  dipole  gives 
rise  not  only  to  x-oriented  currents  but  also  to  y-currents  which  are 
functions  of  azimuth  and  peak  up  in  the  planes  45  degrees  away  from  the 
x-z  plane.  It  is  those  currents  which  are  mostly  responsible  for  the  cross- 
polarized  components  in  the  near-  and  far-fields  and  for  shifting  the  linear 
polarization  direction.  The  concavity  of  the  paraboloidal  dish  induces 
z-currents  which  also  produce  cross-polarized  components,  but  their  effect 
is  small  compared  to  that  of  the  y-directed  currents.  All  the  extraneous 
currents  (y-  and  z-directed)  produce  a  cross-polarized  pattern  which  does 
not  correspond  to  that  produced  by  the  x-directed  currents.  As  a  result, 
cross-polarized  terms  are  most  effective  in  the  "null"  region  where  the 
in-polarized  components  are  weakest. 


* 

This  distribution  is  considered  in:  E.  U.  Condon,  "Theory  of  Radiation 
from  Paraboloid  Reflectors,  "  Research  Report  SR- 105,  24  Sep  1941; 

S.  Silver,  "Microwave  Antenna  Theory  and  Design,  "  p.  419,  Rad.  Lab. 
Series,  Vol.  12,  1947;  L.  Thourel,  "The  Antenna, "  pp.  269-270,  John 
Wiley  and  Sons,  Inc.  ,  New  York,  1960. 
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Figure  1  .  Distributions  on  Spherical  Wavefront 
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When  the  feed  is  circularly  polarized  (with  dipole-like  distribution),  the 
field  on  the  spherical  wavefront  is  elliptically  polarized,  except  on  the  bore- 
sight  axis.  Although  the  ellipticity  arises  both  from  the  fact  that  the  magnitude 
of  each  quadrature  component  varies  with  azimuth,  and  from  the  fact  that  the 
direction  of  each  quadrature  component  is  along  "latitude"  lines,  the  ellipticity 
on  the  spherical  wavefront  is  constant  in  amplitude  with  azimuth  and,  further¬ 
more,  the  directions  of  the  major  axes  of  the  ellipses  present  cylindrical 
symmetry  with  the  boresight  axis.  These  elliptically  polarized  fields  induce 
elliptically  polarized  currents  on  the  paraboloidal  dish.  These  currents  can 
be  considered  as  the  sums  of  two  counter-rotating  circularly  polarized  cur¬ 
rents  which  give  rise  to  counter-rotating  circularly  polarized  fields.  Again, 
because  of  different  distributions  over  the  dish  for  the  counter-rotating  cur¬ 
rent  components,  the  ratio  of  the  two  counter-rotating  circularly  polarized 
fields  at  the  observation  point  is  not  constant  with  elevation,  but  is  constant 
with  azimuth. 

Pseudo-Equibeam  Distribution 

For  the  pseudo-equibeam  distribution,  the  situation  is  somewhat  dif¬ 
ferent  This  distribution  is  like  the  dipole-like  insofar  as  the  direction  of  the 
II-field  is  still  along  parallel  "latitude"  lines.  However,  the  H-field  strength, 
on  the  spherical  wavefront,  is  made  independent  of  azimuth  (Figure  lb). 

Under  linear  polarization,  the  H-field  still  induces  y-directed  current 
components  which  peak  up  in  the  45-degree  planes  and  give  rise  to  cross- 
polarized  components  in  the  near-  and  far -fields.  The  situation  is  similar 
to  that  of  the  previous  case. 

When  the  source  is  circularly  polarized,  the  field  on  the  wavefront  no 
longer  exhibits  constant  ellipticity  with  azimuth.  There  are  now  "preferred" 
planes  (corresponding  to  the  source  quadrature  component  planes)  where  the 
field  is  circularly  polarized,  whereas  it  is  strongly  elliptically  polarized  in 
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in  planes  45  degrees  away.  As  a  result,  the  near-  and  far-fields  exhibit 

4 

an  ellipticity  which  depends  on  azimuth. 

Equibeam  Distribution 

With  the  equibeam  distribution,  the  field  strength  on  the  spherical 
wavefront  is  constant  with  azimuth,  as  for  pseudo-equibeam  distribution. 
However,  the  field  direction  is  no  longer  along  parallel  "latitude"  lines  but 
is  so  directed  that  two  space  orthogonal  sources  produce  fields  which  are 
orthogonal,  not  only  on  the  boresight  axis,  but  everywhere  on  the  spherical 
wavefront  (Figure  lc).  It  turns  out  that  the  projection  of  this  non-parallel 
field  on  the  paraboloidal  surface  induces  parallel  currents.  Thus,  an 
x-directed  source  produces  no  y-directed  current.  Moreover,  the  x-directed 
current  magnitude  is  independent  of  azimuth.  The  only  mechanism  by  which 
ellipticity  comes  into  play  is  through  the  z-directed  currents  which  give  rise, 
in  the  near-  and  far-fields,  to  x-  and  y-components,  as  the  elevation  angle  is 
increased  from  zero.  Again,  this  effect,  which  is  now  small  compared  to  the 
other  two  cases,  becomes  maximum  in  the  "null"  region  where  the  normal 
components  are  small. 

Cassegrain  System 

In  the  Cassegrain  system,  the  three  distributions  are  almost  the  same 
on  the  subreflector  because  it  subtends  a  much  smaller  angle  from  the  source. 
However,  even  with  the  equibeam  distribution  (whose  field  on  the  spherical 
wavefront  is  circularly  polarized  everwhere),  the  currents  induced  on  the 
hyperboloidal  subreflector  have  rather  small,  but  finite,  cross-polarized 
components.  Through  a  process  related  to  the  effective  magnification  of  the 
main  dish  focal  length  (ray  optics  viewpoint),  the  relative  magnitudes  of  the 
various  components  are  essentially  preserved;  whereas,  the  angular  locations 
are  "magnified"  from  the  subreflector  to  the  main  dish. 
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Distribution  Models 


From  a  practical  standpoint,  the  equibeam  distribution  can  be  considered 
as  the  ideal  distribution  (ideal  feed-horn  compensation).  It  can  be  experi¬ 
mentally  verified  by  applying  circular  polarization  to  the  horn  and  measuring 
the  circularity  of  the  primary  illumination  on  the  spherical  wavefront.  The 
dipole-like  distribution  is  a  fair  approximation  to  the  uncompensated  square 
or  circular  horn;  whereas,  the  pseudo-equibeam  distribution  represents  an 
attempt  to  compensate  the  horn  for  equal  amplitude  in  two  orthogonal  planes 
without  regard  to  cylindrical  symmetry.  Those  are  the  three  arbitrary 
mathematical  models  of  primary  distributions  which  are  considered  in  this 
study. 

ANALYTICAL  APPROACH 

In  Section  II,  the  paraboloidal  dish  is  considered  as  a  scatterer  to  the 
primary  illumination.  The  starting  point  of  the  study  is  the  exact  integral 
equation  for  perfect  conductors  (form  of  Kirchhoff's  Equation)  relating  the 
scattered  H-field  to  the  field  on  the  scatterer.  Simplifying  physical  optics 
assumptions  involving  size,  curvature,  etc.  ,  and  certain  range  assumptions 
are  introduced.  The  various  geometrical  vectors  related  to  the  paraboloid 
are  expressed  in  Cartesian  components,  and  each  component  is  expressed 
in  terms  of  elevation  angle,  azimuth  angle,  and  constants  of  the  system. 

The  geometry  of  the  observation  point  is  also  reduced  to  elevation,  azimuth, 
and  range. 

As  the  "geometry"  is  expressed  analytically,  so  is  the  illumination 
function.  Cartesian  components  in  terms  of  elevation  and  azimuth  are 
obtained  for  the  three  distributions  discussed  in  the  preceding  subsection. 
Those  distributions,  in  turn,  can  be  "tapered"  by  appropriate  amplitude 
functions  independent  of  azimuth.  Those  functions  are  set  forth  arbitrarily 
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or  are  taken  from  measurements.  Circular  polarization  is  introduced 
immediately,  since  it  is  a  simple  matter,  when  considering  linear  polariza¬ 
tion,  to  suppress  one  of  the  quadrature  components. 

After  all  the  elements  have  been  expressed  analytically,  the  various 
field  components  can  be  separated  in  terms  of  space  directions  and  time 
phases.  Concurrently,  some  "tools"  are  developed  to  compute  ellipticity 
magnitude  and  orientation  in  terms  of  the  various  field  components  and, 
when  necessary,  to  change  the  coordinate  system  for  a  more  meaningful 
"view"  of  the  field. 

In  Section  III,  the  Cassegrain  system  is  analyzed  by  considering  the 
field  at  the  paraboloidal  dish  as  that  being  scattered  by  the  subreflector. 

The  "geometry"  of  the  hyperboloidal  subreflector  must  then  be  treated 
analytically,  although  the  primary  illumination  functions  are  the  same  as 
those  considered  in  Section  II.  After  it  has  been  solved,  the  field  on  the 
paraboloid  is  treated  as  the  primary  illumination  of  Section  II. 

The  geometrical  problems,  primary  distributions,  and  other  "tools" 
transcending  the  limited  scope  of  this  study  are  relegated  to  the  Appendices. 


COMPUTER  PROGRAMS* 

A  number  of  computer  programs  are  developed  to  perform  the  calcu¬ 
lations  described  in  this  study.  Listings  of  some  typical  programs  are 
reproduced  in  Appendix  K. 

The  programs,  written  in  FORTRAN  IV,  are  run  on  the  IBM  7030 
(STRETCH)  Computer  of  MITRE's  SDL  facility. 


* 

Written  by  J.  L.  Pearlman 
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Since  several  of  the  programs  are  time  consuming,  a  deliberate  effort 
is  made  to  eliminate,  wherever  possible,  the  recalculation  of  often  used 
intermediate  results.  This  is  accomplished  by  devoting  large  amounts  of 
memory  to  the  storage  of  these  results.  In  certain  cases,  the  programs 
are  written  in  two  parts,  with  the  output  of  the  first  part  being  stored  on 
tape  to  be  used  as  input  to  the  second  part,  as  in  the  case  of  the  Cassegrain 
system.  In  this  case,  the  field  scattered  onto  the  parabolic  reflector  from 
the  hyperbolic  subreflector  is  computed  and  stored  on  tape  for  further 
processing  by  programs  which  compute  the  field  components  at  various 
points. 
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SECTION  II 


ANALYSIS  OF  FOCAL  FEED  SYSTEM 


SYMBOLIC  FIELD  EQUATIONS  AND  APPROXIMATIONS 

A  paraboloidal  dish  illuminated  from  a  source  located  at  its  focus  acts 
as  a  scatterer  to  the  primary  illumination  (Figure  2).  and  the  exact  equation 


for  the  scattered  H-field  at  the  observation  point  P,  Hg  i  .  is: 


(1) 


where 

H  =  total  11-field  at  point  Q  on  the  surface  of  the  paraboloid, 
n  =  outward  normal  unit  vector  at  point  Q. 
k  =  propagation  constant,  equal  to  2  n/\  , 
r  =  distance  from  the  observation  point  P  to  point  Q. 
dA  =  differential  area  element  at  point  Q.  and 
S  =  surface  of  the  paraboloid. 

An  exact  expression  for  \  q  (e— r/r)  is: 


It  is  seen  that  when  r  is  greater  than  15  wavelengths,  the  magnitude  of  the 
l/r  term  is  but  one  percent  the  magnitude  of  k.  Furthermore,  if  the 


11 


observation  point  P  is  at  least  15  focal  lengths  from  the  apex  of  the  parab¬ 
oloid,  then,  as  far  as  amplitude  is  concerned,  r  can  safely  be  replaced  by 
the  constant  factor  r^,  the  distance  from  the  observation  point  P  to  the 
focus  (coordinate  center)  of  the  paraboloid.  An  approximate  expression  for 
Equation  (2),  valid  everywhere  but  in  the  very  near  field,  becomes 


VQ 


e 


— jkr 


ar 


(3) 


where  a^  is  the  unit  vector  directed  from  point  P  to  the  running  point  Q. 

Since  it  can  be  assumed  that  the  paraboloid  is  large  compared  to  wave¬ 
length,  and  smooth,  the  approximation 

n  x  H  =  2  (n  x  H.)  (4) 
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essentially  prevails,  where  H.  is  the  magnetic  field  incident  on  the  parab¬ 
oloid  at  point  Q.  The  phase  of  H.  can  be  explicitly  brought  out  according  to 


H. 

1 


H. 
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e 


-jkp 


(5) 


where  p  is  the  distance  from  the  origin-phase  center  to  the  running  point  Q. 
The  magnitude  and  orientation  of  H.  wall  be  subsequently  restricted.  Finally, 
an  implication  of  Equation  (4)  is  that  S  is  to  be  taken  as  only  that  portion 
which  is  directly  illuminated. 

With  all  those  approximations,  Equation  (1)  becomes 

V  ■  2^77  //«*  Hio>  ,  a/e-i^Od  A  (6, 

where  the  quantities  n,  H  ,  p  ,  and  dA  are  functions  of  the  running 
point  Q  only,  where  r  is  a  function  of  the  observation  point  P  only,  and 
where  a^  and  r  are  functions  of  both  P  and  Q. 

The  assumptions  on  which  the  approximations  of  Equation  (6)  are  based 
can  be  summarized  as  follows: 

(1)  Dish  dimensions  are  large  compared  to  wavelength, 

(2)  Observation  point  is  at  least  15  focal  lengths  away, 

(3)  System  is  fed  from  a  phase  center  located  at  the  focus. 

The  next  task  is  obtaining  analytical  expressions  for  the  various  quanti¬ 
ties  in  Equation  (6)  in  terms  of  the  coordinates  of  point  Q  ,  defined  by  the 
elevation  angle  p  and  the  azimuthal  angle  ^  ,  and  of  point  P  ,  defined  by 
r  i  and  by  the  elevation  angle  0  and  the  azimuthal  angle  <j>  . 


GEOMETRY  OF  THE  PROBLEM 

The  "geometry"  of  the  paraboloid  is  illustrated  in  Figure  3  (cross- 
section  through  the  x-z  plane)  and  Figure  4  (coordinates  of  running  point  Q 
on  the  paraboloid).  Briefly,  Q  is  a  general  point  of  coordinates  x,  y,  z. 
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Figure  3.  Geometry  of  the  Paraboloid:  Cross-Section  through  x-z  Plane 


i 

< 


Figure  4.  Geometry  of  the  Paraboloid:  Point  Qq  Rotated  through  (ft  to  Q 
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in  the  x-z  plane,  through  an 


It  can  be  considered  as  the  rotation  of  point  , 
angle  0  .  The  distance  between  the  origin  (focus)  and  is  p  .  The  focal 
length  is  f  .  In  the  final  analysis,  Q  (and  all  the  quantities  dependent  on  it) 
will  be  determined  by  77  .  the  angle  between  the  z-axis  and  p  (the  vector 
from  O  and  Q)  ,  and  by  0  ,  the  angle  between  the  x-axis  and  the  projection 
of  p  on  the  x-y  plane. 

It  is  shown,  in  Appendix  A,  that  because  of  the  geometric  properties  of  the 
the  paraboloid,  the  Cartesian  coordinates  of  n  are 


nx 

—  cos  i/j  cos  77  2 

A 

n 

ny 

= 

—  sin  0  cos  77/2 

nz 

sin  77/2 

In  Appendix  B,  it  is  shown  that  the  expression  for  p  ,  the  vector  from 
the  origin  to  point  Q  ,  is 

p  - » -  (x  sin  77  cos  0  +  y  sin  7 7  sin  0  +  z  cos  77)  (8) 

sin  ^  77/2 

According  to  Appendix  C,  the  expression  for  the  differential  area  dA 
assumes  the  form 

dA  -  f2„sin  2  d  77  d  0  (9) 

sin  }  r//2 

T 

In  order  to  evaluate  the  unit  vector  a^  (Figure  1)  directed  from  the 
observation  point.  P,  to  the  running  point,  Q  ,  and  the  quantity  r  (the 

A 

distance  from  P  to  Q  ),  introduce  first  the  unit  vector  a^  .  directed  from 
the  observation  point  P  to  the  origin  (focus).  The  unit  vector  can  be 
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expressed  in  terms  of  6  (its  angle  with  the  negative  z-axis)  and  of  0  (the 
angle  of  its  projection  in  the  x-y  plane  with  the  x-axis)  as: 


al 


sin  0  cos  0 
sin  6  sin  0 
cos  0 


(10) 


The  reason  for  taking  the  origin  of  0  as  the  negative  z-axis  is  only  for  con¬ 
venience  when  referring  to  measurements. 

Now  r,  being  the  magnitude  of  (a^r  +  *  is,  by  Equations  (8)  and  (10): 


r  magnitude  of 


- ^ - sin  7]  cos  0  —  ri  sin  0  cos  0 

sin~  77/2 

- -  sin  7]  sin  0  —  ri  sin  9  sin  0 

sin  ~  77/2 


f 


sin^  77/2 


cos  77  —  r  j  cos 


0 


which  is  evaluated  as 


f 

V 


2  f2  2frl  ,  .  .  ,  J 

r  =  \  rj  +  - 7 -  —  n  [sin  T)  sin  6  cos  (i Jj  —  cf>)  +  cos  6  cos  77 1 

sin4 17/2  sin^  77/ 2  | 


1/2 


(11) 


The  exact  expression  for  the  vector  is  (a  r^  +  P  )  /  r-  However, 
in  this  expression  which  enters  Equation  (6)  as  a  magnitude  (and  not  as  a  phase), 
the  denominator,  r  ,  can  be  approximated  as  r  ,  provided  that  the  observa¬ 
tion  point  is  at  least  10  focal  lengths  away  from  the  antenna.  The  approximate 
expression  for  a^  becomes(a^  +  p/r^,  that  is 

sin  77  cos  0  —  sin  0  cos  0 


sin ^77/2 


- r -  sin  77  sin  0  —  sin  0  sin  0 

1*1  sinZ  77/2 


f 
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*2  Sln  77/ 2 


cos  77  —  cos  0 
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Except  if  0  is  unduly  large  (larger  than  5  degrees),  the  z-component  ot  a 


can  be  approximated  as  -1  .  Thus,  the  final  expression  for  a  becomes 

2f 

—  cot  n  2  cos  \b  —  sin  0  cos  d> 
rl 


al 


2  f 


cot  77/ 2  sin  0  —  sin  0  sin  c 6 


-  1 


(12) 


This  expression  is  essentially  valid  everywhere  for  amplitude  purposes,  ex¬ 
cept  in  the  very  near  field  (r  less  than  10  focal  lengths). 


PRIMARY  PATTERNS 


The  primary  patterns  are  discussed  in  Appendix  D.  Any  particular 
primary  pattern  is  strictly  a  function  of  the  feed  system,  and  it  can  assume  an 
infinity  of  different  forms.  One  can  consider  the  equibeam  distribution  dis¬ 
tinguished  by  the  fact  that,  when  the  field  is  linearly  polarized,  the  strength 
of  the  field  on  the  spherical  wavefront  has  cylindrical  symmetry.  The 
orientation  of  the  field  on  the  spherical  wavefront  must  obey  clearly  defined 
relationships.  Under  conditions  of  feed  circular  polarization,  the  expression 
for  H  (which  is  circularly  polarized  at  each  point  of  the  spherical  wave- 
front)  is 

2  o 

(—  cosz  0  cos  7]  +  sin"  0)  cos  cot 

+  (—  cos  0  sin  0  [1  f  cost/])  sin  cot 


H jo  =  h  (r/)  sin"  17/ 2 


(—  cos  0  sin  0  [l  +  cos  77])  cos  cot 

2  2  >\ 

+  (—  sinz  0  cos  7]  +  cos  0)  sin  cot 


(13) 


sin  77  cos  ip  cos  cot  +  sin  77  sin  0  sin  cot 


where  sin^77/2  reflects  the  dependence  of  the  field  strength  on  the  distance 
P  .  h(  7/  )  is  an  amplitude  function  independent  of  azimuth,  and  coscot  and 
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sin  iDt  are  temporary  labels  designed  to  tag  the  "vertical"  and  "horizontal" 
components,  respectively. 


It  has  been  found  experimentally  that 

h  (77)  COS1-5  (77  —  T)  ) 


(14) 


fits  the  polar  amplitude  pattern  very  closely  for  an  equibeam  feed,  with  a 
10  db  taper  at  the  edge  of  the  dish  (rio  =  118  degrees). 

The  dipole-like  distribution  is  characterized  by  the  fact  that  when  the 
feed  is  linearly  polarized,  the  orientation  of  the  H-field  is  along  "latitude" 
lines  on  the  spherical  wavefront.  The  field  strength  depends  on  the  azimuth 
angle.  Under  conditions  of  feed  circular  polarization,  the  expression  for 
H  (circularly  polarized  on-axis  only,  but  of  ellipticity  constant  with 
azimuth  angle)  is 


h  (77)  sin“  77/2 


—  COS  T]  COS  CO  t 

—  cos  77  sin  co  t 

sin  77  cos  0  cos  co  t  +  sin  77  sin  0  sin  cot 


(15) 


The  pseudo-equibeam  distribution  is  derived  from  the  dipole-like 
distribution.  The  field  strength  has  been  made  independent  of  azimuth  angle. 
However,  the  H-field  still  lies  along  ’  latitude  ”  lines  on  the  spherical  wave- 
front.  Under  conditions  of  feed  circular  polarization,  the  expression  for  H 
(whose  ellipticity  depends  on  the  azimuth  angle  and  which  is  circularly  polarized 
on  the  ’’vertical"  and  "horizontal”  planes  only)  is 

H|o  h  (77)  sin“  77/2 


—  cos  77  yc  cos  o  t 

—  cos  77  yg  sin  co  t 

(sin  77  cos  0y  )  cos  co  t  +  (sin  77  sin  0  y  )  sin  cot 


(16) 
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where 


(17) 


/  2  -2  2  ,  \~ 1  2 
yc  (cos  r)  +  sin  77  cos  ly ) 

,  2  -2  .  2  m-1/2 

ys  -  (cos  77  +  sin  77  sin  y/ ) 

Linear  polarization  can  be  studied  by  suppressing  either  the  cos  wt 
terms  or  the  sin  at  terms.  Other  non-equibeam  distributions  can  be  approx¬ 
imated  by  appropriate  modifications  of  yc  and  ys  . 


FIELD  FORMULATION  FOR  THE  EQUTBEAM  DISTRIBUTION 


Analytic  expressions  for  each  quantity  in  the  fundamental  equation 
(Equation  6)  are  now  available.  Before  evaluating  the  integral,  it  remains  to 

_ i  Wj*  -f-  p  \ 

calculate  the  exponential  factor,  e  J  v  ,  and  the  cross-product 

[(n  x  Hj0)  x  ']  . 

By  Equations  (11)  and  (173)  of  Appendix  B, 

c-jk  (r+p) 


c  xp .  jk 


1 


f- 


rl2+  ;  4 


2  f 


rl 


1  1/2 


sin  *  7]  2  sin""  rj/2 


[sin  77  sin  0  cos  (0  —  <f>)  4*  cos  6  cos  77  ] 


(18) 

si  11^77/ 2  ( 


By  Equations  (7)  and  (13),  (n  x  H.  )  is  the  vector  whose  components  are 
the  cofactors  of  the  determinant  of 


—  cos  i/j  cos  77  2 

c\  cy 

(—  cos^0  cos  rj  +  sin^0)  cos  cot 
+  f—  cos  0  sin  0  [1  +  cos  77])  sin  cot 


—  sin  0  cos  77/2 

(—  cos  1/7  sin  0  [cos  77  +  ll)  cos  cot 
+  (—  sin^i/f  cos  774-  cos^  0)  sin  cot 


sin  77/2 

(sin  77  cos  0  cos  cot) 

4  (sin  77  sin  0  sin  cot) 


times  h  (77)  sin^7//2  . 


This  expression  simplifies  considerably  to  become 


(n  x  Hj0)  =  h  (77)  sin^  77/2 


—  sin  77/2  sin  cot 
sin  77/2  cos  cot 

sin  0  cos  77/2  cos  cot  —  cos  0  cos  rj/2  sin  cot 


(19) 
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or,  symbolically, 


(n  x  H* 


A  sin  cot 
—  A  cos  cot 
E  sin  cot  —  F  cos  cot 


(20) 


where 


A  =  (— )  h  (77)  siri  77/2  sin  rj/2  A' 

9 

E  =  (— )  h  (77)  sinz  t;/2  cos  rf/2  cos  0  E'  cos  0 

o 

F  =  (— )  h  (77)  sin z  77/2  cos  77/ 2  sin  0  =  E'  sin  0 

It  is  remarkable  that  no  cross-polarized  term  is  present,  that  is, 
vertical  polarization  in  the  feed  (cos  cot  terms  )  gives  rise  to  x-  and  z- 
directed  currents  only  (no  y-component  of  current),  whereas  horizontal 
linear  polarization  gives  rise  to  the  no  x-component  of  current. 

The  cross-product  of  (n  x  H  )  and  a^  (Equation  12)  is  the  vector 
whose  components  are  the  co-factors  of 


A  sin  cot 


—  A  coscot 


~  2f 

—  cot  77/2  cos  0  —  sin#  cos  0 

~  2f 

—  cot  77/2sin  0  —  sin  # sin 0 

/I 

rl 

f  E  sin  cot  —  F  coscot] 
-  1 


that  is,  l(n  x  Hio)  x  a|  /  ] 

/  2f  \ 

A  cos  cot  —  ( -  cot  77/2  sin  0  —  sin  6  sin  0  )  (E  sin  cot  —  f  cos  cot) 

Vri  / 

(2f  \ 

—  cot  77/2  cos  0  —  sin  0  cos  0  ]  (E  sin  cot  —  F  cos  cot) 

rl  / 

(2f  \  /  2f 

—  cot  77/2  sin  0  —  sin  0  sin  6)  +  A  cos  cot  (  —  cot  77  2  cos  0 

ri  /  \ri 


sin  0  cos  0^ 

(21) 
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Consider  now  the  x-component  of  Hgjp  .  The  "in-phase"  terms  are 


r  r  7 

I  j  A  '  exp  H  g  (77)  drj  d<// 


XI 


\3  = 


(22) 


0  Vo 

,2n 


f  f  77  /  2f  .  \ 

1  j  E '  exp  i !  g  (77)  sin  0  ( - cot  77/2  sin  t/<  —  sin  0  sin  0  1  dr;  d 

°  %  '  ,n 


0 

(23) 


The  phase-quadrature  term  arising  from  the  "’horizontal’’  components  (through 
the  z-directed  current)  in  the  primary  pattern  is 

V")  —  J E'  exp  H  g  (r/)  cos  0  cot  *7/2  sin  0  —  sin  0  sin  0^  dr]  d0  (24) 

For  the  y-components,  the  ’’in-phase”  terms  are 


VI  -  XI 


(25) 


\  3 


^ J” F/ exp  |1  g  (77)  cos  0  cot  T)/ 2  cos  0  —  sin  0  cos  0^  dr/  d0  (26) 


The  phase-quadrature  term  arising  from  the  ’’vertical”  components  in  the 
primary  pattern  is 


Y5 


E  '  exp  j !  g  (77)  sin  0 


( 


-  cot  77  2  cos  0  —  sin  0  cos  0 

rl 


d0 


where 

exp  ||  is  given  by  Equation  (18), 

/i  \  1  sin  77  sin  77 

g  (h)  - ^ or  - — - —  (Equation  (9))  , 

sin°77/2  sin0  r//2 

A  '  -  —  [h  (77)  sin"  77  21  sin  rj/2 

E'  —  [h  (77)  sin"  77/2]  cos  77/ 2 

In  the  case  of  small  values  of  |0|  ,  the  z-component  is  neglected, 
although  it  could  be  obtained  from  Equation  (21). 
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FIELD  FORMULATION  FOR  DIPOLE- LIKE  DISTRIBUTIONS 


The  dipole-like  distribution  (Equation  15)  and  the  pseudo-equibeam 
distribution  (Equation  16)  will  be  carried  together  for  the  sake  of  convenience. 
It  will  be  understood  that 

yc  =  ys  =  1 

for  the  dipole  distribution,  and  that 

9  cy  n  1/2 

yc  =  (cos"  77  p  sirr  7/cos^0) 

o  9  9  1/2 

)'  -  (cos“  r)  +  sin*'  T]  sin“  <Jj) 


for  the  pseudo-equibeam  distribution,  as  given  by  Equation  (17). 


The  same  exponential  factor  (Equation  18)  still  applies  as  in  the  pre- 

A  - ► 

ceding  section.  By  Equations  (7)  and  (16),  (n  x  H  )  is  the  vector  whose 
components  are  the  cofactors  of  the  determinant  of 


h  (7/)  sin "  7//2 


X  y 

(—  cos  0  cos  rj/2)  (—  sin  0  cos  77  2) 


(sin  77  2) 


(—  cos  77  yQ  cos  cot)  (  —  cos  77  yg  sin  cot)  (sin  77  cos  0  yQ  cos  cot 

sin  77  sin  0  y  sin  6Jt) 

that  is,  (n  x  Hio) 

sin  77  cos  77,  2  sin  0  (cos  0 yc  coscot  -4  sin  0  yg  sin  cot)  4  sin  77  2  cos  77  y  sin  cot 
H  («?)  sin^/2  sin  77  cos  7/  2  cos  0  (cos  0  yc  cos  cot  -F  sin  0  ys  sin  <ot)  —  sin  77  2  cos  77  yg  cos  cot 

cos  77  cos  77/2  (cos  0  yg  sin  cot  —  sin  0  y  cos  cot) 

or  symbolically, 

(A  +  B)  yg  sin  cot  +  Cyc  cos  cot 

(—  A  +  D)  yc  cos  cot  —  Cyg  sin  cot  (27) 

Eyg  sin  cot  —  Fyc  cos  cot 


("  X  Hio) 
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where 


9 

A  =  h  (77)  sin"-  17/2  sin  r}/2  cos  77  A' 

B  —  h  (77)  sin"- 77/2  sin  77  cos  77/2  sin"- i/f  B'sin^i// 

C  =  —  h  (77)  sin^  77/2  sin  77  cos  77/2  sin  ^  cos  <A  B '  sin  iJj  cos  i/j 

D  =  h  (77)  sin  “  77  2  sin  77  cos  77/2  cos^  i/j  —  B  'cos"  0 

2 

E  h  (77)  sinz  77/2  cos  77  cos  77/2  cos  0  E]'cos  \jj 

o 

1"  h  (77)  sin^  17  2  cos  17  cos  tj/2  sin  (//  -  E '  sin  1 Jj 

Note  that  the  "vertical"  polarization  components  (cos  <ot  terms)  give 
rise  not  only  to  x-  and  z-directed  currents,  but  also  to  y-directed  currents 
(cross-polarized  terms).  In  the  same  manner,  the  "horizontal"  polarized 
components  (sinwt  terms)  give  rise  not  only  to  y-  and  z-cur rents,  but  also 
to  x-currents.  This  effect  does  not  arise  in  the  true  equibeam  distribution 
(see  preceding  sub-section).  It  will  be  appreciated  later  on  how  those  cross- 
polarized  terms  give  rise  to  sizeable  ellipticities  in  the  near-  and  far-fields. 

The  cross-product  of  (n  x  H.q)  and  a^  Equation  (12)  is  the  vector 
whose  components  are  the  co-factors  of 


A 

y 


l(A+  B)  yg  sin  cot  +  Cyc  coscot]  [(—  A  +  D)  yc  coscot  -  Cyg  sin  cot]  [Eyg  sin  cot  -  F  yc  coscot] 
cot  77/2  cos  0  -  sin  6  cos  6  J  cot  77/2  sin  \b  -  sin  6  sin  <i>J  -1 

that  is,  [(n  x  Hjo)  x  a^]  = 

/  2f  \ 

(A-  D)yc  coswt  +  Cygsincot  -  cot  77/2  sin  i>  -  sin  d  sin  d>J  (E  yg  sin  cot  -  F  yc  cos  wt) 

/  2f  v 

(A+  B)  yg  sin  cot  +  Cyc  coscot  +  cot77/2cosi^  -  sindcos^j  (Eyg  sin  cot  -  Fyc  coscot) 

/  2£  \ 

[  (A +B)ys  sin  <ut  +  Cyc  coscot]  coti7/2sin  xb  -  sin  0sin  +  [(A  -  D)  yc  coscot  +  C  yg  sin  rot]  ■ 

( 2f  \ 

—  cot  77  2  cos  \Jj  —  sin  0  cos  c pj 
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Consider  now  the  x-component  of  Hi.  The  "in-phase"  terms  are 

O  ' 


\1  =  +  I  f  A  '  yc  exp  U  g  (tj)  dr)  dip 
°  % 


(29) 

(30) 


X2  =  +  JJ  B'yc  exp  {|  g  (tj)  cos"  ijj  dr)  di ip 

X3  =  +  j'J'  E'yc  exp  {{  g  (77)  sin  1 p  ^  cot  tj/2  sin  ip  —  sin  0  sin  (p^j  dr)  dip  (31) 

The  phase-quadrature  terms  arising  from  the  "horizontal"  components  in  the 
primary  pattern  are 

X4  =  +  JJ B'ys  exp  j  |  g  (77)  sin  ip  cos  iP  dr)  di/i  (32) 

X5  =  —  JJ  E  '  yg  exp  {{  g  (tj)  cos  ip  cot  tj/2  sin  ip  —  sin  0  sin  (6^  drj  d if>  (33) 

For  the  y-components.  the  "in-phase"  terms  are 


Y1 


if- 


ys  exP  U  g  (tj)  dTj  di P 


O 

ys  exp  j(  g  (tj)  sin“  1 p  drj  dip 


(34) 


(35) 


V2  ■  *//,r 

Y3  + ^ ^  E  ys  6X^  ^  ^  (tj)  cos  ^  cot  cos  ^  —  s*n  ^  cos  dr)  dip  (36) 


The  phase-quadrature  terms  arising  from  the  "vertical"  components  in  the 
primary  pattern  are 
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(37) 


^  4  = 


’iS  = 


y  exp  U  g  (77)  sin  0  cos  0  dr]  d0 


/>■ 


-//  E '  yc  exp  H  g  ( V )  sin  i !> 


cot  rj/2  cos  0  —  sin  6 


0  ^  dr;  d0 


where 


exp  1!  is  given  by  Equation  (18), 


/,  x  I  sin  77  sin  7}  ,  ^  /riXv 

g  (h)  - r -  or  - ■= — ‘ —  (see  Equation  (9)), 

sin'1  r\! 2  sin0  rj/2 

A  '  =  +  [h  (77)  sin~  r;/2]  sin  7//2  cos  rj 

B'  -  —  lh  (77)  sin~  r;/2]  sin  77  cos  rj/2 

cy 

E '  =r  +  lh  (77)  sinz  57/2]  cos  77  cos  77  2 

In  the  case  of  small  values  of  6  ,  the  z-component  is  neglected, 

although  it  could  also  be  obtained  from  Equation  (28). 


FAR- FIELD  APPROXIMATIONS 

When  r f  .  then,  from  Figure  1, 


(38) 
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Those  approximations  are  valid  when  r  and  a^  enter  amplitude 
factors,  but  are  still  not  valid  when  they  enter  the  exponential  term  which, 
because  of  its  modulo  2n  characteristic,  is  very  sensitive  to  minor  differ¬ 
ences,  even  if  those  differences  are  minute  compared  to  the  total  argument. 

First,  consider  the  expression  for  r  (Equation  11),  and  re-write  it  as 

.  1 


2 


r  =  rl  V  +  2  ,  4  ~c 


2  f 


o 

r^“  sin*  r)/2  sin“  17/ 2 

We  can  expand  it  as  a  power  series  of  the  form 


‘[sin  17  sin  0  cos  (i Jj  —  <i>)  +  cos  0  cos  //] 


(1  +  X)  2  1  +  7T  - 


\  X2  3  X3 


2  4  8 

and  retain  but  the  first  two  terms,  that  is 


-  rl  I1 


f" 


f 


)  2rj“  s'n4  rj/2  rj  sin“  r;/2 


[sin  rj  sin  0  cos  (i/»  —  q!>)  -t  cos  0  cos  //]  | 


Furthermore,  we  are  assured  that,  for  r  large  enough. 


Fl  y  2 r  j2  sin4  r//2  J 

Consequently,  the  expression  for  r  in  the  far-field  becomes 


r  =  r,  — 


sin2  77  2 


[sin  77  sin  0  cos  (iJj  —  <£)  4-  cos  0  cos  77] 


(39) 


The  expression  for  the  exponential  term  (Equation  18)  in  the  far-field 
becomes 


a-jk(r+p) 


-  exp.  -  jk  rn  - 


f 


*  2  ,n 

sin  7// 2 


[sin  77  sin  6  cos  {iJj  —  <£)  +  cos  0  cos  77  —  1] 


=  exp.  —  jk  •  exp.  j  2k  f 


sin  6  cot  77/2  cos  {ifj  —  <f>) 


(40) 

COS  0  COS  77  —  1 

2  sin"  77/2 
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The  term  e  is  independent  of  the  paraboloid  coordinates  and  can 

be  placed  outside  the  integral  (or  neglected). 

A  further  simplification  results  if  8  is  small  (point  P  near  boresight). 
When  this  is  the  case,  then  Equation  (40)  becomes 

exp.  j 2k f  [sin  0  cot  77/ 2  cos  (0  —  0)  —  1  ] 
or 


exp.  —  j  k  (r  +  p)  -  exp.  j 2 k f  sin  6  cot  77  2  cos  (0  —  0)  (41) 

where  again  the  term  e  J~  has  been  placed  outside  the  integral  (or  neglected). 
From  Equation  (176)  of  Appendix  B,  the  exponent  can  be  expressed  in  terms 
of  the  aperture  D(Q)  of  the  dish  at  point  Q  according  to 

exp.  —  j  k  (r  +  p)  =  exp.  j  - sin  0  cos  (ifi  —  <f>)  (42) 

Thus,  for  small  values  of  8  ,  the  exponent  is  independent  of  f  ,  the 

focal  length,  and  depends  only  on  D(Q)  and  on  8  . 

Aside  from  the  simplification  of  exp.  |j  ,  one  can  discard  the  r  - 
dependent  components  of  X3,  X5,  and  Y3,  Y5  in  the  equibeam  or  dipole-like 
distributions. 

If,  furthermore,  the  study  is  confined  to  small  values  of  8  (6  <  5  degrees), 
the  new  terms  X3,  X5,  and  Y3.  Y5,  which  are  all  products  of  sin  0  ,  will 
be  negligible  with  respect  to  the  dominant  terms  XI  and  Y1  in  the  equibeam 
distribution.  For  the  dipole-like  distribution,  those  terms  will  be  smaller 
than  the  "correction"  terms  X2,  X4,  Y2,  Y4  which  arise  from  the  quadrature 
currents  by  one  (or  more)  orders  of  magnitude. 
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Thus,  under  conditions  of  far-field  and  small  elevation  angles,  the  com¬ 
ponents  of  Hg  j  p  are  as  follows: 

For  equibeam  distribution: 

In-phase  x-component: 


\1 


|  j  cot  t// 2  dr,  di p 


(43) 


In-phase  y-component: 


'i  1 


1  i  cot  r]  '2  dr/  di/j 


(44) 


For  dipole-like  distributions: 
In-phase  x-components: 


XI 


h  (1?)  yc  exp 


H  sin  17/ 2  cos  r/  - = -  dr;  di/> 

si  11  77/  2 


(45) 


X2 


-If 


h  (77)  y  exp  {j  sin  77  cos  77/2  cos~  ifj 


sin  rj 
sin^  77  2 


dr]  d\J; 


(46) 


Quadrature-phase  x-component: 


\4  = 


-// 


h  (7/)  y  exp 


11  sin  77  cos  77/2  sin  iJj  cos  ifj 


dr,  d0 

sin ^  77  2 


(47) 


In-phase  y-components: 


Y 1 


exp 


1 1  sin  77/2  cos  77 


■  SlnT?  ■■  dr,  d0 

•  3  o 

sin  77/  2 


(48) 
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Y2  = 


(49) 


-If 


2 

h  (77)  y  exp  \  \  sin  77  cos  77/2  sinz  0  d 77  d ib 


Quadrature-phase  y-component: 


Y  4  =  —  /  /  h  (77)  y  exp  ||  sin  77  cos  77/2  sin  0  cos  0  - — —  dr/  d«/y  (50) 

sin  77  2 


where  exp  ii  stands  for  Equation  (41). 


FIELD  AMPLITUDE  AND  POLARIZATION 

Consider  the  case  of  small  values  of  0  ,  so  that  the  x  and  y 

coordinates,  at  the  observation  point  P  ,  are  essentially  perpendicular  to 
the  line  of  sight.  (A  coordinate  transformation,  suitable  for  any  value  of 
0%  is  suggested  in  Appendix  E). 

It  is  now*  time  to  substitute,  for  the  real  time  functions  cos  cut  and 
sin  <ot  ,  appropriate  phasors.- 

cos  cot  — ►  1 

sin  cot  — j 


For  a  circularly  polarized  feed  system,  the  total  component  of  II 
is  then 


s|p 


X  =  (XI  +  X2  +  X3)  +  j  (X4  +  X5) 


(51) 


and  the  total  y-component 

'i  =  j  (Y1  +  Y2  +  Y3)  +  (Y4  +  Y5)  (52) 

where  some  of  the  terms  may  be  zero,  depending  on  the  distribution  and 
the  conditions. 
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By  the  considerations  of  Appendix  F,  the  ellipticity  is  obtained  as 


follows:  the  magnitude  of  thr  right-hand  wave  is 

IS*- 1  =  j  [(Re  \  -  ImY)2  +  (1m  X  +  Re  Y)2] 
and  of  the  left-hand  wave 

|S“|  J  [(Re  \  +  ImY)2  +  dm  X  -  Re  Y)2  ] 
The  "voltage"  ellipticity  at  point  P  is  then 


Ellipticity 


1^1  4  |S~| 
I  SI  -  |S-| 


or,  in  db, 


(53) 


(54) 


(55) 


Ellipticity  (db) 


20  log^Q 


M  +  |S~| 

\&\  -  |S"I 


(56) 


If  one  linear  component  in  the  primary  pattern  is  suppressed,  so  that 
the  x-component  only  (cos  cot  )  is  applied,  then  the  x-  and  y-components 
of  H 

S !  P  become  respectively 


X'  =  Xj  +  X2  +  X3  (57) 

Y'  =  V4  +  Y5  (58) 


The  overall  amplitude  is  then 


Amplitude  \J |  X  '  |  2  +  '  |  2 


(59) 


and  the  deviation  angle  from  the  primary  polarization  to  the  polarization 
at  point  P  is 


a  — 


tan 


-i 


XI 

X' 


(60) 
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Equations  (59),  (56).  and  (60)  have  been  computed  and  plotted  under 
various  conditions,  and  the  results  are  found  in  Section  IV.  But  before  con¬ 
sidering  the  actual  quantitative  results,  it  is  proper  to  examine  the  equations 
for  general  qualitative  characteristics.  This  is  accomplished  in  the  next  sub¬ 
section. 

DISCUSSION  OF  FIELD  EQUATIONS 

For  the  equibeam  distribution,  it  is  clear,  from  Equations  (43)  and  (44). 
that  in  the  far-field  and  within  small  elevation  angles,  the  field  at  the  obser¬ 
vation  point,  P  ,  will  be  circularly  polarized,  not  only  on  boresight  but  for 
any  value  of  0  and  4>  .  Thus,  the  behavior  is  similar  to  that  of  a  flat 
circular  plate  illuminated  by  a  plane  wave,  a  case  that  is  discussed  in  Ap¬ 
pendix  H, 

As  the  range  is  decreased  and/or  as  0  becomes  larger,  there  will  be 
a  certain  amount  of  ellipticity  caused  by  the  z-directed  current  originating 
from  both  orthogonal  polarizations. 

In  contradistinction  with  the  equibeam  case,  the  dipole-like  distributions 
will  produce  different  results  with  respect  to  ellipticity.  Components  X4  and 
Y4  arise  from  cross-polarized  currents  which,  in  turn,  are  caused  by  the 
distortion  of  the  "latitude”  lines  when  projected  onto  the  paraboloid.  Com¬ 
ponents  X2  and  Y2  are  extraneous,  in-phase  components  which  arise  from 
the  same  source.  Those  two  sets  of  components  will  tend  to  become  smaller 
as  the  dish  is  "flatter, "  that  is.  as  the  focal  length  increases.  As  in  the  case 
of  the  equibeam  distribution,  the  X3  and  X5  components  on  the  one  hand, 
the  Y3  and  Y5  on  the  other,  arise  because  of  near-field  and  off-boresight 
effects  combined  with  the  z-directed  current.  Except  in  the  very  near-field 
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or  within  large  elevation  angles,  those  terms  are  small  and  provide  but  small 
corrections  to  the  sidelobe  formation  which  is  determined  primarily  by  the 
exponential  factor. 

For  dipole-like  distributions,  on  boresight  (0  0),  both  the  general 

expressions  and  the  far-field  approximations  for  the  exponential  term  (Equa¬ 
tions  (18)  and  (41))  simplify  considerably  insofar  as  they  become  independent  of 
0  .  It  follows  that 


XI 

Y1 

(61) 

\2 

Y2 

(62) 

\3 

Y3 

(63) 

.\4 

Y4  0 

(64) 

\5 

V  5  0 

(65) 

regardless  of  whether  y  -=  ;  s  1 ,  or  whether  Equation  17  applies  since 


It  follows  that,  on  boresight,  the  field  is  circularly  polarized  if  the  feed  is 
circularly  polarized. 
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Now  consider  the  far-field  equations,  which  are  somewhat  simpler 
when  the  observation  point  is  off-boresight  (0  4  0)  .  The  exponential  term 
(Equation  (41))  is 


exp  !1  _  ej2kf  sin  0  cot  77  2  cos  (0  —  c6)  _  ejacos(0-c5) 
cos  [a  cos  (0  —  <f>)  ]  +  j  sin  la  cos  (0  —  a)  ] 

where  a  stands  for  2kf  sin  0  cot  r/  2  cos  (0  -  0). 

For  0 '  =  0  +  ~ 

exp  1 1  cos  la  cos  (i '/  —  0)1  —  j  sin  I  a  cos  (0—0)1 


(69) 


(70) 


Consequently,  when  the  integrands  are  even  functions  of  0  (as  in  the  case  of 
XI,  Yl,  X2,  Y2),  that  is,  when 

f  (0)  f  (0')  , 

then  the  limits  of  integration  0  (0  to  2n)  can  be  replaced  by  0  to  n  , 
and  exp  !i  by  2  cos  I  a  (<//  -  6)1  .  Note  that  e  ~  1  a  cos  -  tj)  w’ould 

give  the  same  result. 

It  is  clear  that,  for  any  value  of  6  , 

XI  =  Yl  (71) 

provided  that  yc  =  ys  .  But  if  Equation  (17)  applies,  then  Equation  (71) 
does  not  hold,  except  when  0  =  +  it/4  ,  since  cos  0  and  sin  0  now 
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have  the  same  weights.  Furthermore 


\2  -  Y2 

when  6  =  ±  n! 4  and  Equation  (17)  prevails. 

Finally,  even  for  the  pseudo-equibeam  distribution. 


X1'<^>  Y1V*77  2 


\2I  •  ■  X2I  jl  o 

“  O  *  >c6  ±  J7  2 


(72) 


(73) 

(74) 


since 


f(sin20)e-iacos  2)^ 


f  (sin20)  e>  a  sin  ^  di/j  (75) 


Regarding  the  quadrature  terms  X4  and  Y4  ,  it  is  shown,  in  Appendix 
G,  that  they  are  zero  for  6  =  0,  n  2  ,  -  and  3  n/2  . 

It  will  be  shown  that,  when  y  y  1  ,  ellipticity  is  constant  with 
the  azimuthal  angle  6  ,  for  a  given  elevation  0  . 

By  Equations  (215)  and  (216),  of  Appendix  F,  the  circularly  polarized 
waves.  S+  and  S  are,  respectively, 
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j  (X4+  Y4)  ] 


S*  =  y  (X  +  j  Y)  =  —  [(XI  +  X2-  Y1  -  Y2)  + 

=  [  (X2  —  Y2)  +  j  (X4+  Y4)  ] 


S~  «  I  (X  -  j  Y)  =  -i  [(XI  +  X2  +  Y1  +  Y2)  ♦  j  (X4  -  Y4)  I 
-  [2X1  4.  (X2+  Y2)  ] 


S  involves  integration  of 


2  2  ,  , 

(cos  tjj  —  sin^1  t/0  —  j  2  sin  cos  U 


—  —  cos  2  0  —  j  sin  2  i/<  = 


so  that 


oc 


e  j  l  a  cos  (0  —  <^>)  +  2 ifj  ] 


e  j  [a  cos  £  + 


=  —  2 77  e  i2<i&  J2  (a) 


which  is  independent  of  <i  ,  except  for  phase. 


(76) 
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S  involves  integration  of  constants,  so  that 

r  2n 

S-  «  /  ejU  cos  (0  —  ci) ]  d(/,  2 3T  J  (a)  (77) 

*/0 

which  is  independent  of  .  Because  the  phase  of  the  circularly  polarized 
wave  S+  depends  on  2 <5  and  the  phase  of  S  is  independent  of  6  ,  the 
orientation  of  the  ellipse  will  vary  directly  as  ,  although  the  ratio  of 
major  axis  to  minor  axis  is  independent  of  4>  . 

In  general,  ellipticity  can  be  inferred  from  a  linear  polarization  pattern 
(linearly  polarized  feed)  by  superimposing  the  same  pattern  rotated  by  n  2 . 
Thus,  ellipticity  in  the  6  0  and  6  v  2  planes  is  attributable  directly 

to  the  difference  in  "beam width"  in  the  ^  0  and  <i>  n/2  planes  (since 

the  quadrature  terms  are  zero).  Where  such  beamwidths  are  equal,  the 
ellipticity  will  be  1  .  Ellipticity  in  the  <t>  ±  »/4  planes  is  attributable  to 

the  phase-quadrature  components  (since  XI  =  Yl,  X2  =  Y2  at  those 
planes).  It  is  remarkable  that,  for  yc  =  ys  1,  all  those  terms  balance 
one  another,  as  shown  above.  Physically,  this  must  be  so,  since  there  is  no 
"preferred"  plane  regarding  amplitude.  This  property  does  not  hold  for  the 
pseudo-equibeam  distribution  because,  as  discussed  in  Appendix  D,  there  are 
preferred  planes. 

Appendix  H  considers  the  case  of  a  flat  circular  plate,  illuminated  by  a 
plane,  circularly  polarized  wave.  It  is  shown  that  in  that  case,  the  far-field 
is  everywhere  circularly  polarized,  if  the  feed  is  circularly  polarized. 
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DEFOCUSING  IN  FOCAL  PLANE 

Before  concluding  the  analytical  study  of  the  focal  feed  system,  the 
effects  of  defocusing  the  source  (phase  center)  in  the  focal  plane  are  considered 
in  this  subsection.  Defocusing  along  the  axis  of  symmetry  will  be  considered 
in  the  next  subsection. 


Let  the  phase  center  be  located  at  point  O'  ,  in  the  focal  plane,  such 
that  O'  is  determined  by  the  vector  d  in  the  x-direction: 


d 


d 

0 

0 


(78) 


where  d  is  small  compared  to  the  focal  length  f  (Figure  5).  It  will  be 
assumed  that,  because  d  <<  f  ,  the  amplitude  of  II.  is  essentially  the  same 
as  if  the  phase  center  were  located  at  the  origin.  Hence,  H  will  be  given 
by  Equations  (13),  (15),  (18).  However,  because  d  is  not  necessarily  small 
compared  to  a  wavelength,  the  phase  of  H^  at  the  running  point  Q  can  be 
substantially  different. 

For  the  far-field,  the  exponential  term  is  given  by 


e\p. 


j  k  (r  +  =  exp  j  2k  f 


sin  0  cot  t]/2  cos  (0  — 


cos  0  cos  7] 

2  sin ^  77/2 


(79) 


according  to  Equation  (39),  where  the  term  e  has  been  neglected 

(it  is  a  constant),  and  where  $  is  the  distance  from  O'  to  Q. 
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From  Figure  5,  it  is  seen  that 


_f _ 

sin^  i]/2 


sin  77  cos  0  —  —  sin"  77/ 2 
sin  77  sin  0 
cos  77 


(80) 


by  Equations  (8)  and  (78).  The  magnitude  of  JR  is 


sin“  77/ 2 


/  d  v  2 

(sin  77  cos  0  —  y  sin"  77/2  )  +  (sin  77  sin  0) 


2  2 

F  COS  77 


sin"  7/  2 


1  + 


,9  .  4 

d"  sin  *  77/2 


2d  sin"  77/2  sin  77  cos 

f 


1/ 
-1  2 


(81) 


Consequently,  the  exponential  term  becomes 


exp  .  —  j  k  (r  +  )  - 


exp.  j 2 kf  sin  0  cot  rj/2  cos  (0  —  0)  f  - ^ - 

j  2  sin^  77/2 


cos  8  cos  77  —  I  1 


d"  sin  ^  77/2 


d  2 

—  2  y-  sin^  77/2  sin  77  cos  0 


(82) 
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Figure  5.  Defocusing  in  Focol  Plane 


When  d/f  is  very  small  (a  criterion  for  "smallness"  might  he  that 
/  2 

2kf  (d/f)  <  0.  1),  Equation  (81)  can  be  expanded  as  a  power  series  with  the 

second  and  higher  degree  terms  neglected.  Equation  (81)  becomes 

,'i\  - r- -  —  d  sin  r)  cos  0  (83) 

sin“  T]  2 

and  Equation  (82)  becomes 

exp.  —  jk  (r  <-  R)  =  exp.  j  2 k f  £sin  0  cot  77.  2  cos  (0  —  <j>)  4  -^y-  sin  r/  cos  V  J  (84) 
if  d  is  small. 

In  order  to  get  an  estimate  of  ,  the  angle  at  which  the  beam  now 

peaks  up.  consider  that  it  arises  from  the  situation  when  all  the  contributions 
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from  the  differential  areas  dA  are  most  nearly  in  phase.  This  situation  will 
prevail  when 


(85) 


and 


(86) 


where  i)m  is  some  "mean’’  value  of  r)  between  n  and  7/ 


For  r)m  -  180  degrees,  an  estimate  of  0Q  would  be 


"o  °-91 


(f) 


Equation  (86)  states  that  as  the  phase  center  is  shifted  "up,  M  the  beam  shifts 
'  down.  " 

The  sum  of  two  beams  arising  from  a  phase  center  located  at  x  -  d 
and  a  second  phase  center  located  at  x  =  -  d  wrould  have  components  of  the 
form  given  by  Equations  (45)  through  (50),  but  with  the  exponential  term 
replaced  by 


(87) 


Amplitude  and  ellipticity  plots  are  presented  in  Section  IV. 

I)E FOCUSING  ALONG  AXIS 

Let  the  phase  center  be  located  at  point  O'  ,  on  the  axis,  such  that 
O'  is  determined  by  the  vector  d  in  the  z-direction: 


0 


cl 


0 


(88) 


d 
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where  d  is  small  compared  to  the  focal  length  (Figure  6).  It  wall  be  assumed 
that,  because  d  <<  f  ,  the  amplitude  of  IF  is  essentially  the  same  as  if  the 
phase  center  were  located  at  the  origin.  Hence,  H  will  be  given  by 
Equation  (13),  (15)  or  (16).  However,  because  d  is  not  necessarily  small 
compared  to  a  wavelength,  the  phase  of  H.  at  the  running  point  Q  can  be 
substantially  different. 

Let  $  be  the  vector  from  O'  to  the  running  point  Q.  From  Figure  6 
it  is  seen  that 


sm~  7] 


sin  17  cos  O 
sin  7]  sin  0 


cos  77 


sin“  7/ 


2  d 


f 


(89) 


by  Equations  (8)  and  (88).  The  magnitude  of  is 


sin  7/ 


sin“  7/2 


0  .  *2 

(sin  77  cos  u)~  +  (sin  77  sin  1!/)^  + 

(cost,  - 
\ 

sin  ~  7)  2  d  j 

2' 

i  ) 

__ 

/ 

i 

•  1  9  .2  ”i  2 

sin  r/  2  d“  2  cos  77  sin“  7/  2  d 

—7  -  f 


(90) 


The  expression  for  the  near-field  exponential  term  (Equation  (18)) 
becomes 
P-jk  (r+S)  = 


’xp-  -  jk 


)L 


9 

rr  f 


o 


rl 


1  9 

sin  7/  1  sin  7]  1 


[sin  7/  sin  0  cos  (u  —  c6)  *  cos  0  cos  7/  1 


d2  . 


2  f  d 


sin  ^  77  2  sin"-  77/ 2 


cos  77 


'•I 


(91) 
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Figure  6.  Defocusing  Along  Axis 
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Figure  7.  Near-Field  Focusing:  Geometric  Construction 


From  ray  optics  constructions  (Figure  7),  it  is  found  that  the  value  of 
d  that  will  result  in  a  focused  condition  at  r^  is 

f2 

d  =  —  (92) 

rl 

This  equation  is  an  approximation  which  becomes  more  exact  as  the  aperture 
becomes  smaller.  From  a  physical  optics  standpoint,  near-field  focusing  can 
be  construed  as  the  condition  when  the  contributions  from  the  various  differ¬ 
ential  areas,  dA  ,  add  up  almost  in  phase.  Of  course,  this  condition  is  im¬ 
possible  to  achieve,  but  an  attempt  can  be  made  by  imposing  certain  restric¬ 
tions.  Let  A  (r/)  be  the  complete  phase  term  (Equation  (91))  when  0  o 
(focusing  on-axis,  at  range  r^)  : 


A  (t?) 


2  f  r 


rf 


•  1  o 
sin  7/  2 


sin  / / 


r -  cos  rf 

■  o 


d2  4 


4  o  •  2  9 

sin  7}  1  sin  7]  ~ 


COS  7j 


(93) 

When  d  A  d»/lr/(Q)  9  ,  that  is,  when  the  phase  associated  with  a  region 

around  point  Q  does  not  vary,  then  the  contributions  from  the  points  in  that 
region  will  reinforce  at  the  observation  point,  P  . 


Thus, 


dA 

dr, 


1  .  \ 

—  sin  r, 


I L 


sin  X  t1  2 


.  9 

sin  “  t] 


COS  7] 


-f2 


*  r  ]  c  os  7] 

+  - -  f 


sin^  t)/ 2  sin  ^  7\  2 


sin  7 1 


d2  4 


!fd 


sin  X  7]/  2  sin""  77  2 


COS  7] 


1/ 

—  72 


-f^ 


f  d  cos  7/  2  f  d 

+  - : - —  i 


sin^  rj/2  sin  X  7/  2  sin~  77 / 2 


(94) 
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In  order  for  this  equation  to  be  zero, 


jO  f2  2  f  d 

Vi 

(  -f2 

fr,  C°s  r)m 

2fri  \ 

d  +  •  4  /o  .2  -0  C0!5 

sin  r,m/2  sin  rjm  2 

\sintj  2 

■  4  9 

Sln  Vm  2 

sin2^,  2/ 

/  -  f2  fd  cos  2fd  \ 

\sinNm  2  sinlr,m/2  sin2  Vm  *) 


(95) 


where  r?m  is  some  weighted  "mean"  value  of  >7  .  If,  arbitrarily,  r/m  160' 


d  1.065 


rl 


96) 


Amplitude  and  ellipticity  plots  are  presented  in  Section  IV. 
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SECTION  HI 


ANALYSIS  OF  CASSEGRAIN  SYSTEM 


GENERAL  FORMULATION  OF  THE  FIELD  INCIDENT  UPON  THE 
PARABOLOID 


The  general  procedure  that  was  applied  to  obtain  analytic  expressions 
for  the  field  of  a  focal  feed  paraboloidal  antenna  is  now  applied  to  the 
Cassegrain  feed  system.  Actually,  this  procedure  must  be  applied  twice: 
first,  the  field  on  the  paraboloid  is  considered  as  the  field  scattered  by  the 
hyperboloidal  subreflector;  then,  this  first  scattered  field  is  substituted  for 
the  ’'primary  distribution1’  of  the  focal  feed  system  in  order  to  obtain  the 
near-field  and  far-field  of  the  complete  antenna. 


Consider  the  hyperboloidal  subreflector  illuminated  by  a  source  located 
at  one  of  its  foci  (point  0,  Figure  8).  The  second  focus  of  the  hyperboloid 
coincides  with  the  focus  of  the  paraboloidal  dish  at  point  C.  An  arbitrary 
point  on  the  paraboloid  is  still  referred  as  Q  ,  of  coordinates  77  and  0, 
but  the  running  point  on  the  subreflector  has  been  designated  as  P  ,  of 
coordinates  p  ,  0  ,  and  ci>.  Distances  from  the  subreflector  to  the  paraboloid 
are  designated  as  r  and  r^  .  It  is  hoped  that  those  quantities  will  not  be 
mistaken  for  the  coordinates  of  the  entire  antenna  field  observation  point, 
since  the  derivations  of  the  field  on  the  dish  and  of  the  space  field  are  con¬ 
ceptually  separate  problems*  * 


First,  one  wants  to  know  what  is  the  field  H 


S  Q 


at  any  arbitrary  point, 


Q  ,  on  the  paraboloid.  Subsequently,  it  will  become  possible  to  compute  the 
far-field  (and  also  the  near-field)  of  the  paraboloidal  dish. 


* 

The  reason  for  not  introducing  new  symbols  is  the  traditional  association 
of  0  with  elevation,  6  with  azimuth,  r  with  distance,  etc. 
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60  6  ‘02-  VI 


Figure  8.  Cassegrain  Feed  System 
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is* 


An  exact  expression  for 


Q 


where 

A. 

n  = 


H  = 
dA  = 


r 


k 


S 


(97) 


the  outward  normal  unit  vector  at  a  point  Q  on  the  surface 

of  the  hyperboloid, 

the  total  H-field  at  point  P 

the  differential  element  of  area  at  point  P 

the  distance  from  the  observation  point  Q  to  P 

the  propagation  constant,  equal  to  2n  /X 

the  gradient,  taken  at  point  P 

the  total  surface  of  the  hyperboloid 


Equation  (97),  as  it  stands,  is  truly  formidable,  and  some  judicious 
approximations  are  in  order.  The  total  field  H  ,  or  rather  the  expression 

A  — ►  A  — ► 

(n  x  H)  ,  will  be  approximated  by  2(n  x  Hi)  ,  that  is 

(nxH)  -  2  (n  x  flj)  (98) 

where  Hj  is  the  incident  H-field  at  point  P.  Furthermore,  S  is  taken  as 
only  that  portion  of  the  hyperboloid  which  is  directly  illuminated.  Those 
approximations  appear  to  be  valid  if  the  diameter  is  at  least  a  few  wavelengths. 

The  expression  for  Vp  r)  is  exactly 


* 

See  Equation  (1). 
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A 

where  a  ' 
then 


is  the  unit  vector  pointing  from  Q  to  P  . 


\ 


P 


But,  if  r  »  A  , 


(100) 


Again,  this  appears  to  be  a  safe  approximation  in  the  practical  cases  one 
is  liable  to  encounter. 

An  approximation  which  is  not  safe,  and  which  must  be  re-examined 
for  a  specific  geometry,  is  one  consisting  in  letting 

A  P  A 

al  -  ai 

A 

where  a^  is  the  unit  vector  directed  from  Q  to  C  ,  the  focus  of  the  hyper- 

A  A 

boloid  (and  of  the  paraboloid).  The  excursion  angle  from  a^'  to  a^  ,  when 
viewed  from  point  Q  for  any  point  P  on  the  hyperboloid,  may  not  be 
negligible. 

The  distance  from  point  Q  to  point  P  ,  r  ,  is  exactly 

r  =  l^i  +  P  I  =  I  «»i  rj  +  p  |  (101) 


Vector 


is  exactly 


a  1  r  1  +  P 

A  * 

al  rl  *■  P 


(102) 


A 

The  quantities  a^  ,  r  ,  and  p  will  be  obtained  in  the  following 
subsection. 

The  approximation  involving  H.  is  that  it  originates  from  a  point 
source  (phase  center)  and  that  it  is  a  spherical  wave,  with  no  longitudinal 
component.  It  will  be  written 

Hj  =  Hioe_jks  (103) 
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where  the  phase,  ks  ,  is  brought  out  explicitly  (  s  being  the  distance  from 
the  origin,  or  second  focus,  to  the  running  point  P  ).  The  l/s  amplitude 
dependency  will  be  taken  into  account  in  the  overall  amplitude  function  of 
II.  which  is  discussed  in  Appendix  D.  In  practice,  the  only  requirement 
imposed  by  the  approximation  is  that  the  subreflector  be  outside  the  very 
near-zone  of  the  feed. 

By  means  of  Equations  (98),  (100),  (101),  (102),  and  (103),  Equation 
(97)  becomes 

>Vq  ffr a * "-> »  V k <”*r) JA  <m« 


A.  A  — ** 

where  r  and  a  '  are  functions  of  both  P  and  Q  ,  and  n  ,  II.  ,  s  ,  and 
1  10 

dA  are  functions  of  P  only. 


The  vectors  will  be  broken  up  into  Cartesian  components,  and  those 
components  -  in  turn  -  will  be  expressed  in  terms  of  elevation  and  azimuth 
angles,  0  and  0  respectively  for  the  hyperboloid,  ?/  and  0  respectively 
for  the  paraboloid. 


GEOMETRY  OF  THE  PROBLEM 


Point  Q  on  the  paraboloid  (Figure  8)  can  be  described  by  the  elevation 
angles  r?  and  the  azimuth  angle  0  of  the  ray  CQ  ,  where  C  is  at  the  focus 
of  both  paraboloid  and  hyperboloid.  The  quantity  r^  ,  a  function  of  rj  and  of 
f  ,  the  focal  length,  is 


f 

sin  “  (7/ / 2) 


(105) 
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The  unit  vector  is  expressed  in  Cartesian  coordinates  as  * 


sin  T]  cos  ifj 
sin  rj  sin  ip 
cos  7] 


The  polar  equation  of  the  hyperbola  is 


a  (1  -  e^) 

1  —  e  cos  0 


(106) 


(107) 


which  is  derived  in  Appendix  J,  where  the  constants  a  and  e  are  also 
defined. 

The  magnitude  of  vector  p  from  C  to  P  is  derived,  in  Appendix  J, 


as 


p  s  —  2  a 


a  [2  e  cos  0  —  (1  +  e“)  ] 
1  —  e  cos  0 


(108) 


whereas  the  vector  p  itself,  in  Cartesian  coordinates,  is 


1  —  e  cos  d 


(1  —  e“)  sin  0  cos  $ 
2 

(1  —  e)  sin  6  sin  cf> 
[(1  +  e^)  cos  6  —  2  e] 


(109) 


Still  in  Appendix  J,  the  unit  vector  normal  to  the  surface  at  point  Q,  n  , 
is  derived  as 


1 

( 1  4-  e~  —  2  e  cos  0 )  2 


sin  6  cos  $ 
sin  6  sin  $ 
(cos  6  —  e) 


(HO) 


* 

See  Equation  (10). 
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and  the  expression  for  the  differential  area,  dA  ,  is  derived  as 

2  2  2 

d\  ~  — — - -  —•  (1  +  —  2 e  cos  0)  sin  OdO  d6 


(1  —  e  cos  0)^ 


(111) 


It  follows,  by  Equation  (101),  that 


r  magnitude  of 


—  f  a  9 

■7: -  (sin  77  cos  i //)  f  - - —  (1  —  e  )  sin  0  cos  c*> 

•  1  —  e  cos  0 


sin~  7/  2 


-f 

.  9 

sin“  T] 

-f 


3  9 

(sin  7?  sin  i h)  +  - - —  (1  —  e^)  sin  0  sin  <£ 

1  —  e  cos  6 


or 


sin^  r]  2  (1  -  e  cos  0)z  sin"1  t?/2  (1  - e  cos  0) 


-  cos  T) 

sin~  tj  2 


a^[2e  cos  0  —  (1  +■  e“)J 
-  +  - - - 


1  —  e  cos  0 


[0  +  e~)  cos  0  —  2e] 


2  f  a 


and  that 


2_/c 


[  (si  n  r]  sin  6{]  —  e“)  cos  («/>  —  <£) 


+  (1  +  e‘)  cos  6  cos  tj  —  2  e  cos  77)] 


(112) 


-  f 


sin2  77/2 


(sin  77  cos  i/j)  4-  — - 


a  2 

- y  (1  —  e^)  sin  0  cos  <f> 


sin  77  2 
-  f 

sin""  77/2 


f  a  9 

- (sin  77  sin  +  - - - - 77  ( 1  —  e~)  sin  0  sin  <b 


1  —  e  cos  0 


cos  77 


—  e  cos  0 


[(1  +  e2)  cos  0  -  2  c] 


f  a  [2e  cos  6  —  ( 1  +  c  )] 


2  f  £ 


sin  ^77/2 


(1  —  e  cos  0)  sin  77  2(1  —  e  cos  0) 


[sin  77  sin  0  ( 1  —  e  )  cos  (ifj  —  <f>) 

cy 

+  (1  +  e^)  cos  6  cos  77  —  2  e  cos  77  ] 


(113) 
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PRIMARY  PATTERNS 


The  primary  patterns  are  discussed  in  Section  II  and  in  Appendix  D. 
When  r)  (of  Figure  2)  is  replaced  by  0  (of  Figure  8),  the  signs  must  be 
changed.  Another  modification  in  the  analytic  expressions  occur  insofar  as 
the  field  strength  dependence  on  the  range  is  now  (e  cos  0  -  1 )  ,  by 

Equation  (107). 

Thus,  the  equibeam  distribution  is  given  by 


.  2  .9 

Icos^  cos  6  +  sin“  cf>)  cos  cut 

+  (cos  <£  sin  <f>  [cos  0  —  l| )  sin  cot 


^io  ~  ^  ^  (e  cos  ^  1)  (cos  cf>  sin  <f)  [cos  0  —  l] )  cos  cot 

2  9 

+  (sin^  <f)  cos  0  f  cos“  c/>)  sin  cot 


—  sin  6  cos  <£  cos  cot  —  sin  0  sin  if)  sin  cot 


(114) 


whereas  the  dipole-like  and  pseudo-equibeam  distributions  are  given  by 


Uj0  =  h  (tf)  (e  cos  6  -  \) 


yc  cos  0  cos  (Ot 
ys  cos  0  sin  cot 


—  [  (sin  0  cos  if))  yc  cos  cot  +  (sin  0  sin  cf>)  yQ  sin  cot 


(115) 


with 


(116) 
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for  the  dipole-like  distribution,  and 


(117) 


for  the  pseudo-equibeam  distribution. 


In  Appendix  J  ,  it  is  shown  that 


h  (0)  -  cos*”  (20) 


(118) 


fits,  very  closely,  an  experimental  pattern  with  a  12-db  taper  at  the  edge  of 
the  sub-dish  (0Q  =  30  degrees). 

FORMULATION  OF  THE  FIELD  ON  THE  PARABOLOID:  EQUIBEAM 
DISTRIBUTION 

Equations  (107)  through  (114)  or  (115)  can  now  be  substituted  into 
Equation  (104).  The  exponential  term  becomes,  by  Equations  (107)  and  (112) 


exp.  —  jk  (s  +  r)  -  exp#  —  jk 


i  1  — ’  e  cos  0 


j  1  —  e  cos 


- ^ - — - - [  sin  7]  sin  0  ( 1  —  e~)  cos  (i/y  —  <£) 

sin  ~  7]  2(1  —  e  cos  6) 


+•  (1  +  e~)  cos  0  cos  7}  —  2e  cos  7]  ] 


'2I 

\ 


(119) 
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The  cross-product  of  n  and  H.  is  the  vector  whose  components  are 


the  co-factors  of 


sin  d  cos  cf> 

9  9 

(cos'- cos  0  +  sin""  <b)  cos  cot 

+  (cos  <f>  sin  4>  [cos  0  —  1] )  sin  co  t 


sin  0  cos  <£> 

(cos  <f>  sin  cf)  [cos  0  —  1  ]) 

9  9 

(sin^^  cos  0  +  eos“<£  ) 


(cos  0  —  e) 

■  sin  0  (cos  <£  cos  cot 
-♦  sin  ci  sin  cot) 


h  ( 0 )  (e  cos  0—1) 

times  - - 

(1  f  e“  —  2e  cos  0)  2 


This  expression  simplifies  considerably  to  become 


[—  sin  0  cos  <£(1  +  e)  (1  —  cos  0)]  cos  cot  +  [(e  —  cos  0)  —  sin~c£  (1  —  cos  0)  (1  +  e)]  sin  cot 
[sin  cos  0  (1  +  e)  (1  —  cos  0)]  sin  cot  —  [(e  —  cos  0)  —  cos ~0  ( 1  —  cos  0)  (1  t-  e)]  cos  cot 
sin  0  cos  0  sin  cot  —  sin  0  sin  0  cos  cot 


times 


h  (0)  (e  cos  0—1) 
( 1  +  —  2e  cos  0)/2 


(120) 
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or,  symbolically 


where 


and 


(A  4-  B)  sin  cot  4-  C  cos  cot 
(  —  A  4-  D)  cos  cot  —  C  sin  cot 
E  sin  cot  —  F  cos  cot 


A  -  +  h  '  (0)  (e  —  cos  0) 

B  —  h  ^  (0)  ( 1  +  e)  ( 1  —  cos  0)  sin"-  A 

C  =  —  h  '  (0)  (1  +  e)  ( 1  —  cos  0)  sin  A  cos  A 

D  h  '  (0)  ( 1  +  eW  1  —  cos  0)  cos^  A 

E  -  +  h  '  (0)  sin  0  cos  A 
\  -  4  h  '  (0)  sin  0  sin  A 


(121) 


(122) 

(123) 

(124) 

(125) 

(126) 
(127) 


h  '  (0) 


4  h  (0)  (e  cos  0—1) 
(1  4  -  2e  cos  0)  2 


(128) 


Since  further  derivations  of  the  expressions  for  the  field  on  the  para¬ 
boloid  are  common  to  the  cases  of  dipole-like  distributions  (after  the  symbols 
A  through  F  have  been  established),  they  will  be  postponed  to  the  next 
subsection. 

FORMULATION  OF  THE  FIELD  ON  THE  PARABOLOID:  DIPOLE-LIKE 
DISTRIBUTIONS 

The  dipole-like  distributions  Equations  (115)  and  (116)  and  the  pseudo- 
equibeam  distribution  Equations  (115)  and  (117)  will  be  carried  together,  for 
the  sake  of  convenience. 
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The  cross-product  (n  x  H.  )  is  the  vector  whose  components  are  the 


co-factors  of 


h  (0)  (e  cos  0—  1 ) 

( 1  + e"  —  2e  cos  0)  2 


sin  0  cos  sin  0  sin  <,£ 


(cos  0  —  c) 


(yc  cos  0  cos  <ut)  (ys  cos  0  sin  tot)  —  [(yc  sin  0  cos  ci)  cos  tot  *  (ys  sin  0  sin  <i)sin  tot  I 


that  is 


h  (0)  (e  cos  0—1) 

(l  e"  —  2e  cos  0)  2 


—  sin"" 0  sin  c6(yc  cos  c6  cos  tot  +  y^  sin  <3  sin  tot) 
sin  0  cos  (3  (y  cos  ch  cos  tot  4-  ys  sin  <i  sin  tot) 
sin  0  cos  0  (ys  cos  <£  sin  tot  —  yQ 


—  )  s  cos  0  (cos  0  —  c)  sill  cot 
t  y  cos  0  (  cos  0  —  c)  cos  tot 


sin  d>  cos  tot ) 


or,  symbolically 


(129) 


(n  x  Hio) 


(  \  B)  y  sin  tot  +  C  )'c  cos  tot 


{ -  A  D)  y  cos  tot  —  C  y  sin  tot 


E  y  sin  tot  —  Fy  cos  tot 


(130) 


where 


\ 

1-  h '  (0)  cos  0  (e  —  cos  0) 

(131) 

B 

—  h  '  (0)  sin"  0  sin"  0 

(132) 

c;  - 

2 

—  h  *  (0)  siii  0  sin  <£  cos  (/> 

(133) 

D 

+  li  '  (0)  sin"  0  cos"  <f> 

(134) 
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f  h  '  (0)  sin  0  cos  0  cos  c 6 


(135) 


t  li  *  (0)  sin  0  cos  0  sin 


(136) 


and 


-  h  (0)  (e  cos  0—1) 
(1  4  e~  —  2e  cos  0)  2 


as  in  Equation  (128). 

*  — *•  A 

The  cross-product  of  (n  x  Hj0)  and  a  '  (Equation  (113))  is  the 
vector  whose  components  are  the  co-factors  of 

A  A  A 

\  y  / 

I ( A  +  B)  ys  sin  o»t  -t-  C  yc  cos  o»t  ]  [(—A  +  I))  yc  cos  c<jt  —  C  yg  sin  cot]  [Hyg  sin  cot  —  I  yc  coscot] 

alx/  aly  ”  nlz' 


that  is 


(n 


*io>  "  *  1 


a  j  /  '  [(—A  4  D)  yQ  cos  cot  —  C  yg  sin  cot]  —  ai  '  (E  yg  sin  cot  -  I  yc  cos  cot) 

-  a |  ? '  [(A  *  II)  ys  sin  cot  4  Cyc  cos  cot]  4  a]x  (E  Ys  sm  <*>t  “  E  y  cos  cot) 
aly  '  [( A  4  B)  yg  sin  cot  4  C  yQ  cos  cot]  —  a  ^  x  '  [( —  A  4  D)  yc  cos  cot  —  C  yg  sin  cot] 


(137) 


or,  by  rearranging  the  terms 


(n 


fa|z  '  (—A  -*  D)  f  ai  '  F]  yc  cos  cot  -  [ajz '  C  4  a  jv  '  F’]  yg  sin  tot 
[-  ajz'  (A  4  ED  4  ajx'  E]  yg  sin  cot  -  [ajz  '  C  4  ajx'  F]  yQ  cos  cot 
[al\  '  C  -  alx'(-A  4  D)]  yc  cos  cot  4  [a]v'(A  4  B)  4  a]x'C]  ys  sin  cot 
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(138) 


where 


*  1  X 


—  f 


My 


,  1 


9 

sin  ?/  2 


-f 


(sin  7f  cos  (//)  + 


(sin  7/  sin  <//)  ■+ 


1  —  e  cos  0 


sin“  7/ 

,  1  \  ~f 


1  —  e  cos  0 


(1  —  e“)  sin  6  sin  0 


r  )  2 

Jsin  "  ' 


cos  7/ 


'/  - 


1  - - 9 - T  [(1  i  c2)  COS  0  _  2c | 

1  —  e  cos  f/  l 


(139) 


(140) 


(141) 


and  r  is  given  by  Equation  (112). 

The  incident  magnetic  field  vector  at  any  point  Q  on  the  paraboloid 
becomes 


8  Q 


oc 


XI  cos  cjt  +  \2  sin  wt 

\\  sin  col  4  Y2  cos  col 
/]  cos  co t  +  Z2  sin  a>t 


(142) 


where 


\,  -  (4 


6  277  0  r  0Q 

J  yc  exp  I  I  2-  lajx '  ( — A  4-  D)  a^.'F]dA  (143) 

<£  =  0  ^  -  0 


X2  =  (_)  ff  Y*  CXp  1 1  7  la)z  C  '  al/ 


E]  dA 


VI  =  (-t)  ysexptl  ~  t-ajz'(A  +  B)  +  ajx'E 


]  dA 


(144) 


(145) 


112  II  yc  exP  1  I  7  +  alxM  1  dA 


(146) 
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(147) 


exp 


a]x '  (—A  4  D)]dA 


Z2  =  (4-)  J*J*  Ys  exP  ^  “-[aiy'(A+B)  ajx  'C]dA 


(148) 


Instead  of  the  form  given  by  Equation  (104),  it  may  be  advantageous  to 

— j  kr  i 

bring  out  and  suppress  the  phase  factor  e  in  order  to  provide  a  direct 

comparision  with  direct  focal  point  feed.  Equation  (104)  must  be  modified  to 
the  form 


n.Iq -  "io 


v  *  ,  -jk(s  +  r-rj) 
;„)  x  a1  e 


d\ 


(149) 


where  the  phase  of  II  1  is  now  referred  back  to  the  paraboloid  focus. 
Equation  (142)  must  be  modified  according  to 


P 


X  1  /  cos  cot  +  \2'  sin  cot 
Y1 '  sin  cot  +  Y  2'  cos  cot 
Z  1 '  cos  cot  +  Z2' sin  cot 


(150) 


where  the  primes  denote  that  the  exponential  terms  have  been  modified,  in 
Equations  (143)  through  (148),  to  the  form 

—  j  k  (s  +  r  — 1*|) 


It  is  understood  that,  for  the  equibeam  distribution,  the  symbols  A 
through  F  are  given  by  Equations  (122)  through  (127),  and  that 


y  c  =  y*  =  1 
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For  the  dipole-like  distribution,  the  symbols  A  through  F  are  given 
by  Equations  (131)  through  (136),  and 

>'c  =  xs  =  1 

For  the  pseudo-equibeam  distribution,  the  symbols  A  through  F  are 
also  given  by  Equations  (131)  through  (136),  but  yc  and  ys  are  given  by 
Equation  (117). 

Symbols  pertaining  to  the  "modified"  Equations  (143)  through  (148), 
where  the  phase  is  referred  back  to  the  paraboloid  focus,  are  summarized  in 
Table  I. 


PARABOLOID  ILLUMINATION  FROM  SUB-REFLECTOR 

Before  considering  the  space  field  (as  will  be  done  in  the  next  subsection), 
it  is  worthwhile  to  transform  the  field  on  the  paraboloid  in  a  form  that  will 
enable  comparisons  with  measurements. 

Let  the  feed  be  linearly  polarized  by  suppressing  the  sinot  component. 

By  Equation  (150),  the  field  on  the  paraboloid,  in  Cartesian  coordinates,  is 


II 


XI' 

Y2' 


Zl' 


But  Cartesian  coordinates  are  not  appropriate  to  measurements  since  the 

A 

measuring  probe  will  generally  be  aligned  with  the  vector  -  a  (Figure  8). 
Coordinate  transformations  can  be  effected  by  means  of  Equations  (208)  and 
(209). 
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Table  I 


Symbols  for  Modified  Equations  (143)  Through  (148) 


(  f2 


sin  ^  77/2 


a"  [2c  cos  0  —  (1  +  e") 
( 1  —  e  cos  0)2 


2  f  c 


d\ 


sin “77  ‘2(1  —  e  cos  0) 

a2(l-eV 


0  9  I 

[sin  77  sin  0(  1  —  e")  cos  (0  —  0)  ^  ( 1  +  e")  cos  0  cos  17  —  2e  cos  77 1 


(1  -  e 


OV 


(1  +  e“  -  2p  cos  0)  sin  0  d0d0 


e'  p  t  |  -  exp.  -  j  k 


!lx 


llv 


llz 


r 

r 


a  (1  -  r*) 

1  —  e  cos  0 

f 


2  9 

sin  7]  1 


sin"  7]  2 

-  f 


(sin  i)  cos  0)  f  -j 


( 1  —  e")  sin  0  cos  0 


•  9 

_sin"  7]  2 


sin  7] 


e  cos  0 

sin  7]  sin  lb)  +  “j - 77  ( 1  —  e“)  sin  0  sin  0 

1  1  —  e  cos  U 

( 

[  ( 1  -  e")  cos  0  —  2e  ]  > 


COS  7] 


1  —  e  cos  0 


h  '09) 


h  (0) 


(1  +  c"  —  2e  cos 


s  0)  2 

DIPOLE-LIKE 


A 

I] 

c 

I) 

K 

I' 

>'c 

)s 


EQUIBKAM 

li  *  (0)  (e  —  cos  0)  +  h  '  (0)  cos  0  (e  —  cos  0) 

h  '  (0)  (1  f  e)  (1  -  cos  0)  sin^  0  -  h  '  (0)  sin"  0  sin"  0 

h'(0)  (1  +  c)  (1  -  cos  0)  sin  0  cos  0  -  h  '  (0)  sin"  0  sin  0  cos  0 

+  h  /  (0)  sin"  0  cos"  0 

+  h  '  (0)  sin  0  cos  0  cos  0 


PSH  DO- KOI  I  BEAM 


II  '  (0)  (  1  f-  c)  ( 1  —  cos  0)  cos"  0 
h  '  (0)  sin  0  cos  0 


h ' (0)  sin  0  sin  0 
1 
1 


+  h  '  (0)  sin  0  cos  0  sin  0 


9  9  9  _  1 

(cos"0  4  sin"  0  cos4  0)  2 

(cos"0  +  sin"  0  sin"  0 )“  2 
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In  particular,  for  an  H-plane  cut  (<£  =  0  or  n  ),  a  rotation  through 
the  angle  (rj  -  n)  about  the  y-axis  is  in  order.  It  is  described  by  the  rotation 
matrix 


[%-n),y  = 

where  must  be  considered  negative  when  </*  =  n  ,  whereas  for  an  E-plane 
cut  =  ±tt/2)  ,  a  rotation  through  the  angle  (//  -  n)  about  the  x-axis  is 
required.  It  is  described  by  the  rotation  matrix 


COS  T]  0  sill  7/ 

0  1  0 

sin  7 1  0  —  cos  7/ 


(151) 


[  J](r,  -n).x 


\  0  0 
0  —  cos  7]  sin  7/ 

0  —  sin  7/  —  cos  7/ 


Consequently,  the  field  on  the  paraboloid  in  the  H-plane  is 

—  \  1  *  cos  rj  Zl/  sin  7/ 

Y2' 

f  X 1 '  sin  rj  —  /A'  cos  rj 

where  the  lower  signs  apply  when  hj  -  n  , 

In  the  E-plane,  the  field  on  the  paraboloid  becomes 

XI' 

—  Y  2'  cos  rj  +  Z 1  '  sin  q 
—  Y2  '  sin  q  —  Z  1  '  cos  q 


(152) 


(153) 


(154) 
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Other  cuts  can  be  obtained  by  the  appropriate  applications  of  [ T ] 

and  [  T]  . 

7 

Plots  relative  to  H  and  are  given  in  Section  V. 

ANALYTIC  FORMULATION  OF  THE  SPACE  FIELD* 


The  procedure  for  finding  the  space  field  when  the  paraboloid  is  excited 
by  a  Cassegrain  feed  system  is  now  the  same  as  for  the  case  of  a  focal  feed 
system.  But  the  distribution  of  the  incident  field  on  the  paraboloid  is  now 
different,  and  adequate  modifications  must  be  provided.  Essentially,  H{ 


S|Q 

(Equation  142)  is  the  field  incident  on  the  paraboloid  and  now  takes  the  place 
of  H  in  Equation  (5). 


Thus,  the  starting  point  is  Equation  (6)  modified  with  regard  to  both 

modulus  and  argument  of  H  .  If  one  uses  the  "modified"  form  of  H. 

i  1 

(Equation  149)  where  the  phase  of  H  is  referred  back  to  the  focus  of  the 
paraboloid,  then  Equation  (18)  is  the  right  expression  for  the  exponential 
factor 


-XP-  -  .)■< 


c  l 

II 


'i2* 


'i 


.  4 

sin  7)/ 2 


•  2  0 
sin  j]  A 


[sin  rj  sin  6  cos  (0  —  (/>)  +  cos  0  cos  77] 


The  field  on  the  paraboloid  is  given  by  Equation  (150).  The  expression 

A  A  - 

for  n  is  given  in  Equation  (7).  Consequently,  (n  x  H  )  is  the  vector 
formed  by  the  cofactors  of 


* 

In  this  and  the  following  sections,  the  symbols  refer  back  to  Figures  2  and  3 
since  point  P  is  no  longer  on  the  hyperboloid,  but  is  now  an  arbitrary 
observation  point  in  space. 
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—  (cos  0  cos  77  2)  —  (sin  0  cos  77  2)  sin  7 )  2 

I  \  1 '  cos  cot  -4  \2'sin  cot]  [A  1 '  sin  cot  4  ^  2' cos  cot]  [  Z 1  '  cos  tot  +  /2 '  sin  <ot) 


or 


(n 


(—  sin  0  cos  77  2)  (71  'cos  cot  +  Z2  'sin  cot) 
(cos  0  cos  77  2)  (Zl  'cos  cot  4  Z2  'sin  cot) 

(  —  cos  0  cos  77  2)  (\  1  'sin  cot  4  Y2'cos  cot) 


or,  in  symbolic  form 


—  sin  7/  2  (\  1  'sin  cot  4  ^1  2 'cos  cot) 

-t-  sin  77  2  (Xl  coscot  +  \2 'sin  cot) 

+  (sin  0  cos  7/  2)  ( \  1  'cos  cot  +  \2'sin  cot) 

(155) 


(n 


\j  sin  cot  4  A 9  cos  cot 
Bj  cos  cot  +  Br>  sin  cot 
Cl -j  cos  cot  ■+■  C9  sin  cot 


(156) 


where 

A  |  ---  —  (\  j '  sin  rj  2  Z2  sin  0  cos  7/  2) 

A  9  -  —  ( Y9 '  si  n  77  2  4  Z  ^ '  sin  0  cos  77  2) 

B|  -  4  (Xj'sin  77/2  4*  Zj'cos  0  cos  77  2) 

B  y  -  +  (\  ■>'  sin  77  2  4  Z0'  cos  0  cos  77  2) 

Cj  -  4  (\  |  '  sin  0  cos  77/2  —  'i  9 '  cos  0  cos  77  2) 

C  9  -  4  (\  9 'sin  0  cos  77  2  —  Vj  '  cos  u  cos  77/2) 

Ak 

An  expression  for  a  •  is  provided  by  Equation  (12).  Consequently,  the 

A  - ►  A 

cross-product  (n  x  H^0)  x  a^1  is  the  vector  whose  components  are  the 
co-factors  of 
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A 

V 


A 

/ 


sin  cot  +  A o  cos  cot) 


(B|  cos  cot  *■  Bg  sin  c/it)  ((  |  coswt  *  C.osiiwoO 


cot  T)  2  cos  k!j  —  sin  6  cos  0^  co^r/^s’n^  “  sinf/sin^ 


that  is,  (n  x  !lio)  -  a |  ' 

IB]  cos  to  t  4.  Bo  sin  C/it)  ■+■  (C  j  cos  cot  +  C o  sin  (at)  ^ —  cot  ?/  2  sin  0  —  sin  $  sin  0^ 

t  \  I  sin  cot  t-  Ao  cos  cot)  ■+-  (Cj  coscot  *■  Co  sin  cot)  cot  rj  2  cos  0  —  sin  0  cos 

(Aj  sin  cot  t  \2c0sc0t)  cot?/  2sin  0  —  sin  #sin  0^ 

-  (Bi  coscot  -f  Bo  sin  cot)  (  —  cot?/  2  cos  0  —  sin  0  cos  0^ 

1  “  \ri  / 


(157) 


Finally,  after  substitution  of  the  expression  for  dA  as  given  by 
Equation  (9),  the  components  of  the  field  at  point  P  assume  the  following 
forms: 


\P  -  coscot 


J  J  (-)  13  |  •  Cj  COt  /,  2  si 


0  % 


sin  0  —  sin  0  sin  0 


expH-^_dr,M 
sin  ?/  2 

(158) 


\Q  sin  cot  J'J' t— )  j^Bo  +  Co  cot?/  2  sin  0  —  sin  0  sin  6^ 

JfbM',  cot/,  2co8*-sin0cos^ 


C\p  II 


sin  ?/ 
sin3?/  2 


dr)  (10 


^  P  sin  cot 


^  0  coscot 


cot  ?/  2  cos  0  —  sin  0cos  0 


) 


c\p  11  — —  dr,  dt'/ 
sin  '/,  2 


exp  ||  -  S‘.n  ’l  dr,  <k 
sin'*  r,  2 


(159) 


(ir,0) 


(161) 


and  z-components  which  are  neglected  for  0  *  5  degrees. 
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Table  II  summarizes  the  symbols  used  in  Equations  (158)  through  (161). 


Table  II 


Polarization  effects  at  point  P  can  be  computed  by  means  of  Equations 


(53)  through  (60),  where  now  the  total  x-component  of  Hc 
polarized  feed  system)  is 


(circularly 


\  (\P  jXQ) 


(162) 


and  the  total  y-component  is 


)  =  (j  YP  +  YQ)  (163) 

whereas,  if  the  primary  pattern  is  linearly  polarized  (for  instance,  only 
x-eomponent  applied), 
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\  \P 


(164) 


V'  =  ^  Q  (165) 

Amplitudes  and  ellipticities  for  the  Cassegrain  sj  stem  have  been 
computed  under  various  conditions  and  plotted.  The  results  arc  found  in 
Section  V. 
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SECTION  IV 


COMPUTED  RESULTS  -  FOCAL  FEED  SYSTEM 


AMPLITUDES 

Figures  9  through  25  are  plots  of  the  overall  field  amplitude  at  the 
observation  point  P  as  a  function  of  elevation  angle,  6  ,  when  the  feed 
system  is  linearly  polarized.  Those  plots  are  computed  values  of  Equation 
(59).  Those  are  plots  of  directivity  rather  than  of  gain,  since  no  attempt 
was  made  to  integrate  and  normalize  the  primary  illumination.  The  peaks 
of  the  beams  are  always  taken  as  the  0  db  level,  although  —  in  certain 
eases  -  gain  differences  are  indicated  (Figure  16.  18  and  20  through  25). 

The  actual  orientation  of  the  total  fit  Id  at  point  P  is  considered  in  detail 
in  Figures  42  through  49.  Quantities  constant  from  graph  to  graph  are  the 
frequency  (co  =  1288  Me),  and  the  paraboloid  aperture  (diameter  =  30  ft). 

Figures  9  through  14  show  amplitude  variations  in  the  far-field,  from 
graph  to  graph,  as  the  primary  illumination  is  varied  in  its  basic  distribution 
(equibeam,  dipole-like,  pseudo-equibeam)  and  the  tapering  function  h (77 )  . 
Figure  15  pertains  to  changes  in  focal  length.  Figures  16  through  19  reflect 
amplitude  variations,  from  graph  to  graph,  caused  by  changes  in  the  range 
of  the  observation  point  P  .  Finally,  Figures  20  through  25  reflect  amplitude 
variation  ehanges  as  the  surface  of  the  paraboloid  is  perturbed. 

2  , 

Figure  9  whose  tapering  function  is  given  by  sin  r\/ 2  h(r])  = 

2  1.5 

sin~7]/2  cos  *  (7r  -  77)  (10  db  taper,  exclusive  of  the  paraboloid  range  factor 
2 

sin~7]/2)  ,  can  be  considered  as  the  ease  resulting  from  some  MstandardM 
primary  illumination:  equibeam,  with  taper  at  the  edge  of  the  dish  optimum 
for  gain,  sidelobe  level,  and  directivity  considerations.  The  far- field  beam 
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ELEVATION,  0  -  DEGREES 


Figure  9. 


Focal  Feed  System:  Amplitude  Related  to  Primary  Distribution 
(Far-Field) 
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is  perfectly  symmetrical  (no  azimuth  variations)  over  the  range  of  elevations 
considered.  The  3-db  half-beamwidth  is  0.  86  degrees.  This  graph  ean  be 
referred  to  Figures  13  and  14  which  pertain  to  the  same  tapering  function, 
but  to  dipole-like  and  pseudo-equibeam  basie  distributions,  respectively. 

In  Figure  10,  the  basic  distribution  is  dipole- like  (H-field  direction 

along  parallel  "latitude1'  lines:  H-field  amplitude  varying  with  azimuth). 

There  is  no  taper,  and  the  factor  reflecting  range  variations  from  the  phase 

2  , 

center  to  the  various  points  on  the  paraboloid  (factor  sin  rj/2)  is  not  taken 
in  consideration.  The  plots  of  this  graph  show  that,  in  the  far-field,  the 
beam  is  not  cylindrical.  Since  the  paraboloid  is  illuminated  more  effec¬ 
tively  in  the  H-plane  than  in  the  E-plane,  the  beam  is  narrower  in  the  H- 
plane  than  in  the  E-plane,  The  H-plane  3-db  half-beamw  idth  is  0.74  de¬ 
grees,  and  the  E-plane  3-db  half-beamwidth  is  0.825  degrees.  Both 
values  are  smaller  than  the  value  of  Figure  9,  because  of  the  absence  of 
tapering. 

In  Figure  11,  the  distribution  is  dipole-like  and  the  tapering  function  is 
2 

implied  in  sin  77/ 2  h(??)  =  cos  r/1 

whore  =  2 0;  '  „  )  <”  +  *  -  2V 

o 

The  illumination  function  is  essentially  the  same  as  in  Figure  10,  except 
that  it  is  multiplied  by  an  overall  eosine  function  which  tapers  the  amplitude 
to  zero  at  the  edge  of  the  dish.  This  tapering  results  in  an  increase  in 
beamwidth.  The  Il-plane  3-db  half-beamwidth  is  now  0.  96  degrees,  whereas 
tlic>  E-plane  3-db  half-beamwidth  is  1.  05  degrees.  Not  shown  in  the  grapli  is 
the  fact  that  the  sidelobe  level  is  down,  as  compared  to  Figure  10  (about 
10  db  down,  for  0  =  4.5  degrees). 
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AMPLITUDE -DB 


Figure  10.  Focal  Feed  System:  Amplitude  Related  to  Primary  Illumination 
(Far-Field) 
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Figure  11.  Focal  Feed  System:  Amplitude  Related  to  Primary  Illumination 
(Far-Field) 
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2  2 

In  Figure  12,  the  distribution  is  dipole-like  and  sin  Tj/2  h (77)  =  cos  r\'  . 
The  primary  illumination  has  the  same  character  as  in  the  preceding  case, 
except  that  it  is  even  more  concentrated  around  the  center  of  the  paraboloid 
at  the  expense  of  the  periphery  region.  The  half-beamwidths  are  even  wider: 
for  the  H-plane,  it  is  1. 15  degrees,  and  for  the  E-plane,  1.24  degrees.  Not 
shown  in  the  graph  is  the  fact  that  the  sidelobe  level  becomes  very  low  (about 
GO  db  down,  for  0  =  5  degrees). 

In  Figure  13,  the  illumination  is  dipole-like  and 

sin"*  r//2  h  (77)  ~  sin^/2  cos^  (n  —  77) 

This  illumination  corresponds  to  a  10  db  taper  (exclusive  of  range  effect)  as 
in  Figure  9,  but,  because  of  the  cos  77  variation  in  the  E-plane,  the  actual 
taper  in  the  E-plane  is  16.  5  db.  The  3-db  half-beamwidth  for  the  H-plane, 
which  is  0.  84  degrees,  is  comparable  to  that  of  Figure  9.  hi  the  E-plane, 
the  3-db  half-beamwidth  is  0.95  degrees. 

In  the  preceding  graphs  dealing  with  the  dipole-like  distribution  (Figures 
10  through  13),  the  beamwidth  increases  continuously  from  the  H-plane  to 
the  E-plane  where  it  becomes  maximum.  The  beamwidth  at  the  45  degree 
plane  is  about  half-way  between  the  two  extremes.  In  the  case  that  follows, 
the  situation  is  different  and  somewhat  like  that  of  Figure  9. 

hi  Figure  14,  the  dipole-like  distribution  has  been  made  "pseudo- 

equibeam"  by  making  the  amplitude  of  H  independent  of  azimuth.  However, 

the  direction  of  the  H-field  is  still  held  along  parallel  "latitude"  lines,  just 

as  in  the  dipole-likc  distribution.  As  explained  in  Appendix  D,  the  45-degree 

planes  on  the  paraboloid  still  contain  "cross-polarized"  current  components. 

2  ,  2  ,  1.5 

The  tapering  function  is  again  sin  r\/2  h(r;)  =  sin  77/2  cos  (77  -  77)  , 

and  now  corresponds  to  a  10  db  taper  (exclusive  of  range  effect)  at  the  edge 
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ELEVATION,  6  -  DEGREES 


Figure  12.  Focal  Feed  System:  Amplitude  Related  to  Primary  Illumination 
(Far-Field) 
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Figure  13.  Focal  Feed  System:  Amplitude  Related  to  Primary  Illumination 
(Far-Field) 
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Figure  14.  Focal  Feed  System:  Amplitude  Related  to  Primary  Illumination 
(Far-Field) 
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of  the  dish,  independent  of  azimuth.  It  is  noted  that  the  H-  and  E-plane 
3-db  half-beamwidths  are  equal  to  each  other  and  to  the  value  obtained  in 
Figure  9.  However,  about  12  db  below  the  peak  of  the  beam  on  boresight, 
the  beam  no  longer  exhibits  cylindrical  symmetry.  It  broadens  more  in 
the  45-degree  planes  than  in  either  the  H-  or  E-planes. 

In  Figure  15,  interest  shifts  to  the  effects  of  focal  length  on  amplitude 
for  the  dipole-like  distribution.  The  plots  pertain  to  paraboloids  having  the 
same  apertures  as  in  the  preceding  cases  (30  ft),  but  with  focal  lengths 
increased  from  12.  5  ft  to  29.  7  ft  and  to  68  ft.  By  Equation  (176),  it  is  found 
that  the  aperture  angles,  rf  ,  are  respectively  151.6  and  167.4  degrees. 

The  tapering  used  in  the  plots  is  of  the  form 

O  O  1  r 

sin“  r]  2  h  (if)  sin“  77/2  cos  (77  —  7/) 

to  correspond  to  Figure  13.  However,  the  77-variation  is  so  small  (from 
151.  G  degrees  and  167.4  degrees  to  180  degrees)  that  Figure  15  is  better 
referred  to  Figure  10.  Because  the  illumination  is  now  almost  constant  and 
equal  in  both  H-  and  E-planes  and  thus  approaches  the  equibeam  condition, 
the  far-field  H-  and  E-planc  patterns  are  essentially  equal  and,  for  the 
68— ft  focal  length  case,  the  beam  exhibits  cylindrical  symmetry  down  to 
30  db.  The  3-db  half-beamwidths,  which  are  0.80  degrees,  fall  between 
the  H-plane  and  E-plane  half-beamwidths  of  Figure  10. 

Figures  16  through  19  exhibit  the  effect  of  range  at  the  observation  point, 
P,  on  the  amplitude  pattern. 

Figure  16  (equibeam,  10  db  taper;  see  Figure  9)  compares  the  ampli¬ 
tude  in  the  far-field  to  the  amplitude  at  a  400-ft  range.  At  400  ft,  the  3-db 
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ELEVATION,  9  -  DEGREES 


Figure  15.  Focal  Feed  System:  Amplitude  Related  to  Focal  Length  (Far-Field) 
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ELEVATION,  6  -  DEGREES 


Figure  16.  Focal  Feed  System:  Amplitude  Related  to  Range 
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half-beamwidth  is  increased  from  0.  86  degrees  to  9.95  degrees.  The  null 
is  filled  in  and  shifted  toward  a  larger  elevation  angle.  Gain  on  boresight 
is  down  1.  65  db. 

Figure  17  corresponds  to  the  same  illumination  as  Figure  10,  but  the 
range  of  point  P  is  taken  at  400  ft,  300  ft,  and  200  ft.  Near  boresight,  and 
roughly  within  the  3-db  points,  the  character  of  the  curves  is  about  the  same 
for  ranges  of  400  ft  and  300  ft  as  for  the  far-field.  However,  the  far-ficld 
nnullsn  fill  in  rapidly  as  the  range  decreases.  For  the  200-ft  range,  the 
character  of  the  curve  is  lost,  and  the  former  ’’nulls”  are  almost  as  high 
as  the  beam  on  boresight,  in  anticipation  of  ’’beam  splitting”  at  a  probably 
slightly  closer  range. 

Figure  18  compares  the  amplitude  pattern  at  a  range  of  400  ft  to  the 
amplitude  pattern  in  the  far-ficld,  for  the  dipole-like,  10-db  taper  illumina¬ 
tion.  The  near-field  pattern  is  characterized  by  a  broadening  of  the  3-db 
beamwidth  and  a  ”filling-in”  of  the  null  region. 

Figure  19  compares  the  amplitude  pattern  at  a  range  of  400  ft  to  the 
amplitude  pattern  in  the  far-field  for  the  pseudo-equibeam,  10  db  taper 
illumination.  Qualitatively,  the  same  characteristics  of  Figures  16  and  18 
apply  to  this  graph:  slight  broadening  of  the  beamwidth,  ”filling-in”  of  the 
null  region. 

Figures  20  through  25  reflect  amplitude  variations  (or  lack  thereof) 
when  the  paraboloidal  surface  is  perturbed. 

In  Figures  20,  21,  and  22,  randomly  distributed  imperfections  on  the 
paraboloidal  surface  are  simulated  by  letting  the  focal  length,  f,  become 
a  function  of  rj  and  ip  .  Thus,  in  the  integration  process,  to  each  point 
(77.  ,  if.)  corresponds 
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ELEVATION,  6  -  DEGREES 


Figure  17.  Focal  Feed  System:  Amplitude  Related  to  Range 
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Figure  18.  Focal  Feed  System:  Amplitude  Related  to  Range 
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ELEVATION,  6  -  DEGREES 


Figure  19.  Focal  Feed  System:  Amplitude  Related  to  Range 
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Figure  20.  Focal  Feed  System:  Amplitude  Related  to  Paraboloid  Imperfections 
(Far-Field) 
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Figure  21.  Focal  Feed  System:  Amplitude  Related  to  Paraboloid  Imperfections 
(Far-Field) 
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ELEVATION,  8  -  DEGREES 


Figure  22.  Focal  Feed  System:  Amplitude  Related  to  Paraboloid  Imperfections 
(Far-Field) 
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The  function 


'i,j  ><„.*,) 


.  (1  .  <|j)  12.5  (t 

has  been  generated  and  stored,  and  it  is  the  same  one  that  is  used  in 
Figures  20  through  22.  The  error  (or  imperfection)  function  <jj  is  randomly 
distributed,  but  limited  such  that 

l«jjl  <  0.005 

Consequently,  the  maximum  deviation  on  f  is 
0.005  x  12.5  ft  *  ±  0.75" 
or  ±  0. 082  A 

or  ±  30  degrees 

The  far-field  approximation  equations  must  be  used  cautiously.  Thus, 

*”  i 

the  term  e  (neglected  from  Equation  (41))  must  be  brought  back  under 

the  integral  sign,  since  it  is  no  longer  a  constant. 

The  illumination  functions  used  in  Figures  20,  21  and  22  are,  respec¬ 
tively,  those  of  Figures  10,  11,  and  12.  The  amplitude  patterns  are 
strikingly  similar  to  those  of  the  unperturbed  surfaces  (same  beamwidth, 
same  null  regions).  However,  the  level  of  the  curves  (the  gain)  is  down 
1. 8  db  from  the  levels  of  the  unperturbed  cases.  * 


*SeeDr.  Ruze's  findings,  NEREM  Record,  1964,  pp.  166-167. 


90 


For  Figures  23,  24,  and  25,  the  focal  length  was  made  to  vary  discon- 
tinuously  as  a  function  of  ip  only,  in  a  random  manner.  Thus,  each  equal 
pie-shaped  slice,  covering  15  degrees,  has  a  random  (but  bounded)  focal 
length,  given  by 

fi  =  f(lpi)  ,i  =  1,2 . 24 

=  [  1  +  €(ipi)]  12.5  ft 

Since  there  are  only  24  values  of  [  1  +  e(i/d)]  ,  they  are  listed  in  the 
following  table: 


Table  III 


Values  of  [  1  +  e(#)]  * 


0.998878 

0.995612 

1.003482 

0. 995667 

0.995607 

0.995452 

0.999624 

0.999905 

0.  995332 

0. 999547 

1. 000148 

0.998926 

0.  998794 

1. 002547 

1. 000484 

0. 998315 

1. 003611 

0.998473 

0.  997937 

0. 998034 

1. 004276 

1.000725 

1. 001686 

1. 002419 

The  illumination  functions  used  in  Figures  23,  24  and  25  are,  respec¬ 
tively  ,  those  used  in  Figures  10,  11,  and  12.  The  amplitude  patterns  show 
gross  distortions  as  the  amplitude  varies  considerably  as  a  function  of  <p , 
although  it  is  noteworthy  that  the  "average"  3-db  half-beamwidths  remain 
about  the  same  as  in  the  corresponding  cases  of  the  perfect  paraboloid.  On 
boresight,  the  amplitudes  are  down  1.  35  db  from  those  corresponding  to  the 
perfect  paraboloid.  Another  salient  feature  of  all  three  plots  is  that  the 
maximum  amplitude  is  reached  for  9  =  0.  25  degrees  rather  than  on  bore- 
sight.  Hence,  a  "pointing  error"  is  introduced  by  the  pie-shaped  segments 
of  unequal  focal  length. *  ** 

*  Read  across. 

**The  reason  is  that  the  24  values  of  Table  I  are  "randomly  "but  not  "evenly" 
distributed. 
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Figure  23.  Focal  Feed  System:  Amplitude  Related  to  Paraboloid  Imperfections 
(Far-Field) 
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Figure  24.  Focal  Feed  System:  Amplitude  Related  to  Paraboloid  Imperfections 
(Far-Field) 
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Figure  25.  Focal  Feed  System:  Amplitude  Related  to  Paraboloid  Imperfections 
(Far-Field) 
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ELLIPTICITIES 


Figures  26  through  42  are  plots  of  ellipticity  at  the  observation  point, 

P  ,  as  a  function  of  the  elevation  angle,  9  ,  when  the  feed  system  is  circu¬ 
larly  polarized.  These  plots  are  computed  values  of  Equation  (56).  Quanti¬ 
ties  constant  from  graph  to  graph  are  the  frequency  (a  =  1288  Me),  and  the 
paraboloid  aperture  (diameter  =  30  ft). 

Figures  26  through  31  show  ellipticity  variations,  in  the  far-field,  from 
graph  to  graph,  as  the  primary  illumination  is  varied.  The  various  illumina¬ 
tion  conditions  correspond,  respectively,  to  Figures  9  through  14.  Figure  32 
pertains  to  a  change  in  focal  length  and  is  the  companion  to  Figure  15,  which 
depicted  amplitude  changes.  Figures  33  through  36  reflect  ellipticity  varia¬ 
tions,  from  graph  to  graph,  caused  by  changes  in  the  range  of  the  observation 
point,  P.  Finally,  Figures  37  through  42  reflect  changes  in  ellipticity  as  the 
surface  of  the  paraboloid  is  perturbed. 

1.  5 

Figure  26,  where  the  taper  function  is  h (77)  =  cos  (tt  -  77)  ,  can 
be  considered  as  the  case  resulting  from  the  "ideal*'  primary  illumination* 
equibcam  distribution,  with  a  10  db  taper  at  the  edge  of  the  dish  (exclusive 
of  the  effect  of  range  variations  from  the  phase  center  to  the  various  points 
on  the  paraboloid).  Ellipticity  in  the  far-field  is  0  db  on  boresight  and 
remains  less  than  0.  03  db  at  the  3-db  half-beamwidth.  Although  there  are 
no  cross-polarized  current  components  on  the  paraboloid,  there  are 
z-directed  currents  whose  origins  partake  of  both  quadrature  components. 

As  a  result,  ellipticityr  increases  as  the  elevation  angle  6  opens  up. 

Although  the  extraordinar}'  circular  component  is  small  and  increases  slowl} 
with  0  ,  ellipticity  is  enhanced  near  the  null  where  the  ordinary  circular 
component  dips  to  low  values. 
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Figure  26.  Focal  Feed  System:  Ellipticity  Related  to  Primary  Distribution 
(Far-Field) 

9 

In  Figure  27.  the  basic  distribution  is  dipole-like,  and  sin“?7/2  h(r?) 

=  1,  as  in  Figure  10.  As  shown  in  Section  II,  the  ellipticity  is  0  db  on  bore- 
sight,  and,  for  a  given  elevation  angle,  is  independent  of  the  azimuth  angle 
op.  The  ellipticity  increases  first  slowly,  and  then  more  rapidly,  away  from 
boresight.  It  is  about  0.  8  db  for  the  average  3-db  half-beamwidth.  (Note 
the  change  of  scale  from  Figure  26. )  Ellipticity  becomes  very  large  and 
tends  toward  infinity  (linear  polarization)  at  the  elevation  region  corres¬ 
ponding  to  the  first  null.  Then,  ellipticity  decreases  and  tends  again  toward 
0  db,  as  the  observation  point  moves  toward  the  peak  of  the  first  sidelobe. 

In  Figure  28,  the  distribution  is  dipole-like,  and  the  taper  function  is 
2 

implied  in  sin~?7/2  h{rj)  =  cos  7 f  ,  as  in  Figure  11.  Again,  ellipticity  is 
0  db  on  boresight,  and  again  it  is  independent  of  <p  for  any  given  value  of 
0 .  The  tapering  of  the  illumination  function  proves  beneficial  to  the 
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Figure  28.  Focal  Feed  System:  Ellipticity  Related  to  Primary  Illumination 
(Far-Field) 
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ellipticity  since  the  region  near  the  periphery  of  the  dish  which  contributes 
most  to  the  ellipticity  (Equations  (47)  and  (50))  is  less  fully  utilized.  Ellip¬ 
ticity  corresponding  to  the  average  3-db  half-beamwidth  is  now  0.  55  db. 

In  Figure  29,  the  distribution  is  dipole-like,  and  the  taper  function  is 
2  2 

sin  77/ 2h(r7)  =  cos  77'  ,  as  in  Figure  12.  Ellipticity  is  0  db  on  bo  resight  and 
is  independent  of  cp  for  any  value  of  6  .  The  tapering  of  the  illumination 
function  being  even  more  pronounced  than  in  the  preceding  graph,  the 
ellipticity  corresponding  to  the  3-db  half-beamwidth  is  only  0.  4  db. 

In  Figure  30,  the  distribution  is  dipole-like  and  the  taper  function  is 
1 . 5 

h(7])  =  cos  (7r  -  77)  as  in  Figure  13.  Ellipticity  is  0  db  on  boresight  and 

is  independent  of  (p  for  any  value  of  6.  The  primary  illumination  being 
better  utilized  than  in  Figures  28  and  29,  but  not  quite  as  well  as  in  Figure 
27,  the  ellipticity  corresponding  to  the  3-db  half-beamwidth  is  0.  7  db.  It 
tends  to  increase  near  the  null  and  decrease  toward  0  db  for  an  elevation 
angle  corresponding  to  the  peak  of  the  first  sidelobe. 

In  Figure  31,  the  distribution  is  now  pseudo-equibeam,  and  the  taper 
1.  5 

function  is  h(77)  =  cos  (n  -  77)  as  in  Figure  14.  The  character  of  the 
ellipticity  versus  elevation  angle  is  changed  drastically.  It  is  now  a  function 
of  cp  and  is  essentially  0  db  in  the  H-  and  E-planes  for  any  value  of  6  . 
However,  it  reaches  a  maximum  value  in  the  45-degree  planes:  there  are 
now  preferred  planes  corresponding  to  the  planes  on  the  spherical  wavefront 
where  the  field  is  circularly  polarized  (quadrature  components  equal  and 
orthogonal)  and  other  planes  where  the  ellipticity  is  maximum  (quadrature 
components  equal  but  of  maximum  space  angle).  It  is  noteworthy  that  the 
maximum  ellipticity  is  essentially  that  of  Figure  30  and  that,  in  particular, 
the  ellipticity  corresponding  to  the  3-db  half-beamwidth  is  also  0.  7  db. 
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Figure  29.  Focal  Feed  System:  Ellipticity  Related  to  Primary  Illumination 
(Far-Field) 


Fiqure  30.  Focal  Feed  System:  Ellipticity  Related  to  Primary  Illumination 
(Far-Field) 
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Fiqure  31.  Focal  Feed  System:  Ellipticlty  Related  to  Primary  Illumination 
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Figure  32  relates  the  ellipticities  of  three  paraboloidal  antennas  of  the 
same  aperture,  but  of  different  focal  lengths  (12.5  ft,  29.  7  ft,  and  68  ft, 
respectively).  The  illumination  is  that  of  Figure  15.  As  the  aperture  angle 
becomes  smaller,  there  is  less  difference  between  the  dipole-likc  and  the 
cquibeam  distributions,  as  noted  in  Appendix  D.  Corresponding  to  the  3-db 
half-beamwidth,  ellipticity  is  only  0.  03  db  for  the  68—  ft  focal  length  system, 
and  0. 15  db  for  the  29. 7— ft  focal  length  system,  as  compared  to  0.  7  db  for 
the  12.  5— ft  focal  length  system. 

Figures  33  through  36  reflect  the  effect  of  the  range  of  the  observation 
point  P  on  the  ellipticity  pattern. 

Figures  32  through  36  reflect  the  effect  of  the  range  of  the  observation 
10  db  taper  and  compares  ellipticities  at  a  range  of  400  ft  to  ellipticities 
in  the  far-field.  Remarks  applying  to  Figure  27  apply  in  the  near-field, 
except  that  since  the  ordinary  circular  component  of  the  field  does  not  dip 
as  much  as  in  the  far-field,  the  ellipticity  peaks  up  to  a  small  value  near 
the  null. 

Figure  34  corresponds  to  the  same  illumination  as  Figure  27,  but  the 
range  of  point  P  is  taken  at  400  ft,  300  ft,  and  200  ft.  On  boresight,  the 
ellipticity  is  0  db  in  all  cases.  Near  boresight  and  within  the  3-db  half- 
beamwidth,  the  ellipticities  are  about  the  same  as  in  the  far-field  case. 

Note,  however,  that  the  ellipticity  for  the  200-ft  range  is  somewhat  higher, 
near  boresight,  in  anticipation  of  the  near-by  "dip"  in  amplitude  (Figure  17). 
In  the  region  of  elevation  angles  corresponding  to  the  far-field  null,  the  near- 
ficld  ellipticities  are  not  as  high  as  in  the  far-field  case.  Indeed,  the  closer 
the  range  and  the  shallower  the  amplitude  "dip,  "  the  lower  the  maximum 
ellipticity. 
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Figure  32.  Focal  Feed  System:  El lipticity  Related  to  Focal  Length 
(Far-Field) 
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Figure  33.  Focal  Feed  System:  El  I i pti city  Related  to  Range 


Figure  35  corresponds  to  the  illumination  of  Figure  30,  but  the  range 
of  point  P  is  taken  at  400  ft.  Ellipticity  is  essentially  the  same  as  in  the 

far-field  within  the  3-db  half-beamwidth,  but,  for  elevation  angles  corres¬ 

ponding  to  the  null  of  the  far-field,  the  ellipticity  at  400  ft  is  far  smaller, 

in  accordance  with  the  shallow  amplitude  dip  of  Figure  18. 

Figure  36  corresponds  to  the  illumination  of  Figure  31  (pseudo- 

equibeam  with  10-db  taper  illumination),  but  the  range  of  point  P  is  taken 

at  400  ft.  As  in  Figure  31,  ellipticity  is  now  a  function  of  c p  and  varies 

from  essentially  0  db  to  a  maximum  in  the  45-degree  planes.  Moreover, 

this  maximum  ellipticity  at  400  ft  exhibits  the  same  characteristics  as  in 

the  preceding  graph  (same  values  as  in  the  far-field,  within  the  beamwidth; 

much  lower  values  than  in  the  far-field,  for  elevation  angles  corresponding 

to  the  far-field  nulls). 
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Figure  34.  Focal  Feed  System:  Ellipticity  Related  to  Range 
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Figure  35.  Focal  Feed  System:  Ellipticity  Related  to  Range 
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Figure  36.  Focal  Feed  System:  Ellipticity  Related  to  Range 
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Figures  37  through  42  reflect  changes  (or  lack  thereof)  in  ellipticity 
when  the  paraboloidal  surface  is  perturbed. 

Figures  37,  38,  and  39  correspond  to  Figures  20,  21  and  22,  respec¬ 
tively,  where  randomly  distributed  imperfections  on  the  paraboloidal  surface 
are  simulated  by  letting  the  focal  length,  f  ,  vary  randomly  from  point  to 
point  according  to 


fi.j  (1  +  «»)  12.5  ft 

where 

!<ijl  <0.005 

Ellipticities  are  no  longer  exactly  0  db  on  boresight,  but  are  still 
extremely  small  (0.011  db,  0.  00G7  db,  and  0.0043  db,  respectively).  The 
patterns  themselves  are  essentially  the  same  as  for  the  unperturbed  cases 
(Figures  27,  28,  and  29),  except  that  they  present  a  small  spread  from  the 
corresponding  ideal  cases.  This  spread,  for  a  given  value  of  the  elevation 
angle,  0,  varies  randomly  as  a  function  of  <p.  Just  as  the  random  imper¬ 
fections  have  no  appreciable  effect  on  the  amplitude  patterns  (except  for  a 
substantial  loss  in  gain),  they  have  no  appreciable  effect  on  the  ellipticity 
patterns. 

Finally,  in  Figures  40,  41,  and  42,  the  focal  length  was  made  to 
vary  discontinuously  as  a  function  of  $  only,  as  in  Figures  23,  24,  and 
25.  The  same  table  (Table  III)  applies  to  these  cases.  Ellipticities  are  not 
0  db  on  boresight,  but  are  still  very  small  (0.  066  db,  0.  04  db,  and  0.  03  db, 
respectively).  However,  the  ellipticity  patterns  spread  considerably  on 
either  side  of  the  curves  pertaining  to  the  perfect  dish,  as  a  random  function 
of  cp.  This  type  of  semi-correlated  perturbation  is  far  more  disruptive  than 
the  random  perturbations  of  Figures  37  through  39,  even  though  their  peak 
amplitudes  are  the  same. 
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Figure  37.  Focal  Feed  System:  Ellipticity  Related  to  Paraboloid  Imperfections 
(Far-Field) 
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Figure  38.  Focal  Feed  System:  Ellipticity  Related  to  Paraboloid  Imperfections 
(Far-Field) 
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Figure  39.  Focal  Feed  System:  Ellipticity  Related  to  Paraboloid  Imperfections 
(Far-Field) 
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Figure  40.  Focal  Feed  System:  Ellipticity  Related  to  Paraboloid  Imperfections 
(Far-Field) 


113 


fO 

cr> 

O 

rg 


0  12  3 

ELEVATION,  Q  —  DEGREES 


< 


Figure  41 .  E 1 1 ipticity  Related  to  Paraboloid  Imperfections  (Far-Field) 
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Figure  42.  Ellipticity  Related  to  Paraboloid  Imperfections  (Far-Field) 


SHIFT  IN  LINEAR  POLARIZATION 


Figures  43  through  49  are  plots  of  the  deviation  angle  from  the  primary 
linear  polarization  to  the  polarization  at  the  observation  point,  P  ,  as  a 
function  of  the  azimuth  angle,  cp  ,  for  various  values  of  the  elevation  angle, 

0  .  They  are  plots  of  Equation  (60),  which  relates  the  far- fie  Id  cross- 
polarized  terms  to  the  Min-polarizedn  terms. 

As  discussed  in  Section  II  and  proved  in  Appendix  G,  the  cross-polarized 
terms  are  always  zero  in  the  H-  and  E-planes  at  the  observation  point,  P. 
Consequently,  regardless  of  whether  the  distribution  is  equibeam,  dipole¬ 
like,  or  pseudo-equibeam,  the  deviation  angle  is  zero  for  <p  =  0,  tt/2,  n , 
and  37r/2.  Since  the  cross-polarized  terms  arc  maximum  near  the  45-degree 
planes,  and  since  the  ratio  of  cross-polarized  to  in-polarized  terms  increases 
from  zero  to  a  maximum  value  in  the  region  of  the  amplitude  null,  the  devia¬ 
tion  angle  is  zero  on  boresight  and  increases  with  0  toward  the  amplitude 
null.  For  a  given  value  of  6  ,  it  is  maximum  in  or  near  the  45-degree  planes. 
Thus,  for  the  equibeam  case  (Figure  43),  the  deviation  angles  are  maximum 
in  the  45-degree  planes  as  a  result  of  the  integration  of  Equation  (33),  but 
the  amplitudes  are  rather  small  (one  degree  or  so),  even  around  the  null. 
However,  for  the  dipole-like  distribution  (Figures  44,  45.  46,  and  47),  the 
deviation  angles  are  maximum  near,  but  not  in,  the  45-degree  planes,  since 
amplitudes  in  the  E-plane  are  larger  than  amplitude  in  the  Il-plane.  Near  the 
nulls,  where  the  cross-polarized  terms  are  much  larger  than  the  "in- 
polarized”  terms,  the  deviation  angles  tend  toward  ±  tt/2.  Note  that,  past 
the  nulls,  the  deviation  angles  decrease  and  change  signs,  for  a  given  value 
of  the  azimuth  angle,  <p. 

Figure  48  applies  to  the  pseudo-equibeam  distribution.  The  curves 
have  essentially  the  same  character  as  for  the  dipole-likc  distribution, 
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Figure  43.  Focal  Feed  System:  Polarization  Deviation  Related  to  Primary 
Illumination  (Far-Field) 


except  that  the  deviation  angles  peak  up  in  the  45-degree  planes,  since  the 
amplitudes  in  the  H-  and  E-planes  are  equal. 

Figure  49  dramatizes  the  fact  that  for  long  focal  lengths  the  cross- 
polarized  terms  are  small,  and  consequently  the  deviation  angles  are  small. 

Those  plots  of  polarization  shifts  can  be  directly  related  to  the  ellip- 
ticity  plots:  no  deviation  on  boresight  (0  db  ellipticity);  largest  deviation  near 
the  nulls  (largest  ellipticity);  low  deviation  and  low  ellipticity  near  the  peak  of 
the  side  lobe. 

Effects  of  range  and  surface  perturbations  have  not  been  plotted, 
since,  for  the  cases  that  were  considered,  those  effects  can  be  inferred 
readily  from  the  amplitude  and  ellipticity  graphs. 
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Figure  44.  Focal  Feed  System:  Polarization  Deviation  Related  to  Primary 
Illumination  (Far-Field) 
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Figure  45.  Focal  Feed  System:  Polarization  Deviation  Related  to  Primary 
Illumination  (Far-Field) 
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Figure  46.  Focal  Feed  System:  Polarization  Deviation  Related  to  Primary 
Illumination  (Far-Field) 
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Figure  47.  Focal  Feed  System:  Polarization  Deviation  Related  to  Primary 
Illumination  (Far-Field) 
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Figure  48.  Focal  Feed  System:  Polarization  Deviation  Related  to  Primary 
Illumination  (Far-Field) 
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Figure  49.  Focal  Feed  System:  Polarization  Deviation  Related  to  Focal  Length 
(Far-Field) 


DEFOCUSING  IN  FOCAL  PLANE:  SUM  OF  TWO  BEAMS 


Figures  50  through  53  pertain  to  the  situation  of  the  paraboloidal  dish 
being  illuminated  by  two  feeds,  located  at  the  intersection  of  the  "vertical" 
plane  and  of  the  focal  plane,  each  one  a  distance  d  =  0. 192  ft  from  the 
boresight  axis,  on  opposite  sides  of  it  (See  Section  II  and  Figure  5).  Both 
feeds  emit  the  same  illumination.  The  plots  describe  the  resultant  far- 
field  beam  for  the  equibeam  and  dipole-like  illumination  functions. 

Figure  50  is  a  plot  of  amplitude  versus  elevation  angle,  0  ,  in  the 

1. 5 

far-field,  for  the  equibeam  distribution,  with  the  taper  function  cos  *  (n  -  rj), 
under  linear  polarization.  The  beam  does  not  exhibit  cylindrical  symmetry, 
but  is  broader  in  the  "vertical”  plane.  Cuts  through  the  vertical  and  hori¬ 
zontal  planes  are  plotted  in  the  graph.  For  the  sake  of  convenience,  the 
focal  feed  case  (Figure  9)  is  redrawn  (dotted  line).  The  change  of  beam- 
width  with  azimuth  is  due  to  the  particular  orientation  of  the  feeds.  The 
distance  d  was  "estimated"  for  a  crossover  angle  6  =  0.8  degrees 

(Section  II).  The  situation  around  the  "nulls"  (not  fully  plotted)  is  very 
complex,  but  of  little  significance. 

Figure  51  refers  to  the  same  situation  as  Figure  50,  except  that  the 
feed  illuminations  are  circularly  polarized.  Ellipticity  is  0  db  on  boresight 
(perfect  circularity)  and  increases  with  elevation,  but  very  slowly.  As 
expected,  ellipticity  is  not  independent  of  azimuth,  since  amplitude  is  not. 
Fundamentally,  the  reason  that  the  ellipticity  remains  so  low  is  that  the  only 
cross-coupling  mechanism  can  still  occur  only  through  the  z-directed  currents 
whose  contributions  in  the  far-field  are  negligible  (Equations  (23)  and  (24)) 
for  small  elevation  angles,  except  when  the  "dominant"  terms  sink  to  a  null. 
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Figure  50.  Defocusing  in  Focal  Plane:  Sum  of  Two  Beams  (Far-Field  Amplitude) 
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Figure  51 .  Defocusing  in  Focal  Plane:  Sum  of  Two  Beams  (Far-Field  Ellipticity) 


Figure  52  corresponds  to  Figure  50,  except  that  the  primary  illumi¬ 
nation  funetion  is  dipole-like.  Figure  13  has  been  redrawn  for  convenience 
of  comparison  with  the  single  foeal  fed  beam.  The  spread  between  vertical 
and  horizontal  patterns  is  not  as  large  as  in  Figure  50,  because  the  vertical 
plane  has  been  made  coincident  with  the  Il-plane.  To  some  extent,  the 
non-symmetry  of  eaeh  individual  beam  tends  to  cancel  eaeh  other.  Had  the 
vertical  plane  been  made  to  coincide  with  the  E- plane,  the  spread  between 
patterns  would  be  larger  than  in  Figure  50. 

Figure  53  refers  to  the  same  situation  as  Figure  52,  except  that  the 
feed  illuminations  are  circularly  polarized.  The  ellipticity  on  boresight  is 
rather  large  (0.41  db).  The  reason  is  that,  on  boresight,  the  extraneous 
components  X2  and  Y2  are  not  equal  to  eaeh  other  (Section  II),  sinee 
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Figure  52.  Defocusing  in  Focal  Plane:  Sum  of  Two  Beams  (Far-Field  Amplitude) 
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Figure  53.  Defocusing  in  Focal  Plane:  Sum  of  Two  Beams  (Far-Field  El  lipHcity) 
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In  the  "vertical"  plane  (plane  of  feeds),  the  ellipticity  keeps  on  rising 
slowly,  whereas  it  first  decreases  in  the  horizontal  plane  before  increasing 
rapidly  (the  minimum  ellipticity  is  not  far  from  0  ). 

In  summary,  the  type  of  primary  illumination  plays  a  fundamental  role 
in  determining  the  character  of  the  ellipticity.  Whereas,  in  the  circularly 
polarized  focal  feed  case,  ellipticity  is  always  0  db  on  boresight,  regardless 
of  the  distribution  over  the  dish,  this  is  no  longer  the  case  under  defocused 
conditions.  For  an  imperfect  illumination  function,  it  would  appear  that  the 
ellipticity  on  boresight  is  about  one  half  of  what  it  would  be,  for  the  focal  feed 
single  beam,  at  an  elevation  angle  corresponding  to  "cross-over",  that  is, 

e  . 

o 

DE FOCUSING  ALONG  AXIS 


Figures  54  through  57  pertain  to  the  situation  where  the  paraboloidal 
dish  is  illuminated  from  a  feed  that  is  displaced  a  distance  d  =  0.  417  ft 
along  the  boresight  axis,  in  an  attempt  to  "focus"  the  field  at  a  range 
r  =  400  ft.  (See  Section  II  and  Figures  6  and  7).  The  displacement  d 
was  computed  from  Equation  (96),  where  the  constant  1.  065  is  "estimated.  " 

Figure  54  contains  plots  of  amplitude  versus  elevation  for  the  equibeam 

1 . 5 

distribution,  with  the  taper  function  cos  '  ( n  -  rj),  under  linear  polarization. 
It  is  observed  that  the  defocused  amplitude  curve  (solid  curve)  for  a  range 
of  400  ft  is  very  close  to  the  far-field  amplitude  distribution,  although  it  is 
still  a  little  broader.  It  represents  a  dramatic  improvement  over  the  near¬ 
field  pattern  of  the  far-field  focused  beam.  This  slight  broadening  may  be 
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Figure  54.  Focal  Feed  System:  Amplitude  Related  to  Defocusing  Along  Axis 
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accounted  for  by  the  facts  that  the  displacement,  d  ,  is  "estimated”,  and  that 
even  for  the  best  estimate  of  d  ,  some  aberrations  will  always  be  present. 
The  gain  is  essentially  that  in  the  far-field,  whereas  the  gain  of  dot-dash 
curve  is  actually  down  1.  65  db. 

Figure  55  refers  to  the  same  situation  as  Figure  54,  except  that  the 
illumination  is  circularly  polarized.  The  defocused  ellipticity  at  400  ft 
closely  resembles  the  far-field  ellipticity.  Since  the  null  is  displaced  to 
a  larger  elevation  angle,  so  is  the  region  of  maximum  ellipticity. 

Figures  56  and  57  correspond,  respectively,  to  Figures  54  and  55, 
except  that  the  illumination  is  dipole-like.  Since  amplitude  varies  with 
azimuth,  the  amplitude  plots  refer  to  cuts  in  the  c p  =  45-degree  planes. 

The  same  characteristics  prevail  so  that  the  defocused  patterns,  at  400  ft, 
are  almost  exact  replicas  of  the  focused  patterns  at  infinity. 


Figure  55.  Focal  Feed  System:  Ellipticity  Related  to  Defocusing  Along  Axis 
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Figure  56.  Focal  Feed  System:  Amplitude  Related  to  Defocusing  Along  Axis 
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Mgure  57.  Focal  Feed  System:  Ellipticity  Related  to  Defocusing  Along  Axis 
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SECTION  V 


COMPUTED  RESULTS  -  CASSEGRAIN  SYSTEM 


AMPLITUDES 

Figures  58  and  59  are  plots  of  the  overall  field  amplitude  at  the  obser¬ 
vation  point  P  as  a  function  of  elevation  angle,  6  ,  when  the  feed  system  is 
linearly  polarized.  Those  plots  are  computed  values  of  Equation  (59) 
wherein  Equations  (158)  and  (161)  have  been  substituted.  The  actual  orien¬ 
tation  of  the  total  field  at  point  P  has  not  been  considered  separately,  but 
cnn  be  deduced  from  the  ellipticity  curves  (Figures  62  through  65). 

Parameters  applicable  to  those  graphs  are  summarized  in  Table  IV. 


Table  IV 

Constants  for  the  MITRE  Antenna 


FREQUENCY  =  1288  Me 

Taper  function  h(0)  =  cos 

2(2  6) 

(12  db  taper) 

Sub-reflector 

Aperture  half-angle 

6 

o 

tl 

CO 

o 

o 

Eccentricity 

e 

=  2.45 

Parameter 

a 

=  1.063  ft 

Main  dish 

Aperture  angle 

\ 

=  118° 

Focal  length 

f 

=  12.  5  ft 

Shadow  aperture  half-angle 

\ 

o 

[> 

d 

i— i 

ll 

Figure  58  compares  the  far-field  amplitudes  when  the  primary  illumi¬ 
nation  is  equibeam  (solid  line)  and  dipole- like  (dashed  line),  with  the  same 
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taper  function  h (6)  =  cos  (26).  It  can  be  observed  that  both  curves  coalesce, 
except  near  the  null.  The  reason  is  that  the  half-angle  subtended  by  the  sub¬ 
dish  is  relatively  small  (0  =  30  degrees)  so  that,  as  noted  at  the  end  of 

Appendix  D,  the  two  distributions  are  about  the  same.  Thus,  for  the  dipole¬ 
like  distribution,  the  difference  in  field  strength  at  the  edge  of  the  subdish 
between  H-plane  and  E-plane  is  only  1.  25  db. 

The  common  3-db  half-beamwidth  is  0.87  degrees. 

The  amplitude  pattern  is  constant  with  azimuth  for  the  equibeam 
distribution.  For  the  dipole-like  distribution,  it  is  essentially  constant, 
except  near  the  null. 

Figure  59  considers  far-field  amplitude  changes  induced  by  the  shadow¬ 
ing  effects  of  the  subreflector  and/or  feed  horn.  Since,  in  the  particular 
case  under  study,  the  two  effects  are  approximately  overlapping,  the  contri¬ 
butions  from  the  main  dish  within  an  angle  77^  =  10.  7  degrees  have  simply 
been  omitted.  This  is  equivalent  to  adding  the  negative  field  of  a  paraboloid 
of  the  same  focal  length,  but  of  aperture  angle  =  (180  -  10.  7)  degrees. 

The  resultant  field  would  be  low  in  amplitude,  but  very  broad.  Thus,  the 
overall  field  is  slightly  weaker  on  boresight  (0.  36  db  loss),  slightly  sharper 
since  the  subtracting  field  is  essentially  constant  near  boresight  (new  3-db 
half-beamwidth  =  0.  84  degrees),  with  a  markedly  deeper  null  and  higher 
sidebands. 

Figure  60  also  considers  the  shadowing  effect  problem,  but  whereas 
the  previous  graph  applies  to  the  equibeam  distribution,  this  one  applies  to 
the  dipole-like  distribution.  The  characters  of  the  two  graphs  are  essen¬ 
tially  the  same. 
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Figure'58.  Cassegrain  System:  Amplitude  Related  to  Primary  Distribution 
(Far-Field) 
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Figure  59.  Cassegrain  System:  Amplitude  Related  to  Shadowing  Effects 
(Far-Field) 
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Figure  60.  Cassegrain  System:  Amplitude  Related  to  Shadowing  Effects 
(Far-Field) 
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Figure  61  relates  far-field  amplitude  to  near-field  amplitude  (400  ft 
range)  for  the  equibeam  distribution,  with  shadow  effects  taken  into  account. 

Remarks  pertaining  to  the  focal  feed  system  apply  here  too.  The 
amplitude  distributions  are  practically  the  same  near  the  peak  of  the  beam, 
although  the  near-field  3-db  half  beamwidth  is  increased  to  0.90  degrees  at 
a  range  of  400  ft.  The  relative  gain  at  400  ft  is  down  1.  67  db  from  the  far- 
field  gain,  as  referred  to  a  standard  horn.  Major  differences  appear  in  the 
null  region  where,  at  400  ft,  the  null  is  "filled-in"  and  the  new  null  is  pushed 
to  a  larger  elevation  angle. 

ELLIPTICITIES 

Figures  62  through  65  are  plots  of  ellipticities  at  the  observation  point, 
P  ,  as  a  function  of  the  elevation  angle,  6  ,  when  the  feed  system  is  circu¬ 
larly  polarized.  Those  plots  are  computed  values  of  Equation  (56)  wherein 
Equations  (158)  through  (161)  have  been  substituted.  The  conditions  of 
Figures  62  through  65  parallel,  respectively,  conditions  of  58  through  61. 

Figure  62  shows  the  somewhat  better  performance  of  the  equibeam 
distribution  over  the  dipole-like  distribution  (far-field).  However,  the 
difference  in  ellipticity  is  very  small,  except  in  the  null  region.  Ellipticities 
at  the  3-db  half-beamwidth:  0.  14  db  for  the  equibeam  distribution;  0.  17  for 
the  dipole-like  distribution.  As  discussed  previously,  at  the  edge  of  the 
subreflector,  the  dipole-likc  distribution  shows  a  difference  of  only  1.  25  db 
in  amplitude  between  the  H-  and  E-planes,  that  is,  an  ellipticity  of  1.25  db. 
On  the  other  hand,  the  equibeam  distribution  does  not  produce  the  perfect 
circularity  exhibited  in  the  focal  feed  system  because  of  the  existence  of 
cross-polarized  currents  on  the  hyperboloidal  subreflector. 
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Figure  61.  Cassegrain  System:  Amplitude  Related  to  Range 
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Figure  62.  Cassegrain  System:  Eliipticity  Related  to  Primary  Distribution 


Figure  63  considers  far-field  eliipticity  changes  induced  by  taking 
shadowing  effects  into  account  (equibeam  distribution).  The  shadow  tends 
to  increase  the  eliipticity,  but  the  effect  is  very  small,  except  in  the  null 
region.  Since  the  amplitude  curve  with  shadow  effects  is  slightly  narrower, 
the  ellipticities  at  the  3-db  half-beamwidths  are  essentially  the  same  for 
shadow  and  non-shadow  conditions,  that  is,  0.  14  db. 

Figure  64  deals  with  the  same  situation  as  Figure  63  except  that  the 
primary  distribution  is  dipole-like.  The  same  remarks  apply.  The  ellipti¬ 
cities  at  the  3-db  half-beamwidths  are  essentially  the  same  at  0.  17  db. 

Figure  65  relates  far-field  ellipticites  to  near-field  ellipticities  (400-ft 
range)  for  the  equibeam  distribution  (shadow  effects  taken  into  account). 
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Figure  63.  Cassegrain  System:  Ellipticity  Related  to  Shadowing  Effects 
(Far-Field) 
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Figure  64.  Cassegrain  System:  Ellipticity  Related  to  Shadowing  Effects 
(Far-Field) 
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Figure  65.  Cassegrain  System:  Ellipticity  Related  to  Range 


Near-field  ellipticity  peaks  up  to  a  low  value  in  the  corresponding  near¬ 
field  minimum  amplitude  region.  Although  this  effect  tends  to  lower  the 
ellipticity  curve,  the  ellipticities  at  the  3-db  half-beamwidth  remain  essen¬ 
tially  unchanged  from  far-field  to  near-field  (400  ft)  at  0.  14  db. 

PARABOLOID  ILLUMINATION 

Plots  of  the  paraboloid  illumination  from  the  subreflector  are  contained 
in  Figures  66  and  67. 

Figures  66  considers  H  Equation  (153),  the  solid  lines,  and  H 

HI  El 

Equation  (154),  the  dotted  line,  as  functions  of  the  elevation  angle,  77  ,  for 

2 

the  equibeam  primary  distribution  (cos  20  taper).  Thus,  H  is  the  field 

HI 

component  in  the  Il-plane  (ip  =  0,  7 r),  at  the  paraboloidal  dish,  along  the 

Tj-direction.  The  H  component  is  zero.  The  HIir>  component  (radial 
HZ  H3 
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Figure  66.  Cassegrain  System:  Paraboloid  Illumination 
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component)  is  not  zero,  but,  on  the  average,  40  db  below  the  peak  of  H  . 

HI 

It  does  go  up  to  33  db  below  the  peak  of  H  over  small  regions.  H  is 

HI  H3 

not  plotted. 

H  ,  the  dotted  line,  is  the  field  component  in  the  E-plane  (d>  =  ±  ?r/2), 
E 1 

directed  along  the  x-direction.  H  and  H  are  zero. 

E2  E3 

The  and  H  components  are  matched  very  closely  to  each 

other,  but  not  exactly,  although  the  primary  illumination  is  equibeam.  The 
difference  is  accentuated  at  the  edge  of  the  dish  where  it  amounts  to  about 
1.4  db.  The  curves,  as  a  whole,  exhibit  shallow'  ripples  which  would  become 
accentuated  in  the  far-field  (note  that  the  curves  refer  to  the  fields  at  the 
paraboloid). 

The  phase  of  H  -  relative  to  the  focal  feed  phase  -  is  also  plotted. 

HI 

It  oscillates  about  ±9  degrees  about  the  value  of  -3  degrees.  This  average 
value  depends  on  the  a- parameter  and  has  no  special  significance. 

Figure  67  applies  to  the  dipole-like  distribution.  The  matching  between 

H  and  H  ,  although  good,  is  not  as  close  as  in  the  preceding  case.  In 
HI  El 

particular,  the  mismatch  at  the  edge  amounts  to  about  2  db.  Those  results 

are  consistent  with  computed  ellipticities  (Figure  62).  Again,  H  is  zero; 

H  is  very  small  (on  the  average,  40  db  below  the  peak  of  H  );  and  Hp,0 
Ho  nl  La 

and  H  are  zero. 

Il<  O 

This  time,  the  phase  II  relative  to  the  focal  feed  phase,  is  plotted 

E 1 

versus  17  .  The  phase  oscillation  amplitude  is  slightly  larger  than  in  the 
preceding  case. 

Figures  66  and  67  can  be  used  to  judge  the  amount  of  approximation 

resulting  from  treating  the  Cassegrain  feed  system  by  ray  optics:  in  the 

latter  case,  H  and  H  would  be  perfectly  matched;  H  would  be 
H 1  E 1  11 3 
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Figure  67.  Cassegrain  System:  Paraboloid  Illumination 
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zero;  the  phases  (referred  to  the  focal  feed  system)  would  be  constant;  there 
would  be  no  amplitude  ripple;  the  tapers  would  be  those  of  the  primary 
illumination  function  (except  for  range  dependency);  and  in  summary,  the 
secondary  illumination  would  be  a  replica  of  the  primary  illumination, 
except  for  angular  magnification.  On  all  counts,  ray  optics  would  give 
quantitatively  inaccurate,  but  qualitatively  valuable  information. 
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APPENDIX  A 


EXPRESSION  FOR  n  IN  CARTESIAN  COORDINATES  (PARABOLOID) 

The  coordinate  system  used  is  shown  in  Figures  3  and  4.  The  origin 
is  at  point  0,  the  focus  of  the  paraboloid,  and  the  z-axis  is  oriented  along 
the  axis  of  symmetry.  Let  be  a  point  on  the  paraboloid,  located  in  the 
x-z  plane,  such  that  the  angle  of  zOQ^  is  V  •  Consider  n^  at  Q  (  .  The 
geometric  properties  of  the  parabola  are  such  that  a  ray,  transmitted  from 
point  0  and  incident  at  ,  will  be  reflected  parallel  to  the  z-axis.  Con¬ 
sequently,  the  angle  of  n^  and  the  z-axis  is  (7 r  -  r\)/2  .  The  components 

of  n  at  Q  are 
o  o 


n 


\ 


n 


z 


(n  - 


(n  - 


—  cos  r/  2 


sin  7]  2 


(166) 

(167) 


The  components  of  n  at  any  point  Q  are  obtained  by  rotating 
about  the  z-axis  through  an  angle  ip  .  Analytically,  the  results  will  be 
obtained  by  operating  on  Equations  (166)  and  (167)  by  a  rotation  matrix 


% 

COS  ifj 

—  sin  0 

0 

—  COS  77/2 

—  COS  0  COS  77/ 2 

\ 

sin  0 

cos  0 

0 

0 

-- 

—  sin  0  cos  ij  2 

nz 

0 

0 

1 

sin  77  2 

sin  77  2 

(168) 
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APPENDIX  B 


EXPRESSION  FOR  p  IN  CARTESIAN  COORDINATES  (PARABOLOID) 

Just  as  in  Appendix  A,  the  coordinate  system  used  is  shown  in  Figures 

3  and  4.  Let  p^  correspond  to  point  in  the  x-z  plane.  The  components 

of  p  are 
o 


p  sin  77 

(170) 

p  cos  77 

(171) 

The  expression  for  p  at  point  Q  is  obtained  by  operating  on  Equations 
(170)  and  (171)  with  a  rotation  matrix  (Appendix  A), 


COS  0 

—sin  0 

0 

sin  77 

sin  77  cos  0 

sin  0 

COS  0 

0 

P 

0 

p 

sin  77  sin  0 

(172) 

0 

0 

1 

COS  77 

cos  77 

Since  the  equation  of  the  parabola  in  polar  form  is 

p  =  f  (sin^  rf/2) 


(173) 


it  follows  that 


sin  rj  cos  0 


sin^  77/2 


sin  77  sin  0 


(174) 


COS  77 
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It  is  noteworthy  that  p  plus  the  distance  from  Q  to  the  focal  plane 


(z  =  0)  is  a  constant  equal  to  2f: 


p  +  p  cos  r)  =  p  (1  +  cos  r\) 


(175) 


- ^ -  (2  sin2  t,/ 2)  2f 

sin“  r)/2 

so  that  any  ray  emanating  from  the  focus  loses  the  phase  2  kf  to  the 
focal  plane. 

If  77^  is  the  value  of  77  at  the  edge  of  the  parabola,  the  aperture  D 
can  be  expressed  in  terms  of  and  f  according  to 

D  -  2p  sin  r)o  -  4f  cot  r]n/2  (176) 
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APPENDIX  C 


DIFFERENTIAL  AREA,  dA,  OF  THE  PARABOLOID 

The  differential  element  of  area,  dA  ,  produced  by  a  differential 
change  dr)  in  r\  and  d^  in  j p  is  (Figure  68) 


d  A 


p  sin  77 


d  \b 


-i  U 


>2  *  (t)' 


At) 


But 


sin~  77/2 


which  is  the  equation  of  the  generating  parabola  in  polar  form.  Hence, 


dp  ^  r  •  — 3  m  o  1  f  cos  77/2 

-7-^—  =  —  2  f  sin  3  77  2  cos  77/2  •  —  —  — L — 


,  -  - 1  sin  77/ z  cos  77/-  0 

d  77  z 


sin2  77/2 


so  that 


(*)■ 


ff _  +  f2  cos2  77/2 


sin  ^  77/ 2  sin^  t]/2 


f  (1  +  cot2  77/2) -  f 


sin“  77/2 


sin2  77/2 


It  follows  that 


d  A 


f~  sin  77 
sinf‘  77/2 


d  77  d0 


(177) 
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Figure  68.  Differential  Area 


APPENDIX  D 


EXPRESSIONS  FOR  PRIMARY  PATTERNS  IN  CARTESIAN  COORDINATES 
EQUIBEAM  CONDITIONS 

— ►  — ►  _ 

In  deriving  an  analytic  expression  for  the  incident  field  H.  =  H.^e  J 
under  equibeam  conditions,  the  following  properties  pertaining  to  the  feed 
horn  and  its  primary  pattern  are  postulated  (those  properties  may  also  be 
used  to  define  the  equibeam  conditions): 

(1)  The  incident  wave  is  spherical  and  issues  from  a  well-defined 
’’phase  center.  " 

(2)  The  H-field  is  perpendicular  to  the  direction  of  propagation,  that 
is,  perpendicular  to  vector  p  (Figure  2). 

(3)  The  ’’vertical”  component  of  H  ,  H^T  ,  has  no  y-component 

(that  is,  no  ’’horizontal”  component)  in  the  H-  and  E-planes 

(ip  =  0,  7T  and  ip  =  ±-  ,  respectively).  The  ’’horizontal” 

component  of  H.  ,  II  ,  present  under  circular  polarization, 
io  II 

has  no  x-component  (that  is,  no  ’’vertical”  component)  in  the 
H-  and  E-planes. 

(4)  Each  linear  component  of  H  is  equibeam  in  the  sense  that 
the  strength  of  H  in  the  H-plane  is  the  same  as  the  strength 
of  H  in  the  E-plane. 

Finally,  it  is  not  clear,  at  first,  what  directions  must  be 

(5)  assigned  to  the  H-field  in  the  other  azimuthal  planes.  It  turns 
out  that  the  supplementary  constraint  is  that  orthogonal  linear 
components  in  the  feed  give  rise  to  orthogonal  linear  components 
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at  any  point  on  the  spherical  wavefront.  Were  it  not  so,  the 
structure  would  exhibit  "preferred"  planes  which  would  be 
inconsistent  with  the  underlying  goal  of  axial  symmetry. 

Thus,  in  summary,  equibeam  conditions  are  such  that  the  field  possesses 
axial  symmetry  in  strength,  but  non-axial  symmetry  in  direction.  Rather, 
directions  are  constrained  in  such  a  manner  that  orthogonal  families  at  the 
feed  give  rise  to  orthogonal  families  on  the  spherical  wavefront. 


Postulates  (1)  and  (2)  will  be  satisfied  if,  in  the  coordinate  system 
77  -  ip  -  p  ,  H  and  are  expressed,  respectively,  as 


/ 

p 

P 

<1 

;  HH  - 

<1 

0 

0 

(178) 


Postulate  (4)  requires  that 

9  9 

p2  +  q“  =  1 

p'2+q'2=l  (179) 

whereas  Postulate  (5)  requires  the  equation 


which  is  satisfied  by 


pp'  *  qq'  0 


p'  =  +  q 


(180) 


q'  =  f  p 


(181) 
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The  simplest  solution  to  Postulate  (3),  consistent  with  Equations  (179)  is 


p  +  COS  l/j 

q  -  +  sin  0  (182) 

where  the  minus  sij^n  corresponds  to  a  positive  value  of  the  x-component  of 
— » 

H.  ,  the  positive  sign  to  a  negative  value  of  the  x-component. 

To  get  an  expression  for  H^T  in  the  Cartesian  coordinate  system,  one 
must  first  rotate  the  77  -  0  -  p  system  about  the  local  0-axis  through 
a  -  ij  angle,  then  -  about  the  z-axis  -  through  a  -0  angle  (Appendix  E). 
The  transformation  matrices  are,  respectively, 


[T] 


—1 


cos  77  0  sin  77 

0  1  0 
—  sin  77  0  cos  77 


(183) 


and 


IT] 


-t b 


cos  0  —  sin  0  0 

sin  \b  cos  0  0 

0  0  1 


(184) 


The  expression  for  in  Cartesian  coordinates  becomes 


—  cos  i /; 

9  0 

—  cos^  0  cos  77  +  sin  0 

II 

r~ 

H 

1 

-CJ 

—  sin  0 

- 

—  cos  0  sin  0  (1  4-  cos  77) 

0 

sin  77  cos  0 
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(185) 


The  projection  of  H.^  onto  the  x-y  plane  is  plotted  semi- graphically  in 
Figure  (69). 

From  Equation  (101),  and  by  using  the  same  procedure,  one  obtains 
for  Hh 


—  sin  p 

—  cos  ip  sin  ip  (1  \  cos  q) 

»h  -  m_*  m_. 

cos  \p 

- 

.  9  ,  2  • 

—  sin~  p  cos  q  -*■  cos“  u 

0 

sin  q  sin  ip 

(186) 


It  may  be  readily  verified  that 


«V  '  «H  0 


(187) 


and  that 

'  nx  | 2  =  |Hh|2  =  1  (188) 

regardless  of  ip  and  rj  . 

As  a  direct  consequence,  if  H  and  H  (as  defined  above)  are  applied 

V  H 

simultaneously  and  in  phase-quadrature,  each  point  on  the  spherical  wavefront 
will  be  circularly  polarized. 


Under  conditions  of  equibeam  circular  polarization,  the  complete 
expression  for  becomes 


=  h  (q)  sin~  q/ 2 


2  9 

(—  cos~  ip  cos  q  +  sin“  ip)  cos  cot 

4-  (—  cos  ip  sin  ip  [l  \  cos  q\)  sin  ot 

(—  cos  ip  sin  ip  [l  +  cos  q])  cos  ojt 
O  9 

+  (  —  sin~  ip  cos  q  +  cos“  ip)  sin  cot 
sin  7/  cos  ip  cos  a>t  +  sin  q  sin  ip  sin  cot 


(189) 
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Figure  69.  Equibeam  Distribution:  Projection  of  H-Field  onto  x-y  Plane 
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where  sin^p/2  is  a  common  factor  reflecting  the  strength  dependence  on 
2 

range  (l/p  cc  sin  p/2),  and  h(p)  is  a  common  tapering  function  dependent 
on  p  only. 

It  has  been  found  that 


h  (77)  cos^’°  ( n  -  7/)  (190) 

fits  the  amplitude  pattern  very  closely  for  an  equibeam  system  with  a  10  db 
taper  at  the  edge  of  the  dish  (p  =  118  degrees). 

DIPOLE- LIKE  DISTRIBUTION 

The  dipole-like  distribution  is  considerably  simpler  than  the  equibeam 
distribution.  Postulates  (1)  and  (2)  still  apply,  but  Postulate  (3)  is  modified 
to  require  that  the  "vertical"  component  of  H  ,  I1^T  ,  contain  no  y-eomponent 
anj'where  on  the  spherical  wavefront,  and  that  the  "horizontal"  component  of 
H  ,  ,  contain  no  x-component  anywhere  on  the  spherical  wavefront. 

By  Postulates  (1)  and  (3),  H  and  H  are  expressed,  respectively, 

v  H 

as 


m 

0 

> 

"v 

0 

hh 

m ' 

n 

n  ' 

(191) 


Postulate  (2)  expresses  the  conditions  that 

H  y  •  p  =  0 

flH  .  p  =  0  (192) 
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But 


sin  7]  cos  \p 


P  P 


sin  7]  sin  \Jj 


cos  rj 


bv  Equation  (172).  Hence,  it  follows  that 


—  COS  7] 

0 

sin  7]  cos  \p 


(193) 


and 


0 

—  COS  7j 

sin  7j  sin  ip 


(194) 


Although  the  projections  of  H  and  H  on  the  x-y  plane  are  equal 

v  H 

and  orthogonal,  it  is  apparent,  from  Equations  (193)  and  (194),  that 


'V  •  nll 


Hii  ^  0 


(195) 


that  is,  and  are  not  orthogonal,  except  when  if 

on  the  H-  and  E-planes.  Also 


=  0  ,  ~  or  ±  7r  /2 


Illy!  fc  of  l // ;  jHj||  =  fc  of  (// 


(196) 
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In  particular,  when  17  =  120  degrees,  |ll^T|  Is  larger  in  the  H-plane 
(ip  =  0)  than  in  the  E-plane  (i p  =  ±7r/2)  by  6  db. 

When  H^r  and  II ^  are  applied  simultaneously  and  in  phase  quadrature 
(circularly  polarized  feed),  H.  becomes 


—  COS  7/  COS  cot 


Hjo  h  (7/)  sin^  7//2 


—  cos  r]  sin  cot 
sin  77  cos  (cot  —  1 //) 


(197) 


where  sin  77/2  is  the  distance  factor  and  h (77)  is  a  tapering  function  inde¬ 
pendent  of  azimuth. 


In  order  to  determine  the  ellipticity  on  the  spherical  wavefront,  when 
the  feed  is  circularly  polarized,  one  can  use  the  coordinate  transformation 
method  of  Appendix  E  and  the  appropriate  equations  of  Appendix  F  to  compute 
the  ellipticity  of  H  on  the  dish.  Replacing  cos  cot  by  1  and  sin  cot  by  j  , 
and  neglecting  the  common  77-dependency  terms,  one  re-writes  Equation  (197) 
as 


—  COS  7] 


H 


io 


OC 


—  j  cos  rj 

sin  7/  (cos  \b  +  j  sin  iff) 


(198) 


By  the  applications  of  Equation  (210)  (but  with  8  replaced  by  (77  -  7r), 
<p  by  ip )  ,  one  obtains 
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H 


io  transf. 


—  cos  rj  (cos  0  [—cos?]]  +  sin  i!j  [—  j  cos  yj\  f  sin  “  rj  (cos  lb  +  j  sin  </v) 

—  sin  ifj  (--  cos  rj)  +  cos  i //  (—  j  cos  if) 

sin  77  [cos  iIj  (  —  cos  77)  +  sin  \Jj  (— j  cos  rf)\  —  cos  77  [sin  if  (cos  1./  +  j  sin  i'j) \ 


cos  i!j  +  j  sin  yj 


los  77  (sin  <jr  —  j  cos  0) 


(199) 


Equation  (199)  represents  the  total  field  on  the  spherical  wavefront  in 
a  right-hand  coordinate  system  oriented  along  77  ,  ip  ,  and  -p  .  As  expected, 
the  p-eomponent  is  zero.  By  Equations  (217)  and  (218), 


2  |SM 

2 

/  2 
(COS 

9  .  9 

4-  cos4*  7/  sin“  i/r)  + 

/  .  2  /  2 
(sin  1/7  V  COS  T] 

cy 

cos'-  iJj)  — 

,2 

/  2 
(cos 

>  2 

(  2  f  2 

(sin  i//  +  cos  77 

2 

2  \S~ 

1  = 

</> 

+  cos“  if  sin~  + 

cos  yj)  ♦ 

oi¬ 

ls4  “^=-(1  -  COS  Tf) 

'S~|  =  (1  4  cos  7 7)  (20°) 

\  2 

so  that 

Elliptic!  ty  -  - - - 

!cos  I/I  (201 ) 

It  is  remarkable  that  the  elliptieity  is  independent  of  azimuth.  The 
elliptieity  goes  to  1  when  77  =  n  (on-axis)  and  to  infinity  (linear  polariza¬ 
tion)  when  77  =  tt/2  . 
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By  application  of  Equation  (221)  to  Equation  (199),  it  is  found  that  the 
orientation  of  the  major  axis  is  at  zero  degree  in  the  local  coordinate  system. 
Consequently,  in  the  fixed  x-y-z  coordinate  system,  the  orientation  of  the 
ellipses  is  actually  radial. 

OTHER  DISTRIBUTIONS 

One  may  attempt  to  remedy  the  shortcomings  of  the  dipole-like  distri¬ 
bution  by  requiring  that  H  and  H  be  constant  with  azimuth.  This  will 

V  H 

be  the  case  for  H^,  if 

jyc  (—x  cos  77  +  z  sin  77  cos  i/»)  |  =  1 

(202) 

regardless  of  ^  .  A  solution  for  this  equation  is 

2  2  2  _  *2 

yc  (cosz  i)  +  sin"  77  cos'1  c ]/) 

(203) 

In  the  same  fashion,  it  will  be  found  that  II  is  constant  with  >p  if 

H 


|ys  (—  x  cos  77  +  7  sin  77  sin  </0  |  =  1 

(204) 


whose  solution  is 


=  ( 


1  -2  .  2  ,  \  -  % 
cos  77  +  sin  77  sin  xp ) 


(205) 


Thus,  the  expression  for 
polarized  feed  system  is 


H. 

10 


in  Cartesian  coordinates  for  a  circularly 
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—  cos  77  yc  cos  a) t 


—  cos  77  ys  sin  6jt 


(sin  7/  cos  lA  y  )  cos  cot  -1  (sin  77  sin  <A  ys)  sin  a>t 


(206) 


A  few  remarks  are  in  order.  Although  I1^T  and  are  equibeam  in 

amplitude  and  can  be  applied  in  phase-quadrature,  the  resultant  on  the 

spherieal  wavefront  will  not  be  circularly  polarized  everywhere,  as  in  the 

ease  of  the  true  equibeam  distribution.  Furthermore,  in  contradistinction 

with  the  dipole-like  distribution,  the  ellipticity  is  not  constant  with  azimuth 

because  the  spaee  angle  of  and  H  varies  with  4  (at  constant  77). 

Only  in  the  vertical  and  horizontal  planes  are  H  and  H  orthogonal  and 

V  H 

H.  circularly  polarized. 


10 


This  phenomenon  ean  be  formally  evidenced  by  applying  the  coordinate 
transformation  method  of  Appendix  E  and  the  appropriate  equations  of 
Appendix  F.  It  is  found  that  the  magnitudes  of  the  circular  waves  are 


s+  *  -  yc  }'s  cos  n 


S  |  1  +  yc  ys  cos  77 


(207) 


that  is,  they  are  functions  of  ip  >  for  a  P^ven  value  of  77  . 

Because  the  ellipticity  is  not  constant  with  azimuth,  this  pseudo- 
equibeam  distribution  cannot  be  produced  by  a  cylindrically  symmetrical 
feed  such  as  a  cylindrically  symmetrical  circular  horn.  Another  point  of 
view  is  that  the  orthogonal  and  phase-quadratured  sources  of  this  pseudo- 
equibeam  distribution  cannot  be  considered  infinitely  small  in  extent  and 
are  not  equivalent  to  a  single,  statically  charged,  rotating  source. 
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Finally,  note  that  when  the  elevation  angle  is  small  (near  on-axis 
region  with  cos  77  «  1,  sin  77  «  0),  all  three  distributions  are  essentially 
the  same  on  the  spherical  wavefront.  Marked  differences  show  up  as  the 
elevation  angle  opens  up. 
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APPENDIX  E 


COORDINATE  TRANSFORMATIONS  AT  OBSERVATION  POINT  P 


For  small  values  of  j  0  ,  the  x  and  y  coordinates,  at  the  observation 

point  P  ,  are  essentially  perpendicular  to  the  line  of  sight. 

If  the  value  of  |  o|  is  not  small,  it  is  convenient  to  deal  with  a  local 
coordinate  system,  x"  ,  y”  ,  z"  ,  such  that  z"  is  along  the  line  of  sight, 
for  purposes  of  computing  amplitudes  and  polarizations. 


Let  H 


S 


be  (X  ,  Y  ,  Z)  ,  as  given  by  Equations  (52)  and  (53) 


(Z  can  be  obtained  from  Equations  (21)  or  (28),  although  this  is  not  essential 
Equation  (211). 


First,  rotate  the  coordinate  system  about  the  z-axis  until  the  vector 
a^  lies  in  the  z-x'  plane  (where  x'  is  the  rotated  x-axis),  that  is,  perform 
a  evrotation  (Figure  70).  This  is  accomplished  by  the  rotation  matrix. 


IT], 


cos  <£>  sin  4>  0 

—  sin  0  cos  0  0 

0  0  1 


(208) 


Next,  perform  a  rotation  about  the  y'-axis  (the  rotated  y-axis),  in  the 
direction  from  z  to  x'  .  This  is  accomplished  by  the  rotation  matrix 


cos  6  0  —  sin  0 


[T1 


0  1  0 

sin  6  0  cos  0 


(209) 
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Figure  70.  Rotations  Related  to  Lirte-of-Sight  Vector  a. 

(1)  Rotation  about  z-axis 

(2)  Rotation  about  y'-axis 
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No  more  rotation  appears  to  be  necessary.  The  new  z-axis,  z"  ,  is  now  lined 
up  with  a^  ,  the  direction  from  the  observation  point,  P,  to  the  focus. 

In  the  new  coordinate  system,  the  field  components  become 


X" 

\ 

cos  X  +  sin  cf)  \ 

Y  " 

Z" 

II 

O' 

Y 

Z 

»  ITlfl 

—  sin  <3  X  +  cos  0  ^ 

z 

cos  0  (cos  <f>  X  +  sin  <*>  \  )  —  sin  0  Z 
-  —  sin  <f>  X  -f 

sin  0  (cos  <3  X  f  sir 

Since  the  z"-component  must  be  0,  it  follows  that 

Z  -  —  \  (cos  <3  X  +  sin  <3  Y ) 

COS  V 

and  that 

X" 

\  " 

Z" 

(212) 


“  (X  cos  <3  +  Y  sin  0) 
(—  X  sin  <3  +  Y  cos  <}>) 


cos  4>  ^ 


n  <3  \  )  cos  0  Z 


(210) 


(211) 
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APPENDIX  F 


COMPUTATION  OF  ELLIPTICITY 

Given  an  x-component,  X  =  Ae^a  ,  and  a  y-component,  Y  =  B  , 
where  the  magnitudes  A  and  B  are  real  and  arbitrary,  and  the  phases  a 
and  $  are  real  and  arbitrary,  what  is  the  elliptic ity  (ratio  of  major  to 
minor  axis)  ? 

}a 

Write  A  e  in  vector  form  as 


(213) 


that  is,  as  the  sum  of  two  circularly  polarized  waves: 


A  eja 
2 


a  right-hand  wave  advanced  by  a,  “ —  a 


and 


A  eja 
2 


,  a  left-hand  wave  advanced  (in  its  own  direction) 


by  a  ,  — eJ a 


In  the  same  fashion,  the  y-component  is 


0 

J_ 

1 

Bei^ 
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(214) 


that  is,  the  sum  of 
1 


B 


J 


(f>  *  t ) 


-J 


.  |  *  ,  ,  j  .  /  j  n  \  B  '  j  +  2  ) 

,  a  right-hand  wave,  advanced  by  (p  ~  1  ,  ~  e 


and 


j  (p  -  f ) 


,  a  left-hand  wave,  advanced  by  ^ ft  —  ,  ~ 


f  (#•  -  f ) 


(Note  how  the  right-hand  wave  is  advanced  by  tt/2  ,  whereas  the  left-hand 
wave  is  retarded  by  tt/2)  . 

The  right-hand  waves  and  the  left-hand  waves  can  be  respectively 
added  up  together,  that  is 


S+  =  T 


A+  eJQ  +  B+ 


)  (^i) 


ilUj  Y] 


S-  = 


A-  eJ a  +  B 


-j  (fi-i) 


(215) 


j  t\  -  j  Y J 

The  magnitudes  of  S+  and  S  are,  respectively, 


(216) 


|S+|  j^(A  cos  a  —  B  sin  f})"  +  (A  sin  a  •  B  cos 

=  4-  [(Re  X  -  l'«Y)2  +  dm  \  +  He  Y)2] 


js- 


j"  (A  cos  a  +  B  sin  fi)~  +  (A  sin  a  —  B  cos  /j)-  J 
r (Re  X  f  Ini  \)2  +  (In.  \  -  Re  A  )2  J 


(217) 


(218) 
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The  ellipticity  is  obtained  by  using  those  equations  in 


Kllipticity 


1S+1  +  |S-j 
IS+I  -  |S-| 


(219) 


The  phase  of  S  with  respect  to  the  datum  direction  (x-axis)  is 


arg  S+  -  tan 


.  ]  /sin  a  +  cos  \  /  B  Ini  X  t  A  Re  ^  \ 

\  cos  a  —  sin  fl  )  *  \  B  Re  X  —  A  Ini  V  / 


(220) 


and,  in  the  same  fashion,  the  phase  of  S  is 


arg 


tan 


-1  /  sin  a  —  cos  ft  \  ^  —  ]  /  B  Im  X  —  \  Re  X  \ 

\  cos  a  4  sin  ft  )  \  B  Re  X  +  A  Ini  X  / 


(221) 

Note  that  if  arg  S  is  positive,  it  denotes  a  right-hand  angle;  if 
arg  S  is  positive,  it  denotes  a  left-hand  angle. 

The  major  axis  corresponds  to  the  situation  when  S  and  S  are  in 
phase,  that  is,  it  lies  half-way  between  arg  S+  and  arg  S  or 


Orientation  of  major  axis 


J  (arg  S+  -  arg  S' 


(222) 
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APPENDIX  G 


THE  FUNCTION 


0  ejm  d0 


This  is  the  integration  with  respect  to  ip  for  the  phase-quadrature 
components  (Equations  (47)  and  (50)). 

First,  let  cp  =  0  or  rr  .  Then 


jm  cost/ 


'  d  0  - 


f  sin  2  0  eir 

‘'o 

/ff/2  s-  n  /*3 77  2 

sin  2  0  eimcos^  d0  +  /  sin  20  eimcos^  d0  +  /  sin  2  0  eim  cos  'A  d0 

u  *4/2  r 

/2?7 

sin  2  0  eimcos^  d0 
jff/2 

/^/2  .  r  ~  f  * 

sin  2  0  eJm  cos"  d0  +  /  sin  2  0  eim  cos  'A  d0  -  /  sin  2  0  eJm  cos  'A  d0 

u  *V/2  Jv/2 


17  /  2 


-  f  sin  2  eJm  cos  P  d0  =  0 

*/0 


since  to  eaeh  value  sin  2  a  e,m  L0S  u  in  the  first  quadrant  is  associated  the 
value  -  sin  2  0  c|,n  los  H'  jn  the  fourth  quadrant,  and  to  each  value 
-  sin  2  6  e|m  cos  ^  in  the  second  quadrant,  is  associated  the  value 
sin  2  6  eim  cos  *A  in  the  third  quadrant. 

Even  faetors  of  </  in  the  integrand  will  not  ehange  the  result. 

Now  let  <p  =  ±  7r  / 2.  Then  the  integral  is  again  zero,  exeept  that  now 
it  is  beeause  values  in  the  first  quadrant  cancel  with  corresponding  values 
in  the  third  quadrant,  whereas  values  in  the  second  quadrant  caneel  with 
corresponding  values  in  the  fourth  quadrant. 
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Hence,  the  far-field  quadrature  terms  become  zero  when  the  observa¬ 
tion  point  is  located  in  one  of  the  two  principal  planes  <p  =  0  or  <p  -  rr/2. 
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APPENDIX  H 


FLAT  CIRCULAR  PLATE  ILLUMINATED  BY  PLANE, 


CIRCULARLY  POLARIZED  WAVE 


This  case  may  be  considered  as  the  limit  of  the  paraboloid  as  the  focal 
length  is  extended  to  infinity  whereas  the  aperture  remains  constant.  As  a 
matter  of  fact,  it  is  often  substituted  for  the  paraboloid  in  the  literature,  * 
but  treated  without  regard  to  polarization. 

The  far-field  equation  is 


(223) 


The  vector  n  is  simply  (see  Figure  72) 


A 


n 


A 

z 


(224) 


A 

The  vector  a  is  still  (see  Equation  (10)) 


sin  0  cos  (f> 


A 

a 


sin  0  sin  0 


cos  0 


(225) 


I 


See,  for  instance,  L.  Thourel,  MThe  Antenna,  M  pp.  260-2G2,  John  Wiley 
and  Sons,  Inc.  New  York,  1960. 
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The  vector  p  ("center"  of  scattering  is  taken  at  center  of  the  plate)  is 


p  p  l  x  cos  tp  +  y  sin  (p 


(226) 


so  that 


jkp-aj  j  ^  p  sin  0  cos  (ip  —  <,’>) 


(227) 


The  differential  area  cLA  is  simply 

d  A  p  d  p  d  (Jj 

The  primary  pattern,  not  necessarily  uniform,  is  of  the  form 


(228) 


II;  h  (P) 


COS  COt 

sin  cot 

0 


so  that  the  eross-product  (n  x  H  )  becomes 


Cofactors 

of 


A. 

A 

A 

X 

y 

z 

—  sin  cot 

Co  factors  . 

of  h  W 

0 

0 

] 

ll  (/>) 

COS  (  >t 

COS  (Ot 

sin  cot 

0 

0 

A 

A 

the  cross-product  of  (n  x 

H  )  and 

becomes 

A 

s\ 

A 

X 

y 

z 

COS  (Ot 

h(P) 

—  sin  <ot 

cos  cot 

0 

h  ( p ) 

sin  cot 

—  sin  0  cos  - 

■  sin  0  sin  <f> 

cos  0 

(sin  0  sin  c/j)  sin  cot  + 

(229) 


(230) 


(231) 
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For  small  values  of  6  ,  the  expression  for  the  far-field  becomes 


(Equation  (77)) 


He  OC  X  COS  (rjt 

^  P 


A 


(k  p  sin  0)  h  (p)  p  d  p 


f  V  s 


f* 

in  on  /  J0  (kp  sin  0)  h  (p)p  d, 


(232) 


Thus,  it  is  seen  that  the  x-  and  y-components  are  equal  in  magnitude,  are 
in  phase-quadrature,  and  are  independent  of  <p  . 


For  large  values  of  0  ,  the  coordinate  transformation  of  Appendix  E 
can  be  used.  Thus,  by  Equation  (211), 


Jq  (k p  sin  0)  h(p)  p  dp 


+  y 


sin  (c  jI  —  c6)  [cos  6\ 


JQ  (kp  sin  0)  h  (p)pdp 


(233) 


Hence,  except  for  the  factor  cos  6  ,  the  far-field  is  everywhere 
circularly  polarized  if  the  feed  is  circularly  polarized.  The  40-location 
of  the  observation  point  P  only  affects  the  phase  of  the  circularly  polarized 
wave. 


1 
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Figure  72.  Geometry  of  the  Flat  Circular  Plate 
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APPENDIX  I 


ANALYTICAL  PROPERTIES  OF  THE  HYPERBOLOID 

Consider  the  locus  of  any  point  Q  (Figure  73),  such  that 

S  p  +  2<i  (234) 

where  a  is  a  constant.  If  a  perfectly  conducting  body  of  revolution  is  gen¬ 
erated  from  such  a  curve,  then  the  phase  of  a  ray  originating  from  a  point 
source  located  at  the  origin  will  have  the  same  phase,  at  point  Q  ,  as  if  it 
were  originating  from  point  C  —  within  a  constant  phase  corresponding  to  2 


i 


Figure  73.  Geometry  of  Hyperboloid 
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Let  the  distance  OC  be  2ea.  Then,  by  the  law  of  cosines 


P 


4e2a2 


1  ^ 

4s  e  a  cos  0)  1 


by  Equation  (234).  This  equation  can  be  solved  for  s  ,  resulting  in 


a  (1  -  e2) 

1  —  e  cos  0 


(235) 


When 


/. 


a  ( l  i  e)  (1  —  e) 
(1  -  e) 


a  f  ac 


which  is  the  location  of  the  apex. 

The  quantity  p  can  be  evaluated  from  Equations  (234)  and  (235)  as 


a  s  —  2  a 


ey 

a  (1  —  e~)  —  2  a  +  2  a  e  cos  (i 
1  —  e  cos  0 


or 

a  [  2  e  cos  0  —  U  +  e~)  1 

P  '  1  -  e  cos  0  (236) 

In  Cartesian  coordinates,  with  Q  in  the  x-x  plane,  the  components 
of  p  are 


px  s  sin  0 
p z  —  |  2  a  e  —  s  cos  0] 


i 
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For  any  azimuthal  angle  <p  ,  the  expression  for  p  is 


COS  <5 

—  sin  ci 

0 

s  sin  0 

s  sin  6  cos  <5 

p 

sin  (£> 

cos  6 

0 

0 

- 

s  sin  0  sin 

0 

0 

1 

s  cos  0  —  2 a  c 

s  cos  0  —  2  a  e 

But  s  depends  on  6  ,  and  it  is  best  to  express  p  in  terms  of  0  and  cp 
(and  constants)  only,  that  is,  by  Equation  (235) 

9 

(1  —  e **)  sin  0  cos 
(1  —  c“)  sin  0  sin  6 

9 

( i  -1-  e.)  cos  0  —  2  e 

(237) 

A. 

The  unit  vector,  normal  to  the  surface  at  point  Q  ,  n  ,  can  be 
obtained,  in  polar  form,  by  taking  the  gradient  of  Equation  (235),  and 
dividing  through  by  its  magnitude,  that  is 


1  - 


IS  0 


A 

n 


Vf 

!\f 


(« 


ILL 

d  s 


A  1  <9f 


a°  S  do 


fl 


where 


f  s  (1  —  e  cos  0)  a  (1  — 


It  becomes 


A 

n 


[as  (1  —  e  cos  0)  4  e  sin  0  I 
[(1  —  e  c  os  0)**  4  (e  sin  0)~  \ 
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But,  if  Q  is  in  the  x-z  plane, 


sin  0 

cos  0 

A 

0 

A 

as 

>  ap  = 

0 

cos  0 

—  sin  0 

so  that,  in  that  plane,  the  unit  vector  n  is 


A 

n 


(1  4  e  —  2e  cos  0) 


sin  0 
0 

cos  0  —  e 


Hence,  for  any  azimuthal  angle  <p  ,  the  expression  for  n  becomes 


Cos 

—  sin  <f> 

0 

i 

sin  0 

sin  4> 

cos 

0 

1 

0 

(1  +  e~  —  2e  cos  0) 

0 

0 

1 

(cos  0  —  e) 

(1  f  e~  —  2e  cos  0) 


sin  0  cos  (£> 
sin  0  sin 
(cos  0  —  e) 


(238) 


It  can  be  shown  that  the  angle  of  incidence  at  point  Q  ,  namely 

“1a  a 

cos  (-a  •  n)  ,  is  equal  to  the  angle  of  reflection,  where  the  reflection 

S  -1  —  *  . 
seems  to  be  originating  from  point  C  ,  namely  angle  cos  (p  •  n/p): 

this  property  together  with  the  property  of  Equation  (234)  constitute  the 
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justification  for  the  choice  of  the  hyperboloid.  Consider  Equation  (238)  and 
a  general  expression  for  a^ 


cos  sin  0 
sin  </>  sin  0 
cos  0 


The  dot  product  becomes 


A 

n 


e  cos  0  —  1 


( 1  +  e~  —  2 e  cos  0) 

By  Equations  (23G),  (237),  and  (238) 

(1  —  e  cos  0)  1 


A  A 

P  ‘  n 


a  [2e  cos  0  —  (1  f  e  )]  (1  f  c"  -  2e  cos  0) 


(239) 


(1  —  e  cos  0) 


9 

(1  —  e~)  sin  0  cos  cfc 

sin  0  cos  0 

(1  —  o^)  sin  0  sin 

sin  0  sin  <f> 

(1  +  e~)  cos  0  —  2e 

(cos  0  —  e) 

(1  +  e2  -  2e  ms  0)  2 


( 1  —  e^)  sin"*  0  +1(1  1  e*9  cos  0  —  2e1  (cos  0  —  c) 

[2  c  cos  0  —  ( 1  f  e~)  i 


The  numerator  can  be  transformed  and  factored  as  follows: 

1(1  —  e“)  sin"*  0  +■  (1  4  e“)  cos^  0]  —  2e  cos  0  —  (1  v  ei  e  cos  0  t  2e~ 
[  1  4  e ^  —  2  e"*  (1  —  cos“  0)1  —  2e  cos  0  —  (1  +  e~)  e  cos  0  +  2e~ 

(1  +  el  +  2e^  cos^  0  —  2e  cos  0  —  (1  +  e“)  e  cos  0 

(e  cos  0—1)  [2e  cos  0  —  (1  +  e2)] 
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Consequently,  the  expression  for  (p  •  n  )/p  becomes 

e  cos  0  —  1 

9  Vo 

(1  f  e~  —  2e  cos  0) 

A  A 

that  is,  it  is  equal  to  -a  •  n  .  Thus,  in  Figure  73,  angles  \  j  and  v.> 

s 

are  equal. 

Finally,  an  expression  for  the  differential  area  cLA  must  be  derived. 
In  polar  coordinates  (Figure  73),  the  general  expression  for  dA  is 

<1  \  l  (s  d<9)2  4  (ds)2],/l(s  sin  0  dci  ) 

1 

2 

s  sin  0(10  (kj 

By  Equation  (235).  the  derivative  ds/dfl  is  evaluated  as 

d  s  —  a  e  ( 1  —  e  _)  s  i  n  0 

d  0  /  j  /i\  2 

U  —  e  cos  v) 


(  pon  substituting  expressions  for  s  and  for  ds/d0  ,  we  obtain 


d  \ 


(I 


e2) 


1  —  c cos  0 


—  a  e  ( 1  —  e“)  sin  0 


(1  - 


or 


(1  -  e“) 


1  - 


sin  0  (I  0  do 


2  j  2 

- - y  (J  *  e“  —  2c  cos  0)  sin  0  d(i  do 

( 1  —  e  cos  Or 


(240) 
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APPENDIX  J 

PRIMARY  ILLUMINATION  FUNCTION  h(0)  FOR  CASSEGRAIN  SYSTEM 

The  amplitude  pattern  has  been  measured  experimentally  as  a  function 

of  0  (Figure  74).  If  the  maximum  subtended  angle  6,0  ,  is  to  be  30 

o 

degrees,  an  analytic  expression  can  be  obtained  over  the  interval.  By  trial 

2 

and  error,  it  is  found  that  the  function  cos  (26)  provides  an  adequate  fit, 
as  shown  in  the  following  table: 


6 

db 

(experimental) 

db 

(cos^  26) 

o 

0 

(J 

0 

10° 

1 

1.  07 

20° 

5 

4.  63 

CO 

o 

O 

12.  5 

12 

The  illumination  function  h(0)  will  be  multiplied  by  l/s  ,  as  given 
by  Equation  (107),  so  that  the  common  amplitude  factor  is  (within  a  constant) 

(e  cos  #  —  1)  cos^  26 
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Figure  74.  Primary  Illumination  Pattern 
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APPENDIX  K* 


COMPUTER  PROGRAMS 

EXAMPLE  OF  FOCAL  FEED  SYSTEM,  FAR-FIELD,  EQUIBEAM  AND 
10  db  TAPER  ILLUMINATION 

Program  "ANFPOLN”  (Figure  75)  is  used  to  compute  the  plots  of 
Figures  9,  26  and  43.  It  computes  the  field  components  at  various  locations 
in  the  far-field  of  a  focal-plane  antenna  with  equibeam  distribution  of  the  feed 
horn,  and  a  cos^*  (7r  -  rj )  taper.  The  notation  used  in  the  program  follows, 
as  far  as  practicable,  that  used  in  the  equations  from  which  it  was  developed. 
Although  this  particular  program  represents  a  far-ficld  case,  very  good  near¬ 
field  results  may  be  calculated  merely  by  changing  the  value  of  range,  e.  g.  , 
the  field  at  a  range  of  600  ft.  could  be  obtained  by  substituting  the  statement 
RANGE  =  600.  for  the  one  shown. 

The  operation  of  the  program  is  as  follows.  A  value  of  elevation  angle 
0  is  read  from  a  data  card.  Azimuth  angle,  q  ,  is  initially  set  to  zero. 

For  this  set  of  field-point  coordinates,  the  individual  contributions  to  the 
field  vector  components  are  computed  for  each  incremental  area  of  the 
reflector  (in  this  case.  1  degree  in  ^  and  77)  and  a  Simpson!s  rule  numeri¬ 
cal  integration  is  performed  over  the  surface  of  the  reflector.  After  the 
computations  are  completed  and  the  results  printed,  0  is  incremented  by  15 
degrees  and  the  computations  are  repeated  for  the  new  field  point.  Upon 
completion  of  the  computations  for  q  =  180  degrees,  the  program  is  ready 
to  read  a  new  value  of  0 . 

The  Simpson  Ts  rule  integration  is  performed  by  subroutine  SRI,  a 
listing  of  which  is  included. 

^Written  by  J.  L.  Pearlman 
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nn  ooooooon-*:»cx>cp 


1 


JOB, ANFPOLN,PEARLMAN  JL , 805 2 . , 08 5 • 0 1 5 , E 4 C2 , S 
TYPE,C0MP I LGO, FORTRAN 
L I S 1  ADDIO.PLBDL757 
SUBTYPE,  FORTRAN 

NEAR  FIELD  POLARIZATION  CALCULATIONS  FOR  J.  ALTMAN  3/22/66 
PER  MEMO  3/11/66 


******  NOTE 


EXPONENTIAL  TERM  MODIFIED  4/7/66 


********** 


MODIFIED  FOR  EQUIREAM  DISTRIBUTION  9/13/66 

DIMENSION  X  Y  IRE  <631,  XY1IE  (63),  X  2RE  (63),  X2IF  (63),  X3RE  (63), 
1X3IF  (63),  XY4RE  (63),  XY4IE  (63),  X5RE  (63),  X5IE  (63), 

2Y2RE  (63),  Y2IE  (63),  Y3RE  (63),  Y3 I E  (63),  Y5RE  (63),  Y5IF  (63), 
3XY1RP  (361),  XY1IP  (361),  X2RP  (361),  X2IP  (361),  X  3RP  1361), 

4  X  3 1 P  (361),  XY4RP  (361),  X Y4 1 P  (361),  X5RP  (361),  X5IP  (361), 

5Y2RP  (361),  Y2IP  (361).  Y3RP  (361),  Y3IP  (361),  Y6RP  (361), 

6Y51P  (361),  PHANG  (25),  Y2R  (25),  Y2I  (25),  Y3R  (26),  Y3  I  (25), 

7  Y5R  (25),  Y  5 1  (25),  ELL  (25),  FM AG  (25),  DEV  (25),  HETA  (63) 

P  I=4.*ATAN(  1. I 
F  3 1  2  •  5 


**  NOTE  -  RANGE 
RANGE3 1 . E  7 


INCREASED  INTO  FAR  FIELD  4/16/66 


C0NSTl=-24./9.22*Pl 
R  SO*  R  ANGE*  *  2 
FSQ=F**2 

TW0FR-2,*F*RANGE 
C2A=.5*F/RANGF 
C  2  8  =  4 • *C  2 A 
D  RC  3  P  (  /  1 8C  • 
ETA0=11R.07*0RC 
N  TH  =  0 

20  READ  10 , THANG 
N I H= 1 +NTH 
THFT  A=THANG*DRC 
S  TH=  S I N (  THETA) 

CTH=COS( THETA) 

PRINT  11, THANG 
DO  23  KC *  1 ,  13 
FKC=KC-1 

PHANG(KC)=15.*FKC 
PHI =PHANG( KC )*DRC 
S  PH  I »  S I N ( PHI  ) 

CPH( *COS( PHI ) 

ST  SP  =  S  TH*  SP  H I 
STCP=STH*CPHI 
no  22  KB* 1 , 361 
F  K  B  =  K  B- 1 
P  S I -DRC*  FK  B 
S  P  =  S  I  N  ( PSI  ) 

CP=COS(PSI  ) 

S°2- SP 

CP2*CP**2 

S2PH=SP*CP 

CONST3=STH*COS( PSI -PHI ) 
DO  21  K A* 1 18, 180 
KAA=KA-1 17 


Figure  75.  Program  ANFPOLN 
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o  o 


F  K  A  =  K  A 
F  T  A  =  DRC*FK  A 
E02=ETA  n. 

S  E=  S I N ( ETA  ) 

CE=COS( ETA  ) 

S  E  0  2  *  S IN{ E  0  2  ) 

C  EU2  =  CG  S  (  EO 2  ) 

CTtU2=CE02/SE02 
C3AR=1./SEU2**2 
C  3 A  =C  3  AR** 2 

H  (ETA)  AS  FITTED  TO  MEASURED  CURVE  4/26/66 
IF  ( (NTH.fcQ.l ) . ANU. (KC.EO. 1 ) .AND. (K0.EQ. 1 ) I 
lHfcTA(KAA)=COS(P I-E  T  A  )  ** (  1  • 5) /C3AR 
Ht=HET A( KAA ) 


GF=HF*SE/SE02**5 

C30=.5*C2A*C3A 

CONSTa=CQNST 1*( SORT ( R  S  O*  F  $U*C  3  A- T  *GF  R*C  3  A  ft* ( St  >CUNS T 3  +  C T H*C E )  ) 
1  ♦  F  *C 3AR ) 

AP=-SF02*GF 

BP=C. 

FP=-CEU2*GE 
C2C=C2B*CTE02 
CC4  =CU  S ( CONST  4 ) 

S  C  4  =  S  I  N  ( C0NST4 ) 

X  Y  1 R  £ (KAA) =AP*CC4 
XYIIE(KAA)  =  A  P*  S C 4 
T  A= ( C2C*SP-STSP )  *  E  P 
X3T=TA*SP 
X  3RE ( K AA ) =X3T*CC4 
X  3  I F (KAA) =  X  3  T  *  SC  4 
X  ST=-TA*CP 
XSRfc (KAA) =X5T*CC4 
X  5  IF (KAA )=X5T*SC4 
T  8  =  ( C2C*CP-STCP )  *E  P 

Y  3  T  =  Trt  +  CP 

V3RF ( KAA ) =  Y  3  T  *C  C  4 
Y3I t ( K AA ) =  Y  3  T  *  SC  4 

Y  5  T  =  - T6+  SP 

Y  5  R  f  (KAA) =YST*CC4 
YS  I  F ( K AA ) =  Y  5  T*  SC4 


21 

CUM  INUE 

CALL 

SRI 

( X  Y1RE 

,  DRC  »  K A  A 

,XY1RP(K0)  ) 

CALL 

SR  I 

(  X  Y  1  I  E 

, DRC , KAA 

, XY1 1 P<  KB  )  ) 

CALL 

SRI 

( X3RF, 

DRC , KAA, 

X3RP ( KB) ) 

CALL 

SR  I 

( X31E, 

DRC  »  K  AA  * 

X 3 1 P ( KH  )  ) 

cal  l 

SRI 

( X5Rt  * 

DRC. KAA, 

X5RP ( KB  )  ) 

CALL 

SR  I 

( X5IE, 

DRC  .  K  A  A  . 

X5 I P<  KB  )  ) 

C  AL  L 

SRI 

( Y  3  RE , 

DRC, KAA, 

Y  3 R P ( KB ) ) 

CALL 

SK  I 

<  Y3IE , 

DRC , KAA , 

Y  3  I  P  (  K  B  )  ) 

call 

SR  I 

( Y5RE , 

DRC , KAA , 

Y  5  R  P  (  K  B  )  ) 

CALL 

SRI 

<  Y5I  t, 

DRC , KA A  , 

Y5IP(Kft) 1 

22 

CUNT  I.NUE 

CAL  L 

SRI 

( X Y1RP 

»  OR  C  •  36  1 

*  X  Y  1R  ) 

CALL 

SRI 

( XY1IP 

,ORC , 361 

.  X  Y  1  I  ) 

CAL  L 

SRI 

( X3RP, 

DRC ,  361  , 

X3K  ) 

CALL 

SR  I 

( X3IP. 

DRC , 361 , 

X  3  I  ) 

(  ali 

S  *  I 

( X  6RP  , 

URC , jol* 

X5R  ) 

CALL 

SR  1 

(  X  b I P, 

DRC. 361, 

X  5  I  ) 

Figure  75  (Continued).  Program  ANFPOLN 
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1 


CALL  SRI  (Y3KP«0RC*361*Y3R(KC) I 
CALL  SKI  ( YJIP,DKC.361,Y3I(KCI ) 

CALL  SRI  (Y5RP*DRC*36l*Y5R(KC) ) 

CALL  SRI  ( Y5IP»DKC'361tY5I(KC) ) 

XR=XYlR+X3R-X5I 
XI=XY1I+X3I*X5R 
YR=-XY1 I  -  Y  3  I  ( KC I +  Y5R I KC ) 

Y  I  =  X  Y  1  R  ♦  Y  3R  (  KC ) +  Y5 I ( KC ) 

SPLU  =  .5*SQRT((XR-Yl)v*2MXI  +  YR)**2) 

SMIN*. 5* SORT (  ( XR  +  Y  I  )  *  *  2  ♦  ( X  I  - YR  )  *  *  2  ) 

ELL ( KC ) =1G.*AL0G10(  (  ( SPLU+SM I N  )  / ( SPLU-SMI N)  )*+2) 

X RN  =  - X 5  I 
X  IN  =  X5R 

Y  R  N=  —X  Y 1 I-Y3I (KC ) 

YIN=XY1R*Y3R( KC  ) 

XHAG  =  SQKT ( XRN**2+X I N**  2  I 
YMAG=SDRT(YRN**2*YIN**2I 
FMAG(KC)=S0KT( XMAG**2+YMAG**2 ) 

0EV(KC)=ATA N( xmag/ymagi /dkc 

PRINT  lCtPHANfi(KCI  • XY1R , XY1 I • X3R, X3 I . X5K, X5  I 

23  CONTINUE 
PRINT  1 2  *  T  HANG 
DO  24  L*lt l 3 

PRINT  10 • PH ANG(L ) , Y3R ( L  K Y3 I ( L  ) . Y5R ( U • Y5  I  ( L I 

24  CONTINUE 

PRINT  1 4 «  THANG  * 

00  25  L  =  1  *  1  3 

PRINT  1  }»PHANO(L I • E LL ( L) • FMAG ( L ) t Ot V ( U 

25  CONTINUE 
GO  TO  2C 
RF  TURN 

1C  FORMAT  (F1G.1.6E15.5) 

11  FORMAT  I ili/llOXfPAGb  1  OF  3$// 

ilLXSFAR  FILL!)  CALCULATIONS  FOR  J.  ALTMAN  4/16/66$/ 
ciCxfbXACT  NEAR  FIFLO  EXPRESSION  *i  I  TH  VERY  LARGE  RANGE  VALUES/ 

3UXSIILUMINAT  ION  FUNCTION  FITTED  TO  MEASURED  DATA  2/26/66$/ 

CIOXSMUUIFIFO  FOR  FQUI3EAM  D I S  T  K I  BUT  I  ON  9/13/66$/ 

4 1 C  X  $  FUR  THETA=$F6.  I  ,  $  Of.GK  E  E  Si  / 

57x$PHl $5X4RF  I  X  1  , Y  1  )  $ 6 X  $  I  V  l  X  1  t  Y  1  )  $  5  X  $R  F  (X3)$8X$lM  1X31$ 

6  8  X  SR  F  1X5) $ 8 X  $  I M  ( X5  I  $  I 

12  FUR.MAT  (  $1  S/1  10X$PAGF  ?  OF  3$// 

1 iJxSFAK  FIELD  CALCULATIONS  FOR  J •  ALTMAN  4/16/66$/ 

2 l C  X  $  E  X  AC  T  NEAR  FIELD  EXPRESSION  WITH  VtRY  L  ARG t  RANGE  VALUES/ 

31 wXSILLUMINATIUN  FUNCTION  FI T  T  F  0  TO  M  E  A  SU  RF  0  DATA  4/26/66$/ 

Cl 3XSMOOIF IFO  FOR  FQUIBEAM  DISTRIBUTION  9/13/66$/ 

410XSFOR  THETA=$F6. 1 , $  DTGRFESi/ 

57X SPHI $5X$RF  (  Y  3 ) $  8X  S I M  (Y3)$8X$RF  IY5)$8X$IM  I Y5 ) S  ) 

13  FURMAT  I  Flu,  1  ♦  3F20.4 ) 

14  FORMAT  (  $l$/ilOXSPA(,E  3  OF  3$// 

UGXiFAR  F  I  ELI)  CALCULATIONS  FOR  J.  ALTMAN  4/16/66$/ 

21CXSEXACT  NEAR  FIELD  EXPRESSION  WITH  VERY  LARGE  RANGE  VALUE  $/ 

3 1  OX  $  I LLUM I NAT  ION  FUNCTION  FITTED  TO  MEASURED  DATA  2/26/66$/ 

ClCXiMOJIFIFD  FOR  E  QU I  BE AM  DISTRIBUTION  9/13/o6$/ 

^rlOXSFOR  THt  T  A  =  SF6.  1  *  $  DEGREESS/  ' 

57X$PHl$3X$cl LIPTIC  ITY  (OB. I $11 XSMAGNI TUDE  $4 X $DF V  I  A T  I  ON  I  DEG. ) $ ) 

PNO 

SUMTYPt  »  0  A  T  A 

•  ? 


Figure  75  (Continued).  Program  ANFPOLN 
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o  o  o  o  n 


T  SUBTYPF, FORT  RAN, LSTRAP,NRL.SEQ 

SUBROUTINE  SRI  (Y,DX,N,Y1I 

NUMERICALLY  INTFGRATCS  Y  OX  BY  SIMPSON'S  RULE 

Y  =  ARR A  Y  TO  BE  INTEGRATFD 
DX  =  X  INCRFMfcNT 

N  =  NO •  OF  POINTS  IN  Y.  IF  N  NOT  ODD,  SUBROUTINE  CAUSES  ABfcfJJ 
WITH  PRINTOUT  OF  CURRENT  VALUE  OF  N. 

DIMENSION  YIN) 

IF  ( (MODIN* 2) ).Nt. 1 >  GO  TO  21 
SUM=  Y  (  lDY(N)f4,*Y(2) 

N  A  =  N - 1 

DU  20  K  =  3  *  N  A  »  2 
SUM=SUM*2. *Y ( K)*4.*Y(K+1) 

20  CONTINUE 

Y  I  =SUM*DX/3. 

RETURN 

21  PRINT  10, N 

sru^ 

1  vj  FORMAT  (SCN  =  fUS,S,  *hlCH  IS  AN  iLLtGAL  VALUE.  S> 

FND 


Figure  75  (Concluded).  Program  ANFPOLN 
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EXAMPLE  OF  CASSEGRAIN  SYSTEM:  FIELD  ON  PARABOLOID, 
DIPOLE- LIKE  ILLUMINATION 


Program  "ACFSHB"  (Figure  76)  computes  the  vector  components  of  the 
field  scattered  onto  a  parabolic  reflector  from  a  hyperbolic  subreflector  and 
stores  the  results  on  tape  for  further  processing. 

The  operation  of  this  program  is  similar  to  that  of  the  previous  pro¬ 
gram  except  that  the  coordinates  of  the  field  points  are  determined  by  the 
geometry  of  the  parabolic  reflector.  The  points  on  the  reflector  for  which 
the  field  is  calculated  are  separated  by  15  degrees  in  1/  and  1.9375  degrees 
in  r\  (this  value  being  the  closest  to  two  degrees  which  is  an  exact  odd 
divisor  of  the  angle  subtended  by  the  parabola).  The  components  of  the 
total  field  at  each  point  on  the  parabola  are  integrated,  using  Simpson's 
rule,  over  the  surface  of  the  hyperbolic  subreflector  using  an  incremental 
area  of  the  one  degree  in  0  and  0  . 

A  dipole-like  distribution  at  the  feed  horn  is  assumed. 

Again,  subroutine  SRI  is  used  to  perform  the  numerical  integration. 
EXAMPLE  OF  CASSEGRAIN  SYSTEM:  FAR- FIELD 

Program  "APCFA"  (Figure  77)  computes  the  components  of  the  field 
vector  in  the  far-field  of  a  Cassegranian-feed  antenna. 

The  operation  is  similar  to  that  of  "ANFPOLN"  except  that  the  field 
components  on  the  surface  of  the  reflector  are  not  computed  in  this  program, 
but  are  read  in  from  a  binary  tape  generated  by  another  program  (such  as 
"ASFSHB").  Also,  the  incremental  area  on  the  surface  of  the  reflector  was 
1.9375  degrees  in  V  and  15  degrees  in  . 

A  recently  added  facility  for  handling  complex  numbers  was  used  in 
this  program. 
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R 

S 

RL  I  R 
T 

B3R 

T 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 


JOB»ACFSHR,PEARLMAN  JL*  RO  52  ♦  ,085,015, FhO 2, S 

TYPE, COMP iLGOt FORTRAN 

ADOIO,  PLB01674 

SUBTYPE, F (00 

100, TAPE ,  , , ,  ,  , S  AVE 

END 

SUBTYPE, FUKTKAN 

CALCULATION  OF  COMPONENTS  OF  H(S/P)  AS  PER  MEMO  BY  J.  ALTMAN 
4/29/66 


ORIGINAL  PROGRAM  WRITTEN  5/6/66  (J.  L.  P.  ) 


VtRSION  TO  WR  (  T  E  BINARY  TAPF  FUR  C ASSFGRAN ( AN  FEED  PROGRAM  7/19/66 

‘  ** ** **** **** 


01  POLE  -  LIKE  0 ( STR ( BUT  I ON 


OIMFNSIUN  PSANG  (26),  PS(  (25),  SPS  (251,  CPS  (25), 

1 X  Y  7  (12),  PHI  (361),  SPH  (361),  CPH  (361),  XlPR  (361),  X1PI  (361), 
2  X  2  PR  (361),  X2PI  (341),  Y i PR  (361),  Y1PI  (361),  Y  2  P  R  (361), 

3Y2PI  (361),  7 1 P  R  (361),  Z1PI  (  361  ),  72PR  (361),  Z2P(  (3ol), 

hTHETA  (31),  STH  (31),  CTH  (31),  TTH(31),  CSOTT  (31),  X1TR  (31), 

5  X  1 T I  (31),  X  2  TR  (31),  X2T(  (31),  Y 1 T  R  (31),  Y 1 T (  (31),  Y  2  T  R  (31), 

6  Y  2  T (  (31),  7 1 T  R  (31),  Z  1  T  I  (31),  72TR  (31),  Z?T(  (31) 

P I =4 • *  AT  AN (  1.  ) 

T  WOP  I =  2 ♦ *P ( 

F  CON  =  12.5 

CONST 1-12.*T WO PI/9.22 
ACON  =  1  •  G  6  3 
f  XC  =  2.45 
hXSO=EXC **2 
OPEXS-l.^EXSO 
GMtxS-l*-ExSO 
ADMFX=OMEXS*ACON 
ORC  =  P ( /  18C  . 

DR  C  2= 31.*0KC/ln. 

T  F  A=2.  *FC0.\1*AC  IN 
T  WE  =  2 • *EXC 
ASQ= AC0N**2 
A  (J  ME  XS  =  AOMt  X**2 
00  3C  NE_=1 , 33 
FNE=33-NE 
F  T  A  =  P  I  -  F  NE  *■'  ORC  2 
E02  =  F  TA/2 • 

SF=SIN( ETA) 

C. E  =  C 0 S  (  ETA  I 
SF02  =  S  I IM(  EO?  ) 

SSEO?=SF 02**2 
F0M-FC0N/SSF02 
«  T  A= FUM**2 
RTC1 =TFA/SSEU2 
RTC2A=SE*0MEXS 
RTC3A=Ct*UPEXS 
k TC4  =  T  WE*CF 
A  1  X ]  TA  =  -E04*Sr 
A  1  7  1  T  =  -Ei)m*CF 
WR  I  Tb  (  3  )  HA 


Figure  76.  Program  ACFSHB 
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00  29  KC*1 ?  2 5 
IF  (NE.F0.1)  GU  TO  21 
GO  TO  22 

21  FKC=KC-1 

PSANGI KC ) =15,*FKC 
PSKKC)  =PS ANG  (  KC  )*DRC 
SPS(KC)*SIN(PSI  (KC)  ) 

CPS(KC)=COS(PS  I  (KC  )  ) 

22  SECP=SE*CPS(KC ) 

SESP=$E*SPS(KC) 

A1X1T  =  A1X1TA*CPS(KC  ) 

A1Y1T=A1X1TA*SPS{KC) 

00  28  K 8= 1 1  361 

IF  ( (NE.EQ. 1) . ANO. IKC.EO.l) )  GO  TO  23 
GO  TO  24 

23  F K H  =  KB  — 1 

PHI { KB ) =  FK  B ♦ORC 
SPH( KB )=SIN( PHI  {  KB  )  ) 

C  PH { KB) =  COS( PHI  ( KB )  ) 

24  SPHSO=SPH( KB ) **2 
CPHSQ=CPH(KB)**2 
SCPH=SPH(Kb)*CPH(KB) 

C  Af)  =  CO  S  (  P  S  I  (  KC  ) -PH  I  (KB  )  I 
RTC28=RTC2A*CA0 
CiMECP*OMExS*CPH(  KB  I 
AOMESP=AOMEX*SPH(KB> 

AOMFCP=ACON*OMtCP 
00  27  K  A= 1  *  3 1 

IF  (  (NF.EO.  1).  ANO.  (KC.EU,  1  )  .AiMD.  (KB.EU.l)  )  GO  TO  25 
GO  TO  26 

25  FKA=KA-1 

T  HE  T  A  { KA ) =FKA*ORC 
STH(KA)=SIN( THETA(KA) ) 

CTHl KA ) =CUS {  THt  T A{ KA )  I 
T  TH{ KA ) =2. ♦THE  T A(K A ) 

CSOTTIKA)»C(JS(TTH(KAI  )**2 
2<>  ECTH  =  EKC*CTH(KA) 

OMFCT=  l.-ECTH 

R=SORT( RTA+ASQ*( (2 • *FC T H-OPF X S ) / OMFC T I * *2 -R TC 1 /OMFC T 
1*1 RTC28*STH(KA) +KTC3A*(TH(KA)-PTC4) ) 

RR=  1. /ft 

A  1  X  P  =  KR ♦{ A1 X  IT* AOMECP*STH(KA ) /OMECT) 

A1YP-RR*!  A1Y1T+AU*ESP*STH(KA)  /UMhCT  ) 

A1ZP»RR*(  AlZlT+ACUN/OMFCT*(OPFXS*r.THf  KAI-TWEI  ) 
TECTH*TwE*CTH(KA) 

HTH=-UMFCT ♦CSQTT (KA ) 

HTHP  =  HTH/SOkT  {  0  P  F.  X  S-  T  fc  C  T  H  ) 

STHSQ=STH(KA)**2 

SCTH=STH(KA)*CTH{KA) 

bcot=hthp*sthso 

EFT  =  HTHP*SC  T  H 

A  =  HTHP*CTH(KA)*(EXC-CTH(KAM 

B=-bCDT*SPHSO 

C  =  —  BCO  T *  SC  PH 

0  =  BC  OT ♦CPHSO 

F=EFT*CPH(KB) 

F  =  tFT*SPH(KB) 


Figure  76  (Continued).  Program  ACFSHB 
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EXP=-C0NST1*(AUMEX/0MECT-F0M+R) 
DAT=A0MEXS/0MECT**3* <  GPEXS-TECTH ) *  STH ( K  A ) 
FXPR=COS( EXP ) 

FXPI=-$IN(EXP) 

A  P  B=  A  +  B 
DM A  =  D- A 


27 


D  AQR  =  RR  *DAT 

X 1T=DA0R*( A 1ZP*DMA+A1 YP*F ) 
X2T=-DAUR*<  AlZP*C+AlYP*t) 

Y 1 T  =  D  AOR* (A1XP  +  E-A1ZP*APB) 

Y2T  =  -D  A(JR*  (  A1ZP*C+A1XP*F  ) 

Z  1  T  =  D ADR*  (  A  1YP*C  — A  1  X  P*  0  M A  ) 

Z  2  T  =  D  ADR* ( A1 YP* A PB ♦ A  1  XP * C ) 
X1TR(KA)=X1T*EXPR 
X  1  T  I  (KA)=X1T*EXPI 
X  2 T  R ( K  A ) =  X2  T*E  XPR 
X  2  T  I  (KA )=X2T*EXPI 

Y  1  TR ( K  A ) =Y1 T*F  X  PR 
Y1TI(KA)=Y1T*EXPI 
Y?TR(KA)=Y2T*EXPR 

Y  2  T I ( KA ) = Y2T*FXP I 
7  1  T  R  (  K  A  )=Z1T*FXPR 
7  1  T  I  (KA)=Z1T*EXPI 
Z2TR(KA)=Z2T*EXPK 
7  2 T I (KA  )  =Z2T*FXPI 
CUNT INUE 

(  X1TR*  DRC  »  31  . X 1  PR (  Kb)  ) 
(  X  IT  I  *  ORC  *  h.  X1PI  (KB  M 
(X2TK.DRC.31.X2PR(KB) ) 


CALL 

call 

CALL 
CALL 
CALL 
(  ALL 

call 

call 

CALI 

call 

f  Al  L 
CALL 


SR  I 
SR  I 
SR  i 
SR  I 
SR  I 
SR  I 
SR  I 
SRI 
SR  I 
SR  I 
SR  I 
SR  I 


31 

( X  2 T I ♦ DRC  v 31  *  X2PI  (KB)  ) 
( Y1TR*0RC.31*Y1PR(KB» ) 
31 


( Y1TI  *  DRC , 
(Y2TR.ORC. 31, 
( Y  2  T I  .DRC *  31 , 


YIP  M  KB  )  ) 
Y^PR (KB)) 
Y2PI  (KB)  ) 
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CUNT INUE 
CALL  SRI 


CAL  L 
CALL 
CALL 
CALL 
CALL 
C  ALL 
C  A  l  L 
CA!  L 
LAI  L 
C  A  L  L 
C  Al  L 
ViR  I  TF 


SRI 
SRI 
SR  l 
SR  I 
SRI 
SR  I 
SR  I 
S  R  I 
SK  I 
SRI 
SR  I 
(  3 ) 


(  Z 1  TR  ♦ DRC  «  31 1  Z 1  PR ( KB )  ) 
( Z1TI *  ORC  t  31.Z1PI  (Kn)  ) 
( 7 2  T  K  *  ORC .  *1  .  Z2PR ( KB )  ) 
( Z  2  T I . ORC  *  31 ♦ Z2P I KB)  ) 

( X 1PR. ORC • 361 ,XYZ ( 1 ) ) 

( X1PI .DRC. 361,XYZ<2) ) 

<  X2PR. DRC  *  361  .  xY7 (  3 )  ) 

^6i  ♦ ; 

361  , 


( X2PI , 
( Y1PR, 
(YiPI, 


DRC  i 
DRC  , 
DRC  , 


,  XYZ ( 4  )  ) 
.  XY7 ( 5  )  ) 
361 , XYZ ( 6) ) 


(  Y2PR .DRC . 361 . XYZ<  7  )  ) 
(  Y2P I  . DRC . 361 , XYZ ( 3  )  ) 
(  Z 1  PR  *  OKC *36l*XY7 (4)  ) 


( 7  1PI  , 
( Z2PK, 
( Z2PI , 


ORC  , 
ORC, 
ORC  , 


36 
3  61* 


,  XYZ  (  10)  ) 
XYZ  (ID) 


36l  .XYZ  <12  )  ) 


P  S I ( KC  ) , XYZ 


20 

30 


CUNT  I  N'JF 
CUNT INUE 
LND  FILE 
REWIND  3 
RETURN 
f  NO 


Figure  76  (Concluded).  Program  ACFSHB 
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B 

8 

BLIB1 

T 

B3B 

ft 

T 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


c 


JOB,  APCF  A,  PEARL  HAN  J L , 80 52 . , 085 , 0 1 5 , E 4 02 , S 

type. comp  it  go, Fortran 

ADD  1 0*  PL  BD1 9^4 
SUBTYPE . f I QD 
100, TAPE 
REEL »  PLBD1 7  92 
END 

SUBTYPE, FORTRAN 

FAR  FIELD  POLARIZATION  CALCULATIONS  FOR  J.  ALTMAN 


C ASSEGRAN  IAN  FEFO  8/3/60 


COMPONENTS  OF  FIELD  SCATTERED  ONTO  PARABOLIC  RE  ELECTOR  FROM 
HYPEBOLIC  SUriREFLECTUR  READ  INTO  ARRAY  XYZ  FROM  BINARY  TAPE  (3) 
GENERATED  BY  PROGRAM  ACFSHB 


i  **  ** ** **** 


*****  NOTE  -  DELTA  ETA  =  2  DEG.  AND  DELTA  PS!  =  15  DEG.  ********** 


EXACT  EXPONENTIAL  TERM  USED 


DIMENSION  X  PR  fc  (  3  3)  ,  XP I E  (33),  XQRE  (33),  XQIE  (33),  Y  PR  E  (33), 
1YPIE  (33),  YORE  (33),  YQIE  (33),  SE  (33),  CE  (33),  SE02  (33), 

2CE02  (33),  SE02F  (33).  ECHK  (  3  ;> )  ,  XPRP  (25),  XPIP  (?5)t  XORP  (25), 
3X0IP  (25),  YPRP  (25),  YPIP  (25),  YQRP  (25),  YO ! P  (25),  SP  (25), 

4CP  (25),  PSI  (25),  PSCHK  (25),  PHANG  (13),  YPR  (13),  YPI  (13), 

5Y0R  (13),  YQ!  (13),  S^HI  (13),  CPHI  (ii),  PH!  (13),  ELL  (13), 

SUE V  (13),  F  MAG  (13),  X  Y  Z  (12,25,33) 

COMPLEX  EXPT,  Xi,  X2,  Yl,  Y2,  Zi,  12 ,  ZA1,  ZA2 •  ZB1,  /  82  ,  ZC1, 

1  Z  C  2 ♦  XPCT,  XOCT,  Y  PC  T ,  YQCT 
REWIND  3 
P  I  =<♦ .  *  aTan  (  1 .  > 

E  =  1 2 . 5 


***jJr********X***t**  ************** 


k ***************** 


*****  RANGE  EXTENDED  T(J  EAR  FIELD  a******************************* 
RANGE* 1 • E  7 


CONST  1 *-24. /V. 22  *P I 
R  SO  =  R  ANGF**2 
ESQ=F**2 

TW0ER=2.*E*RANGE 

C2A=.5*F/RANGE 

C?B=4.*C2A 

DKC  =  P I / 1 80 • 

DRC2=URC*31./16. 

DRC 15=15. *URC 
N  T  H  =  0 

20  READ  lO.THANG 


NTH=1+NTH 
THET  A=THANG*DKC 
S !H=SIN( THETA ) 
CTH=COS( THETA) 

PRINT  ll.THANC, 

DO  30  K C  =  1 ,  13 
IE  (NTH.EO.l) GO  TO  21 
GU  TO  22 


I 


Figure  77.  Program  APCFA 
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21  FKOKOl 
PHANG(KC)=15.*FKC 

PHI ( KC )  =  PHANG ( KC )  *  OR  C 
SPHI (KC)=SIN(PHl (KC) ) 

CPHI (KC)=COS(Ph I (KC ) ) 

22  STSP=STH*SPHI ( KC ) 

STCP=STH*CPHl (KC) 
no  29  K 8  =  1 t  2 5 

I F( ( NTH.EO. 1 ) .ANO. ( KC.tO. 1 ) )  GO  TO  23 
GO  TO  24 

23  F  K  R-K  8 - 1 

PSI ( K  3 ) =  FK  B*DRC i 5 
SP(KB) =SIN(PSI(K8) ) 

C  P ( K  B  )  =COS ( PSI(KB)  ) 

24  SP2=SP(KR)**2 
CP2=CP(KB)**2 

SP  2PH=  SP ( KB ) ♦C  P ( KR ) 

CONST3=STH*COS( PSI CKRI-PHI CKCI ) 

DO  2  H  K AA  =  1 « 33 

IF  (  (NTH.EO. 1). AND. (KC.tO. 1 ) . A  NO . ( KB . 1 0. 1 >)  GO  TO  25 
GO  TO  27 

25  FKA=33-KAA 

READ  (3)  ECHK(KAA) 

DO  26  K  P H  = 1 ,25 

READ  (3)  PSCHK(KPH),(XYZ(L,KPHfKAA),L=lfl2) 

2o  CONTINUE 

F T A=P I-FK A*DRC2 
E02=FTA/2. 

SE(KAA)=SIN(ETA) 

C  E ( K  A  A  )  =COS ( ETA ) 

St02( KAA) =S IN( tU2) 

SEO?M  KAA  )  =  SEU2  (  KAA  1**5 
C  c  nr ( KAA ) =COS ( EG  7 ) 

27  C  3AR  =  1  • / SEU  2 ( KAA 1**2 
C  3 A  =  C  5  AR  **  2 

Cr.NST4=C0NST  1*  (  SQR  T  (  RSQ«-PSu*C  j>  A-*TW0FH*C3AR*  (  SF  (  K  A  A  )  *CONS  T  3 
1 ^CTH*CE (KAA)  )  )  +F  *C  3  AR  ) 

CTtU7=CF02(KAA)/SE02(KAA) 

C  7C  =C2tt*CTE02 
CCs=C0S(C0NST4) 

SCa^SINICONST*) 

GE=SE( KAA ) /5E0  2MKAA ) 

SPctH=SP( KB ) *C  E02 (  KAA) 

CPCFH=CP(K3)*CE02(KAA) 

T A=C2C*SP(KB)-$TSP 
TH  =  C2C*CP(Ke>)-STCP 
E  XPT=CMPL  X ( CC4, SC4 ) *Gt 

X 1*CMPLX ( AYZ  U  t  KB*  KAA)  *XYZ (?  t KB.  KAAI  ) 
x2=CMPLX( X  Y  Z (  3, KB, KAA)  • X Y Z ( 4 , K 3 , K AA )  ) 
Yi=CMPLX(XYZ(5,KB,KAA),XYZ(6,KB,KAA)) 

Y2  =  CMPL  X ( XYZ (  7, KB, KAA)  #  X YZ ( R  t  KB»  KAA ) ) 

Z 1=CMPLX ( XYZ (9, KB, KAA) ,XYZ(10,KRfKAA) ) 

Z?=CMPLX(XYZ(11,K0, KAA )  , X YZ ( 1 2  •  KB » K AA )  ) 

/ Al=-Yl*SE02( KAA )-Z  2  *  S  P  C  E  H 
ZA2=-Y2*SE02( KAA )-Z1*SPCFh 
ZB1=X1*SE02( KAA ) +Z 1 *CPC  E  H 
ZB2=X2*SE02 (  K  A  A  )  +Z  2  *C  PC  F  H 
ZC1=X1*SPCEH-Y2*CPCEH 
ZCc=X2*SPCtH-Yl*CPCEH 
*PfT=-tXPT*(Z3WCl*TA) 

Figure  77  (Continued).  Program  APCFA 
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XOCT=-EXPT*(  ZB2+ZC2*TA  ) 
YPCT=EXPT*(ZA1*ZC2*TBI 
Y OC  T  =  E  XP T* ( ZA2+ZC1*TB) 
XPRF  ( K  A  A ) =  R  E  AL ( XPCT  ) 


XPIE(KAA)=A l MAG ( XPCT) 
XQRE (KAA) =REAL  (  XOCT  ) 
XUIE ( KAA) =  A  l  M  A  G  ( XOCT) 
YPK  E ( K  A A ) =REAL ( YPCT ) 
YPI E( KAA) =AIMAG(YPCT ) 
YQRE(KAA)=REAL  < YOCT ) 
Y01F(KAA)=AIMAG( YOCT) 


2  B 

CONTINUE 

CALL 

SRI 

CALL 

SR  I 

CALL 

SRI 

CALL 

SRI 

CALL 

SRI 

CALL 

SRI 

CALL 

SR  I 

CALL 

SRI 

29 

Cunt inue 

call 

SRI 

CALL 

SRI 

CALL 

SRI 

call 

SRI 

CALL 

SRI 

CALL 

SRI 

CALL 

SRI 

CALI 

SR  I 

X P KP  f  ORC 15*25 
X  P I P  »  ORC 1 5  t  2  5 
XORP  t  DRC 1  b »  25 
XOl  Pt  ORCiWb 
YPR P  f  DRC 15*25 
YP  I  P  *  DRC  15*25 
Y  OR  P  «  ORC  1  5  »  2  5 
YOlP,DRC15,2b 

X  K=  X  PR- X Q I 


> XPRP(KB) ) 
»  X P l P ( KB  )  ) 
. XORP ( KB )  ) 

>  XQ I P ( KB  )  ) 

>  YPR P ( K B  )  ) 
i Y P I P ( KB  )  ) 
, YORP( KB)  ) 
, Y  0  I P ( KB)  ) 

*  X PR  ) 

*  XPI  ) 

,  X  OR  ) 

*  XOl  ) 

t  YPR ( KC  )  ) 
* YPI  ( KC  )  ) 

*  Y Ok ( KC  )  ) 

*  YOI C  KC  )  ) 


X  I  =  X P I ♦ XOR 
YR=YOR(KC)-YPI (KC) 

Y l  =  YQ I  ( KC ) ♦  YPR ( KC) 

SPLU  =  .  5*  SORT  (  CXR-Yl  1**2  MX  I  *YR  )**2  I 
SMIN=.  5  ♦SORT  (  (  XK  +-Y  I  XI-YK)^^2) 

ELL(KC)=1C.*ALUGU(  (  (  SPLU+SMIN)  /  (  SPLU-SMl  N)  >**2) 
XkN=  XPR 
X  IN- XP  I 
YR  N  =  YQR ( KC  ) 

YIN=YOl (KC) 

XMAG=SJRT ( XRN**  2  +X I  N**  2 ) 

YMAG=SORT  (  YRN**2*Y  I  N**  2  ) 

FMAG(  KC  )  =SORT  (  XMAG**2*YMAG**2 ) 

OEV( KC  )  =  A  T  A  N  ( XMAG/YMAG) /URC 
PRINT  1C*PHANG(KC)  tXPRtXPI «  X  OR  «  X  0 1 
3L  CONTINUE 

IF  (NTH.EO.l)  CALL  FPRtE  (J) 

PRINT  1  2  *  T  HANG 
DO  31  L  =  1  *  13 

PRINT  lGtPMANGd  )  »  YPR(  l  )  ,  YPI  (L  )  •  YOR  (U  t  YOI  (  L) 

31  CONTINUE 
PRINT  1 3  *  TH  ANG 
DO  32  L  =  1  i  13 

PRINT  1  »  PH  ANG (  L  ),ELL(L),FMAG(L) 

32  Continue 
G(i  TO  2C 
RETURN 

10  FORMAT  ( F10.1 , 7  b 1 5  -  5 ) 
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11  FORMAT  (  SI  S//1  1CXSPAGE  1  OF  it// 

1 10XSFAR  FIELD  POLARIZATION  C AL CUL T I ONS t / 
210XSCASSEGRANIAN  F  Ft 0  ANT  t  NN A  t / 

3 1  OX  S  FQR  THE TA= tFfc.  1 f S  DEGREES*/ 

47XSPHI SRXSRF  (XP)tHXtlM  (XPJtHXtRF  (XQ)SRXSIM  <  XO )  S  I 

12  FORMAT  (  tlt//UCXtPAGt  2  OF  j>S// 

1 lCXtFAR  FIELD  POLARIZATION  C AL CUL T I  ON S S / 

2  1CX  tCASSEGRAN I  AN  FttO  ANTENNAS/ 

310XSF0R  THtTA=SFo. l.S  UFGREESS/ 

47XSPHI S5XSPE  CYPItRXtIM  (YPJS8XSRE  (Y0IS8XSIM  ( YQ  )  t  ) 

13  FORMAT  ( S1S//11CXSPAGF  3  OF  3t// 
lloXSFAR  FIELD  POLARIZATION  C AL CUL T I  ON S S / 
210XSCASSEGRANIAN  FEED  ANTENNAS/ 

31CXSFOR  THE TA=SF6. It S  DFGR  E  E  S  S / 
47XSPHlS3XSbLLlPTICITY  {  D  8  •  )  SI  1  XSMAGN  I  TUOt  S  ) 

14  FORMAT  (FlC.lt3F20.4l 
END 

T  bUHTYPt *  DAT  A 

.  2 


C.  5 
0.7 
1  .C 
1.  8 

1 . 7 
2 .  G 
2.0 


Figure  77  (Concluded).  Program  APCFA 
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